













PREFACE 


This research was conducted in cooperation with the NASA Langley Research 
Center (Fluid Mechanics Division - Theoretical Flow Physics Branch) and the Institute 
for Computational and Applied Mechanics (ICAM) of Old Dominion University during 
the period 1984 through 1990. The work on basic formulations and computational 
procedures was completed by the end of 1985 and was published as a progress 
report “Radiative Interactions in Transient Energy Transfer in Gaseous Systems,” 
NASA CR-1 76644, December 1985. 

The formulations and procedures developed in this study were applied to solve 
several realistic problems during 1985-1990. This resulted in various publications in 
forms of technical reports, technical papers (presented at national and international 
conferences), and journal articles. Selected articles are included in the appendices of 
this report to demonstrate the applications of the developed techniques in various 
fields involving radiative interactions in molecular gases. The procedures developed 
are being used at present to investigate radiative interactions in other challenging 
problems. 
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ABSTRACT 


Analyses and numerical procedures are presented to investigate the radiative 
interactions in various energy transfer processes in gaseous systems. Both gray and 
nongray radiative formulations for absorption and emission by molecular gases are 
presented. The gray gas formulations are based on the Planck mean absorption 
coefficient and the nongray formulations are based on the wide-band model 
correlations for molecular absorption. Various relations for the radiative flux and 
divergence of radiative flux are developed. These are useful for different flow 
conditions and physical problems. Specific plans for obtaining extensive results for 
different cases are presented. The procedure developed has been applied to several 
realistic problems. Results of selected studies are presented in the appendices of this 
report. 
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1. INTRODUCTION 


In the past two decades, a tremendous progress has been made In the 

field of radiative energy transfer in nonhomogeneous nongray gaseous systems. 
As a result, several useful books [1-18] and review articles [19-26] have 
become available for engineering, meteorological, and astrophysical 

applications. In the sixties and early seventies, radiative transfer analyses 
were limited to one-dimensional cases. Multidimensional analyses and 
sophisticated numerical procedures emerged in the mid-to-late seventies. 
Today, the field of radiative energy transfer In gaseous systems is getting an 
ever increasing attention because of Its applications In the areas of the 

earth's radiation budget studies and climate modeling, fire and combustion 

research, entry and reentry phenomena, hypersonic propulsion and defense- 
oriented research. 

In most studies involving combined mass, momentum, and energy transfer, 
the radiative transfer formulation has been coupled only with the steady 
processes. The goal of this research is to include the nongray radiative 
formulation in the general unsteady governing equations and provide the step- 
by-step analysis and solution procedure for several realistic problems. The 
specific objective of the present study Is to Investigate the one -dimensional 
transient radiative transfer In a nongray gaseous system. In the future work, 
the present analysis will be extended (in a systematic manner) to the problems 
of combined transfer processes in chemically reacting flows. 

For the present study, the Information on band absorption and correlation 
is summarized in section 2 and fundamental radiative flux equations are 
presented in section 3. The basic formulation for the transient radiation Is 
given in section 4, and this is applied to a special case in section 5. The 
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solution procedures are described In section 6, and plans for obtaining 
specific results are presented In section 7. 

2. BAND ABSORPTION AND CORRELATIONS 

The study of radiative transmission In nonhomgeneous gaseous systems 

requires a detailed knowledge of the absorption, emission, and scattering 

characteristics of the specific species under Investigation. In absorbing and 

emitting mediums, an accurate model for the spectral absorption coefficient Is 

of vital Importance In the correct formulation of the radiative flux 

equations. A systematic representation of the absorption by a gas, In the 

Infrared, requires the Indentlflcatlon of the major Infrared bands and 

evaluation of the line parameters (line intensity, line half-width, and 

• • 

spacing between the lines) of these bands. The line parameters depend upon 
the temperature, pressure and concentration of the absorbing molecules and, in 
general, these quantities vary continuously along a nonhomogeneous path In the 
medium. In recent years, considerable efforts have been expended in obtaining 
the line parameters and absorption coefflcents of Important atomic and 
molecular species [27-30], 

For an accurate evaluation of the transmittance (or absorptance) of a 
molecular band, a convenient line model Is used to represent the variation of 
the spectral absorption coefficient. The line models usually employed are 
Lorentz, Doppler, and Voigt line profiles. A complete formulation (and 
comparison) of the transmittance and absorptance by these line profiles Is 
given In [22-26], In a particular band consisting of many lines, the 
absorption coefficient varies very rapidly with the frequency. Thus, It 
becomes very difficult and time-consuming task to evaluate the total band 
absorptance over the actual band contour by employing an appropriate line 
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profile model. Consequently, several approximate band models (narrow as well 

as wide) have been proposed which represent absorption from an actual band 
with reasonable accuracy [22-26, 31-40]. Several continuous correlations for 
the total band absorption are available In literature [22-26, 36-40]. These 
have been employed In many nongray radiative transfer analyses with varying 
degree of success [22-26, 41]. A brief discussion Is presented here on the 
total band absorption, band models, and band absorptance correlations. 

The absorption within a narrow spectral Interval of a vibration rotation 
band can quite accurately be represented by the so-called "narrow band 
models." For a homogeneous path, the total absorptance of a narrow band is 

given by 

A n = / [l-exp(k w X)] dw (2.1) 

Au> f\ 

where k^ Is the volumetric absorption coefficient, w Is the wave number, 
and X = py Is the pressure path length. The limits of Integration In Eq. 
(2.1) are over the narrow band pass considered. The total band absorptance of 
the so-called "wide band models" Is given by 

ao 

A - / [l-expf-kJO] d(u.-u) 0 ) (2.2) 

—CD 

where the limits of Integration are over the entire band pass and u Q Is the 
wave number at the center of the wide band. In actual radiative transfer 
analyses, the quantity of frequent Interest Is the derivative of Eqs. (2.1) 
and (2.2). 

Four commonly used narrow band models are Elsasser, Statistical, Random 
Elsasser, and Quasi-Random. The application of a model to a particular case 
depends upon the nature of the absorbing emitting molecule. Complete 
discussion on narrow bands models, and expressions for transmittance and 
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integrated absorptance are available in the literature [22-26, 31-33]. 

Detailed discussions on the wide band models are given in [22-26, 34-40]. The 
relations for total band absorptance of a wide band are obtained from the 
absorptance formulations of narrow band models by employing the relations for 
the variation of line intensity as [22-26, 37-40] 

Sj/d * (S/A o )exp{[-b 0 |w-(u o | ]/A q ) (2.3) 

where Sj is the intensity of the jth spectral line, d Is the line spacing, S 
Is the integrated intensity of a wide band, A Q is the band width parameter, 
and b 0 * 2 for a symmetrical band and b 0 * 1 for bands with upper and lower 
wave number heads at u Q . The total absorptance of an exponential wide band. 
In turn, may be expressed by 

*(u.p) = A(u.3)/A 0 - f w1de [A M (u.3)]d(«-w o ) (2.4) 

0 band 

where u * SX/A 0 is the nondimensional path length, p = 2ny^/d is the line 
structure parameter, y^ is the Lorentz line half-width, and A^fu.P) repre- 
sents the mean absorptance of a narrow band. 

By employing the Elsasser narrow band absorptance relation and Eq. (2.3) 
the expression for the exponential wide band absorptance is obtained as 
[25,16] 

71 

A(u,p) - Y ♦ (1/n) / [In 4 + Ej (4>)] dz (2.5) 

0 

where 4> * u sinh p/(cosh p- cos z), y * 0.5772156 is the Euler's constant, 
and Ej(4>) is the exponential integral of the first order. Analytic solution 
of Eq. (2.5) can be obtained in a series form as [25, 26] 

co 

A(u,p) - E {-(A) n [SUM(mn) ]/[n(B+l) n n! (n-1) ! ]) 
n*l 


( 2 . 6 ) 
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where «, 

SUM (mn) = E [(n+m-1)! (2m-l) !C m ]/(2 m (m! T3 
m=o 

A = - u tanhp, B ■ l/coshp, 

C = 2/ ( 1+coshp) = 28/ (B+l ) . 

The series In Eq. (2.6) converges rapidly. When the weak line approximation 
for the El sasser model Is valid (l.e. p Is large), then Eq. (2.5) reduces to 

X(u) = y + 1 n ( u ) + Ej(u). (2.7) 

In the linear limit, Eqs. (2.5) and (2.6) reduce to A = u, and In the 
logarithmic limit they reduce to X 3 Y + ln(u). It can be shown that Eq. 
(2.5) reduces to the correct limiting form In the square-root limit. Results 
of Eqs. (2.5) and (2.6) are found to be Identical for all pressures and 
pa thlengths. For p > 1 atm, results of Eqs. (2.5)-(2.7) are In good agreement 
for all path lengths. 

By employing the uniform statistical, general statistical, and random 
Elsasser narrow band models absorptance relations and Eq. (2.3), three 
additional expesslons for the exponential wide band absorptance were obtained 
In [25, 26]. The absorptance results of the four wide band models are 

discussed In detail In [26]. The expression obtained by employing the uniform 
statistical model also reduces to the relation (2.7) for large p. 

Several continuous correlations for the total absorptance of a wide band, 
which are valid over different values of path length and line structure 
parameter, are available In the literature. These are discussed, In detail. 
In [22-26, 37-40] and are presented here In the sequence that they became 
available In the literature. Most of these correlations are developed to 
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satisfy at least some of the limiting conditions (nonoverlapping line, llenar, 
weak line, and strong line approximations, and square-root, large pressure, 
and large path length limits) for the total band absorptance [23-26]. Some of 
the correlations even have experimental justifications [22,35]. 

The first correlation for the exponential wide band absorptance (a three 
piece correlation) was proposed by Edwards et al. [34, 35]. The first 
continuous correlation was proposed by Tien and Lowder [22], and this is of 
the form 


S(u,p) * ln(uf ( t) (u+2)/[u+2f ( t) ]}+l) (2.8) 

where 

f(t) => 2.94[l-exp(-2.60t) ] , t * P/2. 

This correlation does not reduce to the correct limiting form In the square- 
root limit [23,26], and Its use should be made for p > 0.1. Another 
continuous correlation was proposed by Goody and Belton [39], and in terms of 
the present nomenclature, this Is given by 

A(u,P) = 2 ln(l+u/[4+(*u/4t)]^}, P = 2t. (2.9) 

Use of this corelation Is restricted to relatively small p values [23-26]. 
Tien and Ling [40] have proposed a simple two parameter correlation for 
A(u,p) as 

X(u) * slnh - * (u) (2.10) 

which Is valid only for the limit of large p. A relatively simple 

continuous correlation was introduced by Cess and Tlwarl [23], and this is of 


the form 



7 


A(u,p) = 2 ln(l+u/(2+[u(l+l/p)] 1/2 }) (2.11) 

where p 3 4t/n * 2p/n. By slightly modifying Eq. (2.11), another form of 
the wide band absorptance Is obtained as [25, 26] 

A(u,P) = 2 ln(l+u/(2+[u(c+n/2p)] 1/2 }] (2.12) 


where 


1 0.1, p < 1 and all u values 

0.1, P > 1 and u < 1 

0.25, p > 1 and u > 1. 


Equations (2.11) and (2.12) reduce to all the limiting forms [23]. Based on 
the formulations of slab band absorptance, Edwards and Balakrlshnan [37] have 
proposed the correlation 

A(u) = ln(u) + Ej(u) + y + j - Ej(u) (2.13) 

which Is valid for large p. For present application, this correlation 
should be modified by using the technique discussed In [25, 26]. Based upon 
the formulation of the total band absorptance from the general statistical 
model, Felske and Tien [38] have proposed a continuous correlation for 
A(u,p) as 

A(u,p) = ^E^tpJ + Ej(p u / 2) - E 1 [(p u /2)(l+2t)] 

+ ln[tp u ) 2 /(l+2t)] + 2y (2.14) 

where 

p u = Ut/u)[l + (t/u)]}" 1/2 
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The absorptance relation given by Eq. (2.7) is another simple correlation 
which is valid for all path lengths and for t * (P/2) > 1. The relation of 
Eq. (2.6) can be treated as another correlation applicable to gases whose 
spectral behavior can be described by the El sasser model. In [26] Tiwari has 
shown that the Elsasser as well as random band model formulations for the 

total band absorptance reduce to Eq. (2.7) for t > 1. 

Band absorptance results of various correlations are compared and 
discussed in some detail in [25, 26, 41]. It was found that results of these 

correlations could be in error by as much as 401 when compared with the exact 

solutions based on different band models. Felske and Tien’s correlation was 
found to give the least error when compared with the exact solution based on 
the general statistical model while Tien and Lowder's correlation gave the 

least error when compared with the exact solution based on the Elsasser model. 
The results of Cess and Tiwari' s correlations followed the trend of general 
statistical model. Tiwari and Batki's correlation [Eq. 2.6 or 2.7] was found 
to provide a uniformly better approximation for the total band absorptance at 
relatively high pressures. The sole motivation in presenting the various 

correlations here is to see if their use in actual radiative processes made 
any significant difference in the final results. 

In reference 41, use of several continuous correlations for total band 
absorptance was made to two problems to investigate their influence on the 
final results of actual radiative processes. For the case of radiative 
transfer in a gas with internal heat source, it was found that actual center- 
line temperature results obtained by using the different correlations follow 
the same general trend as the results of total band absorptance by these 
correlations. From these results, it may be concluded that use of the Tien 
and Lowder's correlation should be avoided at lower pressures, but its use is 
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justified (at moderate and high pressures) to gases whose spectral behavior 
can be described by the regular Elasasser band model. For all pressures and 
path length conditions, use of the Cess and Tlwarl’s correlations could be 
made to gases with bands of highly overlapping lines. In a more realistic 
problem Involving flow of an absorbing emitting gas, results of different 
correlations (except the Tien and Lowder's correlation) differ from each other 
by less than 6? for all pressures and path lengths. Use of Tien and Lowder's 
correlations Is justified for gases like CO at moderate and high pressures. 
For gases like COj, use of any other correlation is recommended. While Felske 
and Tien's correlation Is useful for all pressures and path lengths to gases 
having random band structure. Tlwarl and Batkl's simple correlation could be 
employed to gases with regular or random band structure but for P > 1.0 atm. 

3. RADIATIVE FLUX EQUATIONS 

For many engineering and astrophyslcal applications, the radiative 
transfer equations are formulated for one-dimensional planar systems (Fig. 
3.1). For diffuse boundaries and In the abscence of scattering, expressions 
for the radiative flux and Its derivative are given as [8] 

< W T \ ) = 2 B 1X E 3 {t X ) " 2 B 2X E 3 (t oX " T X ) 


+ 2 C/ Tx e bx (t) E 2 (t x - t)dt - / T ° X e bx (t)E 2 (t - x x )dt] 
0 \ 


and 


(3.1) 


dq 


R\ 


2 B 1X E 2 (x X ) + 2 B 2X E 2 (t oX " V 


+ 2 C/ T ° X e bX (t)E 1 (|x x - t|)dt 
0 


(3.2) 
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where 


y l 

T x a / k x d * • x 0 x = / k x d y 

0 0 

E n (t) =* J u n 2 e ^ d|i 
o 


(3.3a) 

(3.3b) 


In the preceding equations, E n ( t ) are the exponential Integral functions, and 
and t ^ represent the optical coordinate and optical path, respectively. 
The quantitltes B^ and 8^ represent the spectral surface radiosltles and 
for nonreflecting surfaces, Bj^ = e^ = e bi\ etc. Thus, for non- 

reflecting boundaries, Eqs. (3.1) and (3.2) are expressed In terms of the wave 
number as (see Appendix A) 

q Rto {x J “ e lu " e 2w 

+ 2 [/*" F . ( t) E 2 (t w -t)dt - /^ 0W F 2 (t) E, (t-xjdt] (3.4) 

O U) c 

and 

dq R«u 

- = - 2 [F. (t ) + F, (t )] 

dT^ la) 0J 2(A) 0) 

♦ 2 £/'" F la, (t) E 1 ( V t)dt + ^T° W F 2a, (t) E l (t ~\> )dt J (3 ' 5) 

o io 

where 

F l» (t) * e J t] • e lo, ; V * e u (t) ' e 2» 


Equations (3.4) and (3.5) are the general equations for one-dimensional 
absorbing-emitting medium with diffuse non-reflecting boundaries. For nongray 
analyses. It Is often convenient to replace the exponential Integrals by 



11 


appropriate exponential functions [6, 8]. Upon employing the exponential 
kernal approximation [8] 

E 2 (t) = | exp (- t); Ej(.t) * | exp (- | t) 

Eqs. (3.4) and (3.5) are expressed In physical coordinate tes as 
^Ro) (y) * e lo, - e 2w 

+ 7 % F l(jj (z) k a, ex P t-7 k w (y " 2)] dZ 
• 1 f f 2u (z) k w exp c_ ? k u, (z -y )] dz {3 * 6) 

" ‘ 2 CF lw (y) + F 2a) (y);1 

+ T /Y F Ioj {Z) k « ex P<>l k u (y_z)]dz 
0 

+ I J F 2u> ( z) k w exp[ " 7 k a) (z " y)]dz (3 - 7) 

where 2 Is a dummy variable for y. However, by differentiating Eq. (3.6) 
directly, there Is obtained 

' ‘ar^ * ' I k to CF lu> (y) + F 2w (y)] 

+ T /V F 1 W (z) k ! expl> 7 k J y ' z)]dz 
0 

+ */ y F 2u> k ! expC- | k w ( 2 -y)]dz 


(3.8) 
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The slight difference In Eqs. (3.7) and (3.8) should be noted. This Is a 
consequence of using the exponential kernal approximation. If one has to make 
a decision as which equation to use. It Is recommended to use Eq. (3.8). 

The total band absorptance, as given by Eq. (2.2), can be expressed In a 
slightly different form as 

<o 

A(y) * / [1 - exp(- k u y)] du» ~ cm" 1 (3.9a) 

o 

where both k^ and w have units of cm"*. Differentiation of Eq. (3.9a) gives 

00 -? 

A' (y) = / k w exp(- k^ y) doj ~ cm (3.9b) 

o 

and 

oo 

A" (y ) * / - k* exp(- k^ y) dw ~ cm” 3 (3.9c) 

o 

Equations (3.9) are employed to express Eqs. (3.6) and (3.8) In terms of the 
band absorptance. 

The total radiative flux Is given by 



(3.11) 
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Upon substituting Eq. (3.6) Into Eq. (3.10) and Eq. (3.8) into Eq. (3.11) 

there is obtained for a multi band gaseous system 
Q„(y) - e, - e 2 


* 7 L , {/y F l»l <z) eXp[ " 7 k u i (y - z)]dz 

1 = 1 Awl 0 
L 1 

-I F 2ui (z) k w 1 expC ' 7 k u>1 (z " y):Idz}d<i, 1 (3.12) 


dq R (y) 

3y 


It should be 
obtaining Eqs. 


7 I , k o,|C F l^ ( 7> + l W>'» d ” t 

1=1 Awl 

I , r i„) (z > k »l expC- 4 k „ ( «y-z>3<Jz 

1 = 1 Atol 0 

+ / F 2a)i (z) k f>i e *pc-l k «i( 2 -y) 3dz >d- f (3.13) 

y 

pointed out that the following relations have been used in 
(3.12) and (3.13) 

09 00 

/ e lu d» * e t ; / e 2u du * e 2 
0 0 

/ {/ y F iJ z) k u expC- | kjy-z)]dz}du> 

0 0 

* / . </ y F l„1 (z) k „l «p [ - 1 k u l (y - z) dz,d “t 

1 = 1 Awl 0 


where n represents the number of bands in a multiband system. 
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By util zing the definitions of the band absorptance and Its derivatives 
as given by Eqs. (3.9) and evaluating the value of the Planck function at the 
center of each band, Eqs. (3.12) and (3.13) are expressed as 

q R (y) ■ 6] - e 2 


*4 " (/ F,„ (2) a; ci (y-z)Jcfz 

* 1=1 o 1 ol 1 * 


- / F- (z) a; d (z-y ) ]dz) 
y %1 1 1 


(3.14) 


k “< du '’ 


+ !■ " U ^ (2) A*; t| (y-z)]dz 


1=1 0 A ol 


+/ F- (z) A” d (z-y)]dz> 

y 1 * 


(3.15) 


where u> 0 ^ represents the center of the 1th band. 

Equations (3.14) and (3.15) are In proper form for obtaining the nongray 
solutions of molecular species. However, In order to be able to use the band 
model correlations, these equations must be transformed In terms of the 
correlation quantities defined In Eq. (2.4). The following quantities, 
therefore, are needed for the transformation 

u = (S/A q ) py; u q » (S/A q ) PL; PS * J k^ do> (3.16) 

Au> 

Now, by using the definition X = A/A q , Eq. (3.9b) Is written as 
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A'(y) -^ [*<»>] • TT ^ ■ A'(y)/A 0 - cm' 1 

Thus, 

A*(y) - A o ^M= A o C^iHl^] = P S(T) A' (u) (3.17a) 


Similarly 

A* * (y) = [P S(T)] 2 (l/A o ) A”(u) (3.17b) 

The dimensions of both sides In Eqs. (3.17a) and (3.17b) agree with the 
dimensions given In Eqs. (3.9b) and (3.9c). By employing the definitions of 
Eqs. (3.16) and (3.17), Eqs. (3.14) and (3.15) are expressed as 

q R (u) - ej - 


+ 


3 

7 


n 

E 

1=1 


dq R (u) 

iilu 



n 

E 

1-1 


n 

l 

1-1 


u, 

*0) </ F 


- 1 


oi 


A . [F. 
oi lw. 


A „. </ 


* / 


'ot 


i», (u i> *i [ 1 ,u i ' u i ,] d “i 

(u i > l \ c l ,u i - V 1 du i> 

(u) * F 2 c 1 (u)] 

F l„ | ' u 'l A", [j (u, - uj)] duj 

F 2u. 1 lu ' ) l '\ [ 7 ,u i - “i ),du i ) 


(3.18) 


(3.19) 


where u‘ Is the dummy variable for u and A*(u) = dA/du. It should be noted 
and F 2w In Eqs. (3.18) and (3.19) represent the values of 
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Fj u and at the center of the 1th band, and dq R /dy = (dq R /du) (du/dy ) * 
[P S(T)/A q ] (dq R /du) . 

By defining the new Independent variables as 

l - u/u Q » y/L; V = u'/u o = z/L (3.20) 

Eqs. (3.18) and (3.19) can be expressed as 
q„ (V - - e 2 



(3.21) 


(3.22) 


where again X'(u) dentoes the derivative of A(u) with respect to u, and 
dq R /du * (dq R /dO (d£/du) = (l/u Q ) (dq R /d£). 

Equations (3.18) through (3.22) allow us to make use of the band model 
correlations for the wlde-band absorptance because these correlations are 
expressed In terms of u and p. However, It Is often desirable and 
convenient to express the relations for q R and dlv qp which only Involve 
X(u) and X'(u) but not X''(u). This Is accomplished by Integrating the 
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Integrals In the expressions for qq and dlv qq by parts. This results In 
simpler Integrals. Upon performing the Integration by parts on the Integrals 
In Eqs. (3.21) and (3.22), the equations can be expressed In alternate forms 
as (see Appendix B) 

q„U> • «! - 


+ 


dq R (0 
— 3T" 


ss 


n 

S 

1=1 


3 

7 


5 . •, 

A q i (/ Cde^ (5’)/dC'] [|u oi U-£') d5‘ 

0 1 

1 . . 

+ / [de w (5')/d5‘] A^ [| u q1 U'-OMe'} 

( ", »ot %, J E«V 5 ' )/«'J *i [7 

1 = 1 0 1 

1 _ , 

- / [de^U*)/d 5 '] a; [7 


(3.23) 


(3.24) 


It should be noted that Eq. (3.24) can be obtained directly by differentiating 
Eq. (3.23) with respect to t, using the Leibnitz formula. This Is shown In 
Appendix B. 

Equat1ons(3.21) , (3.23), and (3.24) are the most convenient equations to 
use when employing the band-model correlations In radiative transfer analyses. 
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4. BASIC FORMULATION FOR TRANSIENT PROCESSES 

The Interaction of radiation In transient transfer processes has received 
very little attention In the literature. Yet, the transient approach appears 
to be the logical way of formulating a problem In general sense for elegant 
numerical and computational solutions. The steady-state solutions can be 
obtained as limiting solutions for large times. 

A few studies available on radiative Ineractlons reveal that the 
transient behavior of a physical system can be influenced significantly In the 
presence of radiation [42-45], Lick Investigated the transient energy 
transfer by radiation and conduction through a semi-finite medium [42]. A 
kernal substitution technique was used to obtain analytic solutions and 
display the main features and parameters of the poblem. Doornlnk and Herlng 
studied the transient radiative transfer In a stationary plane layer of a 
nonconducting medium bounded by black walls [43]. A rectangular Milne- 
Eddlngton type relation was used to describe the frequency dependence of the 
absorption coefficient. It was found that the cooling of the layer Initially 
at a uniform temperature is strongly dependent on the absorption coefficient 
model employed. Larson and Vlskanta Investigated the problem of transient 
combined laminar free convection and radiation In a rectangular enclosure 
[44]. It was demonstrated that the radiation dominates the heat transfer In 
the enclosure and alters the convective flow patterns significantly. The 
transient heat exchange between a radlatldlT plate and a high- temperature gas 
flow was Investigated by Melnikov and Sukhovlch [45]. Only the radiative 
Interaction from the plate was considered; the gas was treated as a non- 
participating medium. It was proved that the surface temperature Is a 
function of time and of longitudinal coordinate. 
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The objective of this study Is to investigate the Interaction of nongray 
radiation In transient transfer processes In a general sense. Attention, 
however, will be directed first to a simple problem of the transient radiative 
exchange between two parallel plates. In subsequent studies, the present 
analysis and numerical techniques will be extended to Include the flow of 
homogeneous, nonhomogeneous, and chemically reacting species In one- and 
mul tl -dimensional systems. 

The physical model considered for the present study is the transient 
energy transfer by radiation In absorbing-emitting gases bounded by two 
parallel gray plates (Fig. 4.1). In general, Tj and T 2 can be a function of 
time and position and there may exist an initial temperature distribution in 
the gas. It Is assumed that the radiative energy transfer In the axial 
direction Is negligible In comparison to that In the normal direction. 

For radiation participating medium, the equations expressing conservation 
of mass and momentum remain unaltered, while the conservation of energy, In 
general, is expressed as [8] 

p c p JR « div (k grad T) + pT + p 4 - dlv q R (4.1) 

where p Is the coefficient of thermal expansion of the fluid and <t> Is the 
Rayleigh dissipation function. For a semi-infinite medium capable of 
transferring energy only by radiation and conduction, Eq. (4.1) reduces to 

■>%«•-$ «- 2 > 

where q Is the sum of the conductive heat flux q c * - k OT/3y) and the 
radiative flux q^. For the physical model where radiation Is the only mode of 
energy transfer, the energy equation can be written as 
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P c 


QT = dq R 
p cST Sy 


(4.3) 


Use of this simplified equation Is made to Investigate the transient behavior 
of a radiation participating medium. 

As pointed out In the previous section (Sec. 3), Eqs. (3.21) and (3.23) 
are convenient equations for the radiative flux. Equations (3.22) and (3.24) 
are two expressslons for the dlv q R (y), but Eq. (3.24) Is preferred because It 
only Involves the first derivative of A and avoids singularities In the 
large path length limit. 

Upon defining nondlmenslonal radiative heat flux by 


Q U,t) - q R (C, t)/[e x ( t) - e 2 (t)] 


(4.4) 


Eq. (3.21) can be written as 

. n 


Q U,t) =1+| £ A Qi u q1 {/ (5\t) Aj [2 U Qi U-5')]d5' 

1=1 o 


1 


/ c 2i (s\t) a; [| u Qf (r-^]d^'} 


(4.5) 


where 




Nat<. : 

^ A ol C- -V n 0n <},~rr*;^ ^ p ^ ' 

- I 


w<tt 


C 1 (5 * t) * F w i U.tJ/Cejtt) - e 2 (t)] 


A** ^ c 

^ ^ ^ vJq t-t * ~ 

Equation (4.5) provides the general expression for the radiative flux in the 
nondlmenslonal form. A similar nondlmenslonal form can be obtained also from 
Eq. (3.23). 


By defining 4 (5,t) = T(5,t)/T Q with T Q representing some constant 
reference temperature, Eqs. (4.3) and (3.24) can be combined to yield the 
energy equation In nondlmenslonal form as 


- <?*- 





- a* U,t)/at » \ " {/^ <1^ U‘.t) a; [* u Ql U-S')d5' 
1=1 o 


1 - i 

- / * u1 u\t) a; [^u ot (^ , -5)]d5*} (4.6) 

where 

•^(S.t) = <PS i <T)Ca e ui (e,t)/dg]/(p C p ^/t^)} 3= 

The time t In Eq. (4.6) is defined as t* * t/t m with t,,, representing some 

characteristic time scale of the physical problem; however, for the sake of 

convenience, the asterisk Is left out here as well as In further developments. 

From the definitions of 4»U,t) and (I^U.t). It should be noted that Eq. 

(4.6) is a nonlinear equation in T(5,t). Equation (4.6), therefore, 

i.ransi*t 

represents a general case of the transient energy by radiation between two 
semi-infinite parallel plates. A similar expression can be obtained also by 
combining Eqs. (4.3) and (3.22). 


5. A SPECIAL CASE OF TRANSIENT INTERACTION 


As a special case. It is assumed that the entire system Initially Is at 
the fixed (reference) temperature T Q . For all time, the temperature of the 
upper plate Is maintained at the constant temperature equal to the reference 
temperature, l.e., T 2 = T Q . The temperature of the lower plate Is suddenly 
decreased to a lower but constant temperture, l.e., Tj < T 2 . The problem, 
therefore, Is to Investigate the transient cooling rate of the gas for a step 
change in temperature of the lower plate. 


Since small temperature differences have been assumed and the absorption 

coefficient has been taken as Independent of temperature one may ^mploy ’ 

✓ 


additionally the linearization. 


( 
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* u1 < T > - W * ,d e »i /dTl r (T - r »' ,5 -‘i 

w 

where again the subscript 1 refers to the 1th band such that <*>j is the wave 
number location of the band and T w represents the temperature of the reference 
wall which could be either Tj or T 2 . For the special case considered, since 
we are Interested In investigating the transient behavior of the gas because 
of a step change In temperature of the lower plate, T w Is taken to be equal to 
Tj. Thus, 

• (d e uf /dT) Ti (T-Tj) (5.2a) 

e w1 (l.t) - e ui (0,t) = (d e u1 /dT) T (T^Tj) (5.2b) 

e w1 (5.t) - e ui (l,t) * (d e w( /dT) T (T-T 2 ) (5.2c) 

Mote that Eq. (5.2c) Is obtained by subtracting Eq. (5.2b) from Eq. (5.2a). 
Also, for linearized radiation, 

T 4 * 4 Tj T - 3 Tj (5.3) 

Thus, ej * a T 4 , e^ * o T 4 ■ a (4 T^ T 2 - 3 T 4 ) such that 

e i ~ e 2 “ 4 a Tj (T j- T 2 ). 

It should be pointed out that for a single-band gas, the linearization Is 
not required because the temperature distribution can be obtained either by 
combining Eqs. (3.22) and (4.3) or from Eq. (4.6) and the radiative heat flux 
can be calculated from Eqs. (3.21), (3.23), or (4.5). However, for the case 
of multiband gases and for systems Involving mixtures of gases. It Is 
convenient to employ the linearization procedure In order to use the 
Information on band model correlations. The following definition are useful 
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in expressing the governing equations in linearized forms: 


9 = (T-T 1 )/(T 2 -T 1 ) 


(5.4a) 


N 11 3 {Pt m /p c p ) K li • K li = S i (T) (d e u>i /dT) T, (5 * 4b) 


■l 3 (p V p V K 1 • K 1 3 A K ii 


i=l 


(5.4c) 


M li 3 (t m /L p c p> H 11 * H li = A oi (T) (d e a>i /dT) Tl (5 * 4d) 


M 1 3 ( V L p C p> H 1 • H 1 = H li 


( 5.4e) 


M li u oi “ N li ’ u oi H li " PL K li 


(5.4f ) 


where H^, Kj, Nj and Mj represent the values of H, K, N and M evaluated at the 
temperature T^. As explained in Refs. 8 and 23, these quantities represent 
the properties of the gaseous medium. 

By employing the definitions of Eqs. (5.2) - (5.4), relations for the 
radiative flux, as given by Eqs. (3.21) and (3.23), are expressed as 

0 


Q(£,t) - 1- (3/8 a Tj) E u q1 {/' e(S'.t) Aj[| u Qf (5-5 ' ) W 


1 


( fleie 

/(vr*) 


+ / [1 - 0(£',t)] Aj [-| u q ^ (£' -£) ]d£' ) <t (5.5a) 


and 


QU,t) = 1 - (1/4 a Tj) E H 21 {/ C 59( |^f^ A^l u q1 [l-V ) ]d5' <&= *W, 


t rr r r\ f ♦ /»< J 
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- J 1 u of (5'-5)]d5'> 


(5.5b) 


Thus, the expressions for the heat flux at the lower wall are given by 


1 


Q(0,t) » 1 
and 

Q(0,t) 


(3/8 o Tj) E u Qf H lf / [1-0(5* ,t)]A f (| U Qf 5*)d5’ 
1 =*1 0 


1 ♦ (1/4 o Tj ) H 1( / s )( | « ol ?•>«• 


(5.6a) 

<4= t\Joi t ; Hii /(& ) 

IS non dt m fsiZt 'on**. 4 

(5.6b) 


It should be pointed out that Eqs. (5.5a) and (5.6a) are convenient forms for 
the optically thin and general solutions while Eqs. (5.5b) and (5.6b) are 
useful for solutions In the large path length limit. Once the solutions for 
9(5»t) are known from the energy equation, the appropriate relations for the 
heat flux can be obtained from Eqs. (5.5) and (5.6). 


By employing the definitions of Eqs. (5.2) - (5.4), a combination of 
Eqs. (3.22) and (4.3) provides one form of the energy equation and Eq. (4.6) 
Is transformed to obtain another form; these are expressed as 


S9 ^ ,t) + 3 Nj 0(5, t) - | Nj = 


A Jote : 

yl i ‘ is 


T 1 M 11 u of Aj'[| u o1 (5-5’)]d5' 

1 * 1 0 


+ J [0(5'. t) - 1] A'j[| u oi (5'-5)]d5‘> aW > mms i «•* (5.7a) 


and 


M(t.t) . 3 " „ 89(5'. t) ;, r 3 .. 

—ST" 2 "ll (/ if Vf U ot <5 5 ,]de 
1-1 0 


fJoif 


^ /V // ft 
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-/ -- — - i - A;[| u ,(5‘-5)]d5‘} (5.7b) 

o£ 1 t ol 

The initial and boundary conditions for Eq. (5.7) are specified as 

9(5,0) = 1 ; 9(0, t) = 0 ; 0(1, t) * 1 (5.8) 

The parameters in Eq. (5.7) are Nj and u Q . For a given gas, the parameters 
are the gas pressure and the temperature of the lower wall. Equation (5.7b) 
is the convenient form for solutions in the large path length limit. 

6. METHOD OF SOLUTIONS 

For the general case, the temperature distribution Is obtained from the 
solution of the energy equation, Eqs. (5.7). Once 9(5, t) Is known, the 
radiative heat flux Is calculated by using the appropriate form of Eq. (5.6). 
Before discussing the solution procedure for the general case, however. It is 
desirable to obtain the limiting forms of Eqs. (5.5) and (5.7) in the 
optically thin and large path length limits and investigate the solutions of 
resulting equations. 


'T 


6.1 Optically Thin Limit 

In the optically thin limit [8, 23], A(u )J = u, A‘(u) » 1, 

A*'(u) * 0. In this limit, therefore, Eq. (5.7a) reduces to 


and 


<$: — 


+ 3 Nj 9(5, t) - | Nj * 0 


(6.1a) 


From an examination of Eq. (6.1a) along with the definitions given In Eq. 
(5.4), it is evidlent that in the optically thin limit the temperature distri- 
bution In the medium is independent of the 5~ coordinate. This Is a 
charateristic of the optically thin radiation [8, 23]. Thus, Eq. (6.1a) can 
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be written as 


d0(t) 

dt 


+ 3 N (t) 0(t) - | M 1 (t) 


0 ; 0 ( 5 , 0 ) = 1 


Since gas properties are evaluated at known reference conditions, 
essentially constant, and solution of Eq. (6.1b) is found to be 


(6.1b) 
N x is 


0(t) = j[l + exp(- 3 Njt)] 


( 6 . 2 ) 


In the optically thin limit, Eq. (5.7b) reduces to 
90(5. t) _ 3 # ; u wf 5 90(5* »t) / 90(5\t) Jp( , 

- -at— 2 "k 1 {/ — 9T“ \ —W~ 61 } (6 ' 3a) 

i=l o 5 

A differentiation of Eq. (6.3a) with respect to 5 (by using the Leibnitz's 
rule) results in 

jg [ d9 ^ ,t) ♦ 3 N x 0(5. t)] = 0 (6.3b) 


or 

+ 3 Nj 0(5, t) * C(t) (6.3c) 

The constant of integration C(t) is evaluated through the combination of Eqs. 
(5.8) and (6.3a) and is found to be C(t) * Nj. A substitution of this in 
Eq. (6.3c) gives Eq. (6.1a) for which the solution is given by Eq. (6.2). 
Thus, as would be expected, both general forms of the energy equation reduce 
to the same equation in the optically thin limit. 

In the optically thin limit Eqs. (5.5a) and (5.5b) respectively reduce to 
i 5 1 

Q(5,t) * 1 - [3/{8aT, )] (PLK.) {/ 0(5',t)d5' + / [1-0(5' ,t) ]d5' ) (6.4a) 

1 1 o 5 


and 
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5 1 

Q(5.t) = 1 - [3/8 otJ)] (PLK.) {/ (5-5’)d5' + / Ur (5'-Od5'> (6.4b) 

1 1 o5 _ ot, 

o 5 

Through Integration by parts. It can be shown that Eq. (6.4b) reproduces Eq. 

(6.4a). 8y noting that. In the optically thin limit, 9(5, t) = 9( t) , Eqs. 

(6.4) can be expressed as tibi' : (PLK X /(TT* ^ ' 


nondt ^7 rslS' * 

Q(5,t) = 1 - [3/ (8 cjTj ) ] (PLKj) [(1-5) + (25-1) 9(t)] (6:5) 

It should be pointed out that Eq. (6.5) can be obtained directly from Eq. 
(6.4b) without performing the Inegratlon by parts. The heat transfer from the 
lower surface In the optically thin limit, therefore. Is given by 




Q(0,t) - 1 - [3/(8orp] (PLKjJCl - e( t) ] •' 


A 


( 6 . 6 ) 


The result of Eq. (6.6) can be obtained directly by letting 5=0 In either 
of Eqs. (6.4). The relation for 0(t) In Eq. (6.6) Is obtained from Eq. 
(6.2). Thus, evaluation of the temperature distribution and radiative heat 
flux In the optically thin limit does not require numerical solutions. 


6.2 Large Path Length Limit 

In the large path length limit (l.e., for u Q ^ >> 1 for each band), one 
has A(u) = Jln(u), fl'(u) ■ 1/u, and A* ’ (u) * - 1/u^ [8, 23], Thus, In the 


large path length limit, Eq. (5.7a) reduces to 
a9 ^^ t) + 3 9(5, t) - | Nj 


tfokt \ Af; S***/ N I 4ft 

-1 


{/^ 0 ( 5 * , t) + / [9(5*. t) - 1] 7 ) 


(6.7a) 


(5-5' r 5 


(5' -5) 


It should be noted that for any fixed value Tj and a given ga$, Nj and Mj are 


/ 


/V, = [PL 2 CS^/Acy) ] M) = [ 2 Ai, 
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constants; but, 0(5.t) does depend on 5. For a given gas and with known 
values of T 2 and t^,, the solution of Eq. (6.7a) can be obtained by specifying 
Tj. Equation (6.7a) Involves singular Integrals with Cauchy type kernals and, 
therefore, a closed form solution does not appear to be possible; numerical 
solutions, however, can be obtained by the variation of parameter technique. 
Because of the singular nature of Integrals, Eq. (6.7a) Is not a convenient 
equation for the large path length limit solutions. 


In the large path length limit, Eq. (5.7b) reduces to 

30(5. t) _ _ „ r 1 99(5'. t) d5' 

— 5t ' M l L — 5T TT^FT 


(6.7b) 


Equation (6.7b) Is a convenient form for solution In the large path length 
limit. An analytical solution of Eq. (6.7b) may be possible, but numerical 
solution can be obtained quite easily. 


In the large path length limit, Eqs. (5.5a) and (5.5b) reduce 


respectively to 


Q(5.t) ■ 1 - U/«rJ> H U W.ti - I jgrjl 

1 0 5 


tide-. ri,/(<r-r; 3 ) 

i s fiend 

(6.8a) 


cT 


and 


Q(5,t) * 1 - ( 1/ 4oTj ) E H (/ In [| U q1 (5-5‘ ) ]d5' 

1*1 0 ° 


A /<"> 7 c. • d 1 Ctnn/ 77 ! <ZfC 
by Sa 


+ / An[| u o! (5'-5)]d5'> (6.8b) 


The expressions for dimensionless radiative heat flux from or to the wall 
are obtained by setting 5*0 In Eqs. (6.8) as 
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Q(0,t) » 1 - (Hj/40^) 


1 Jr i 

/ [1 - 0(5\t)] p 


and 


(6.9a) 


No it : 


«. l<rr, ’ 

*5 nondS/fieo s ' 



0(0, t) - 1 (1/4 oTj) ! Hjj J 59( ^: t) *n(j u oi 5')d5’ (6.9b) 

Thus, once the temperature distribution is known from solutuions of Eq. (6.7), 
the wall heat flux can be calculated by using the corresponding form of Eqs. 


(6.9). 


6.3 Numerical Solutions of Governing Equations 

General solutions of Eqs. (5.7a) and (5.7b) are obtained numerically by 
employing the method of variation of parameters. For this, a polynomial form 
for 9(5, t) Is assumed in powers of 5 with time dependent coefficients as 


0(5, t) 


n 

E 

m=o 


c <t> 5 m 
m 


( 6 . 10 ) 


By considering only the quadratic solution In 5, and satisfying the boundary 
conditions of Eq. (5.8), one finds 

9(5, t) = 5 2 + g(t) (5-5 2 ) (6.11) 


where g(t) represents the time dependent coefficient. At t = 0, a combination 
of Eqs. (5.8) and (6.11) yields the result 

g ( 0) - (l-5 2 )/(5-5 2 ) (6.12) 


Also, from Eq. (6.11) there is obtained 


90(5, t) _ , 2 dg( t) 

dt * ^ dt 


(5-5 2 ) g'(t) 


(6.13) 
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and 

d9(5,t) 

• 2 5 + g(t) (1-2^) (6.14) 

o5 

Equations (6.11) - (6.14) are employed to obtain specific solutions of Eqs. 
(5.5) and (5.7). 

By substituting Eqs. (6.11) and (6.13) In Eq. (5.7a), there Is obtained 

g’(t) + G j(5) g(t) = G 2 U) (6.15) 

where the Integral functions Gj(5) and G 2 U) are defined In Appendix C. The 
solution of Eq. (6.15) Is given by 

g ( t ) * c exp[- Gj(5)t] + G 2 (5 )/Gj( 5) (6.16a) 

Since at t s 0, g(t) * g(0), then c ■ g(0) - G 2 ( 5 )/Gj( 5 ). Thus, Eq. (6.16a) 
becomes 

g(t) = Cg( 0) - G 2 (c)/G 1 (^)3 exp[- GjUJt] + G 2 (c)/G 1 (^) (6.16b) 

where g (0) Is given by Eq. (6.12). The Integrals In functions GjU) and 
G 2 (5) can be evaluated easily by numerical means, after substituting the 
relation for A' ' (u) . 

A substitution of Eqs. (6.11), (6.13) and (6.14) Into Eq. (5.7b) results 
In 

g ' ( t) + G 3 (5) g(t) = g 4 (5) (6.17) 

where the Integral functions G 3 (5) and G^(5) are defined In Appendix C. The 
solution of Eq. (6.17) Is found to be 

g(t) * [g(0) - G 4 (0/G 3 (5)3 exp[- G 3 (5)t] + G 4 (5)/G 3 (5) 


(6.18) 
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where again g(0) Is given by Eq. (6.12). 


The solutions of Eqs. (5.7a) and (5.7b) can be expressed In a convenient 
form as 


9(5, t) 


q,(t). 


- 


(6.19) 


In Eq. (6.19), gj { t) Is given by Eq. (6.16b) anti Is used for the solution of 
Eq. (5.7a) and g£(t) Is given by Eq. (6.18) and is used In obtaining the 
solution of Eq. (5.7b). The both approach should result In the same final 
solution. 

For the steady state case, the solution again is given by Eq. (6.19), but 
functions gj(t) and g 2 ( t) are no longer a function of time and are given by 

g x - G^SJ/GjU); g 2 = G 4 (5)/G 3 (5) (6.20) 

The solutions for the steady case are available In the literature and are 
useful In comparing the results of this study In the limit of t -*• ®. 

The expressions for the nondimenslonal radiative flux are obtained from a 
combination of Eqs. (5.5), (6.11) and (6.14) such that 

QU.t) - 1 - G 5 (5) gj(t) - Gg(5) (6.21a) 

and 

Q(5,t) - 1 - G ? U) g 2 (t) - Gg(5) (6.21b) 

0 

where 65(5) through Gg(5) are defined In Appendix C, and qj(t) and,Qj?(t) 
are given respectively by Eqs. (6.16b) and 6.18). Consequently, the 
expressions for the radiative heat flux at the lower wall are obtained as 
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0(0, t) = 1 + G g qj(t) - G 10 (6.22a) 

and 

Q(0, t) - 1 + G n g 2 (t) + G 12 (6.22b) 

where Gg through Gj 2 are defined in Appendix C and are not function of t,. 
It should be noted that the solutions presented In Eqs. (6.21) and (6.22) 
require the solution of the energy equation as given by Eqs. (6.19) 

6.4 Numerical Solutions of Large Path Length Equations 

As mentioned earlier, Eqs. (6.7b), (6.8b) and (6.9b) are the most appro- 
priate equations to use in the large path length limit. However, numerical 
procedure is presented for both forms of the energy and radiative flux 
equations. Once again Eqs. (6.11) through (6.14) providt the basis for 
numerical solutions also in the large path length limit. For this limit, the 
solution given by Eq. (6.19) is expressed as 

eu,t) * r + [ g ^ (t) ] u.r) (6.23) 

where ( ~t ) is used for the solution of Eq. (6.7a) and g^t) for Eq. (6.7b). 

A substitution of Eq. (6^23) into Eq. (6.7a) results In 

g 3 (t) + G 13 m g 3 (t) = G 14 m (6.24) 

where Integral functions G 13 (£) and G 14 (5) are defined in Appendix C. The 
solution of Eq. (6.24) is found to be 


g 3 (t) - Cg(0) - G 14 (5)/G 13 (^)] exp[- G 13 U)t] + G 14 U)/G 13 (5) 


(6.25) 



where g ( 0) Is defined again by Eq. (6.12). Equation (6.23) along with Eq. 

(6.25) provides the solution of the energy equation, Eq. (6.7a). 

A combination of Eqs. (6.23) and (6.7b) results In 

g; (t) + G 15 (S) g 4 (t) = G 16 U) (6.26) 

The Integral functlns G 15 (5) and G lg (5) appearing In Eq. (6.26) are defined 
In Appendix C. These, however, can be evaluated easily with the results 

Gi 5 U) = [M 1 /(5-5 2 )](2 + (25-1) An[(5-l)/5)> (6.27a) 

G 16 (5) = 2 + 5 *n[(5-l)/5]} (6.27b) 

The solution of Eq. (6.26) Is found to be 

g 4 (t) = Cg ( 0) - G 16 (5)/G 15 (5)] exp[- G 15 (5)t] + G 16 U)/G 15 (5) (6.28) 

where again g ( 0 ) is defined by Eq. (6.12). A combination of Eqs. (6.23), 

(6.27) and (6.28) provides the solution of the energy equation, Eq. (6.7b). 
The only parameter appearing In the solution of Eq. (6.7b) Is Mj. 

The expesslons for the nondimenslonal heat flux In this case is obtained 
from a combination of Eqs. (6.8) and (6.23) as 

Q(5.t) * 1 - G 1? (5) g 3 (t) - G lg (5) (6.29a) 

and 0(5, t) = 1 - G lg (5) g 4 ( t) - G 2Q (5) (6.29b) 

where Gj 7 (5) through G 2 q(5) are defined In Appendix C and can be evaluated 
In closed forms. The corresponding expressions for the radiative heat flux at 
the lower wall are found from Eqs. (6.9) and (6.23) as 


Q( 0 , t) = 1 - G 21 g 3 ( t) - G 


(6.30a) 
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and Q(0,t) = 1 - G.^ g^(t) - G24 (6.30b) 

where again G21 through G24 are defined in Appendix C and are not functions of 

5 . 
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7. RADIATIVE INTERACTION IN LArilNAR FLOWS 
The physical system considered is the energy transfer in laminar, in- 
compressible, constant properties, fully-developed flow of absorbing-emit- 
ting gases between parallel plates (Fig. 7.1). The condition of uniform 
surface heat flux for each plate is assumed such that the temperature of the 
plates varies in the axial direction. Fully developed heat transfer is 
considered, and axial conduction and radiation is assumed to be negligible 
as compared with the normal components. Consistent with the constant prop- 
erties flow, the absorption coefficient is taken to be independent of temp- 
erature and radiation can be linearized. Extensive treatment of this prob- 
lem is available in the literature [23, 41]. The primary motivation of 
studying the problem here is to investigate the extent of radiative inter- 
action for high temperature flow conditions. 

7.1 Basic Formulation 

For the physical problem considered, the energy equation, Eq. (4.1), 
can be expressed as [8J 

peril + U II + v II) = k ifl + bTu + g(iH) 2 -di vqn (7.1) 
p 3 1 3x 3y 3y 2 dx 3y 

where u and v denote x and y components of velocity, respectively. In de- 
riving Eq. (7.1) it has been assumed that the net conduction heat transfer 
in the x direction is negligible compared with the net conduction in the y 
direction. This represents the physical condition of a large value of the 
Pec let number. By an analogous reasoning, the radiative heat transfer in 
the x direction can be neglected in comparison to that transferred in the y 
direction. If, in addition, it is assumed that the Eckert number of the 
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Figure 7.1 Physical model and coordinate systems for flow of 
radiating gases between parallel plates. 
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flow is small, then Eq. (7.1) reduces to [8] 


3t 3x 3y 3y 2 pC p 3y 


(7.2) 


where a = (k/pCp) represents the thermal diffuslvity of the fluid. 

For a steady fully-developed flow, v * 0, and u Is given by the well- 
known parabolic profile as 


u = 6 u^ (£ - K 2 ); 5 = y/L (7.3) 

where u m represents the mean fluid velocity. Also, for the flow of a per- 
fect gas with uniform wall heat flux, 3T/3x is constant and is given by 

3 T/3 x = (2aq w )/(u m L/k) (7.4) 

Now, by combining Eqs. (7.2) - (7.4), the energy equation is expressed in 
nondlmens lonal form as 

4 - 

— ('J 12(C - K 2 ) = (7.5) 

3t ' 3C 2 q w 3£ 

where 

r = at/L 2 ; 0 • (T - TJ/^L/k) 

The expression for Sq^/SS In Eq. (7.5) is obtained from either Eq. (3.22) or 
Eq. (3.24). 

By assuming that the initial temperature distribution in the gas is 



some uniform value T 0 = Tj , the initial and boundary conditions for this 
problem can be expressed as 


e(e, o) = o 

(7.6a) 

9(0, t) = 0 (1, t) = 0 

(7.6b) 

(C = 1/2) = 0 ; 0^(C ■ 0) - - 0^(5 ■ 1) 

(7.6c) 


It should be noted that all the boundary conditions given in Eqs. (7.6) are 
not independent and any two convenient conditions can be used to obtain 
solutions. Also, the initial temperature distribution can be any specified 
or calculated value of 0(c, 0) = f(S). 

For flow problems, the quantity of primary interest is the bulk temper- 
ature of the gas, which may be expressed as [41] 

9b = (T b -Ti)/(q w L/k) = 6/* e(€, r) (e - K 2 ) dC (7.7) 

The heat transfer q w is given by the expression, q w = h,. (Ti-T^), where h c 
is the convective heat transfer coefficient (W/cn^-K). In general, the heat 
transfer results are expressed in terms of the Nusselt number Nu = h c D^/k. 
Here, 0^ represents the hydraulic diameter, and for the parallel plate geom- 
etry it equals twice the plate separation, i.e., 0^ = 2L. Upon eliminating 

the convective heat transfer coefficient h from the expressions for q and 

c w 

Nu, a relation between the Nusselt number and the bulk temperature is 
obtained as 
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The heat transfer results, therefore, can be expressed either in terms of Nu 
or e b . 

7.1.1 Steady Laminar Flow 

For steady-state conditions, 90/ar = 0 and Eq. (7.5) becomes 

0" - 12 (C-5 2 ) = (l/q w ) dq R /dC (7.9) 

By integrating Eq. (7.9) once and using the conditions that at c = 1/2, 
q R (C) and (d6/d£) are equal to zero, one obtains 

0' - 2(35* - 2C 3 ) + 1 - d R (e)/q w (7.10) 

The expression for q R (c) in Eq. (7.10) is obtained from either Eq. (3.21) or 
Eq. (3.23). 

For the present physical problem, ej * e 2 and F^ = . Thus, for 

the case of linearized radiation, a combination of Eqs. (3.21), (5.2a), and 
(7.10) results In 

0 ' - 2 (3c 2 - 2c 3 ) + 1 

»|(L/k) ^I H 11 u (J1 {/* 0(C) A;[lu 01 (C-C-)] dC 1 

- I* e(c') J\ [1 u Qi (c'-c)] dc 1 } (7.11a) 

C 2 u 

A combination of Eqs. (3.23), (5.2a), and (7.10) gives an alternate form of 
the energy equation for the steady case as 
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6' - 2(3 5 2 - 2*3) + 1 
n 

(L/k)^ Hjj (/^ (de/de*) V~ U 01 (H')] d*' 

+ \ (d0/dr) ^ [jU oi (r-€)J d*'} (7.1 

Note that this equation can be obtained directly by integrating the left- 

hand side of Eq. (7.11a) by parts. Equations (7.11) provide two forms of 
the energy equation for the steady-state conditions. 

For the case of negligible radiation, Eqs. (7.11) reduce to 


O' = 2 (3e* - 2*3) - 1 


The solution of Eqs. (7.12) is found to be 


(7.12) 


0(0 - * (2* 2 -*3-i) (713 

Thus, a combination of Eqs. (7.7) and (7.13) gives the result for the bulk 
temperature for the steady case with no radiation as 


- 0 b 3 17/70 (7.14) 

This result Is useful in determining the extent of radiative contributions. 
Transie nt Radiative Interactions 

For the transient case, a combination of Eqs. (3.22), (5.2a), and (7.5) 
gives the energy equation for the linearized radiation as 
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e cc - e T - 3Ne - 12 U-e 2 ) 

’ f (L/k> |il f o 9U '’ T) [ 7 u 01 (£ - {,) 1 dC 

+ f l 0 (€ ’ , t) A" u Qi (c-e)] de'l (7.15a) 

where 

n 

N = (PL2/k) K. = (PL 2 /k) E S,(T) (de /dT) T 

1=1 1 “i T i 

Note that this definition of N Is slightly different than the definition of 
Nj in Eq. (5.4c). The dimensionless gas property N characterizes the rela- 
tive importance of radiation versus conduction within the gas under opti- 
cally thin conditions [23, 41]. Also, by combining Eqs. (3.24), (5.2a), and 
(7.5) another form of the transient energy equation is obtained as 

e ?e - e T - 12 U-e*) 

= 2 (L/k) 1-1 HliUoi * A 1 u 0 i (c " 5 ‘ )] d£ ' 

- (ae/ac • ) Aj u Qi (c'-e)] dc * 1 (7.15b) 

Note again that Eq. (7.15b) can be obtained directly by integrating the 
left-hand side of Eq. (7.15a) by parts. Quite often, Eq. (7.15b) is the 
convenient form to use In radiative transfer analyses. 

For the case of negligible radiation, N = 0 and both forms of Eq. 

(7.15) reduce to 
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e 5C - e T = 12 u-s 2 ) (7.16) 

By employing the product solution procedure, the solution of Eq. (7.16) can 
be obtained and the result can be expressed in terms of the bulk temperature 
through use of Eq. (7.7). 

The solution of Eq. (7.16) is assumed to be of the form 

0(C,t) = g(e) + h(c.r) (7.17) 

From Eqs. (7.16) and (7.17), there Is obtained two separate equations as 

g" = 12(C-€ 2 ) (7.18) 

h^-h T = 0 (7.19) 

The solution of Eq. (7.18) Is obtained by direct integration as 

gU) * C(2C 2 -C 3 -1) (7.20) 

This Is the same result as given by Eq. (7.13) for the steady case if gU) 

Is replaced by 6(c). The solution of Eq. (7.19) Is found to be (see Sec. 
7.2) 


n 

h(c,r) * E C sin (nwc) exp[-(n») 2 T] (7.21a) 

n=l n 


where 


42 


C — 9 r ^ 

n " " 2 \ sin(mre) d£, n = 1,2,... 

^us, the complete solution of Eq ( 7 i fi \ 1e « 

cq * I'-lo; Is given by 


(7.21b) 


»(e, r) = 5 (2 c 2 -?3 


+ n ij C n s1n W exp(-a 2 T ) ; 


3 = m r 


(7.22) 


expression for C is oht>in^ , 

n Ulned frm E ’ s - (7-?0) and (7.216) as 

C " ' <4/aS) C1? ' “■* * a ‘> c„ S <a, - 24J , „ . 1>2 p 

g p> (7.22). B, Eqs (7 7, and (7 

expression for tde 6o, k tanperatore „ obta)ned „ ' <7 ' 2?) ' *' 


23) 


b / ° 6 n =l ° n ^ 1/a ^ + H/a 3 )] exp (- a 2 T ) ( 7> 

Where C n is given by Eq . (7 . 23)< 

7,2 Optically Thin Limit 


24) 


M.i «/>, (;25a) 

6' - 2 (3f2-2r3)+ 1 a 3 ^ 

' 2 c/ 0 (c ' f,) ( de/d C')de' 
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+ f 1 (e'-C) (de/de'JdC' (7.25b) 

C 

The differentiation of Eqs. (7.25a) and (7.25b) yields the same energy 
equation for the optically thin limit as 

0 " - 3N0 = 12 (C-e 2 ) (7.26) 

The solution of Eq. (7.26) satisfying the boundary conditions 0(0, x) = 0 and 
0 (1,t ) = 0 is found to be 

0(C) = (16/3N 2 ) [s 1 nh(- V"3N/2)/s1nh (^3N)] cosh[V~3MC-l/2)] 

+ (4/N)(e 2 -C + 2/3N) (7.27a) 

Alternately, the solution of Eq. (7.26) is written as 

0(C) = C! exp (>[m C) + C 2 exp {-y[m 5) 

+ (1/m 2 ) (24-12 me + 12 m c 2 ) ; m * 3N 

The constants Ci and (4 are obtained by using the boundary conditions 0 (O) = 
0 and 0 * ( 1 / 2 ) a 0 , and the solution for 0 (c) is found to be 

0 (C) = (1/m 2 ) (-24/(1 +e’V^)] (e'VvV™" K + e'V* C ) 

+ 24 - 12 m C + 12 m C 2 } (7.27b) 

Equations (7.27a) and (7.27b) should produce Identical results. The 
expression for the bulk temperature, in this case, is obtained by combining 
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Eqs. (7.7) and (7.27b) as 

0 b a ^77 r [(l-e’V"S/(l+e“V™) J - + ?i . li 

m 2 iii 3 m2 5m 

or 

0 b • [l/(3N) 3 j (576(3N)" 1/2 (NEXP) - 21.6N2 + 72N - 288} 

where 

NEXP = f 1-exp [-(3N) 1/2 ]}/fl + exp[-(3N) 1/2 ]} 

In the optically thin limit, the transient energy equations, 
(7.15a) and (7.15b), reduce to 

e «~ 0 T ‘ ^ = 12 

0 €e - 0 T - 12(€ - e 2 ) 

= ~ N C/q ( 30 / 0 C')dC- /*( 30/3C) dc'] 

Note that Eq. (7.29b) Is identical to Eq. (7.29a). 

The solution of Eq. (7.29) Is assumed to be of the form 

0 (5 » T ) 3 9(0 + h (c,t) 

Thus, Eq. (7.29) can be written as 

- h T - 3Nh » - g^ + 3Ng + 12(C - C 2 ) 


(7.28a) 

(7.28b) 

Eqs . 

(7.29a) 

(7.29b) 

(7.30) 


Consequently, 


(7.31) 
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g" - 3Ng = 12(c - C 2 ) ( 7 . 32 ) 

and 

h CC - h T ' 3Nh = 0. (7.33) 

The conditions for Eqs. (7.32) and (7.33) are obtained from Eq. (7.6) as 

® (°» T ) 3 h (°» T ) + 9(0) 3 0; h(0,r) = 0, g(0) = 0 (7.34a) 

9(l,r) = h(l, T ) + g(l) 3 0; h(l, T ) = 0, g(l) = 0 (7.34b) 

, 0(5,0) = h(c ,0) + g(C) 3 0; h(c ,0) = - g(c) (7.34c) 

The solution of Eq. (7.32) satisfying the boundary conditions given by 
Eqs. (7.34a) and (7.34b) is identical to the solution of Eq. (7.26) as given 
by Eq. (7.27) if 9(c) is replaced by g(e), i.e., 

g(5) = (16/3N 2 ) [sinh (-^3N/2)/ sinh(^llN)] cosh [ ^3N (c-1/2)] 

+ (4/N > (5 2 - 5 + 2/3N ) (7.35) 

The solution of Eq. (7.33) is obtained by using the product solution 
procedure and implying the conditions h(0,t) = h(l, r) = 0 and h( c ,0) = 

-g(€). For the product solution, it Is assuned that 

h(5,r) 3 F(C) G(t) (7.36) 

By using Eq. (7.36), Eq. (7.33) is separated into two ordinary differential 
equations which are expressed along with appropriate conditions as 


F" + X 2 F = 0 ; F ( 0) 3 0, F(l) = 0 


(7.37) 
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6 + (3N + > 2 ) G = 0 ; h(C ,0) = G(0) = gU) 

The solution of Eq. (7.37) Is given by 

F n (0 * sin (n»e), n = 1,2,... 

and the solution of Eq. (7.38) Is found to be 

G n (r) = C n exp {- [3N ♦ ( nir ) 2 ] t} 

Thus, the complete solution of Eq. (7.33) Is 

00 

h(c,r) * E C sin (nir?) exp f - [3N + (nir ) 2 ]t} 
n=l n 


where 


(7.38) 


(7.39) 


C n - -2 /* g(C ) sin (nir?) de, n = 1,2,... (7.40) 

n o 

Now, the solution of Eq. (7.29), as expressed by Eq. (7.30), Is written as 


eU.r) = ( 16/3N 2 ) [slnh (-V3N/2)/ slnh(VlN)] cosh [ ^3fT U-l/2)] 

+ (4/N)(5 2 - C + 2/3N) 

00 

+ E C n sin (mr$) exp{-[3N + (nw ) 2 ]r) (7.41) 


From Eqs. (7.35) and (7.40), It follows that 
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n » for n even 

= 32[3N+(nn)2]/[3N2(nw)3] +2(n* )/f 3N2[3N+(n,r)2]} , f or n odd 

8y combining Eqs. (7.7) and (7.41), the expression for the bulk 
temperature Is obtained as 

°b = 6 ( (16/3N 2 )[s1nh(-Vii/2)/s1nh(^lN)][(l/3N) cosh( yfsH/Z) 

- (4 + V3N) (3N) -3/2 sinh( V"3N/2)J + (4/N)[-l/30 + 1/(9N)] 

OB 

+ r c n Cl/(mr) + 4/(nn)3] exp[-(3N + n 2 w 2 ) T ]1 (7 43) 

n=l ' * ' 

where C n Is given by Eq. (7.42). 

7.3 Large Path Length Limit 

In the large path length limit, the steady-state energy equations, Eqs. 
(7.11a) and (7.11b), reduce to 

0' - 2(3£ 2 - 2£ 3 ) + 1 = M / 1 0(C) de'/($ - £') (7.44a) 

0’ - 2 ( 3C 2 - 2? 3 ) + 1 ° 

n 

■ (L/k) ^ H lf |/' {de/dC ) m[iu 0) (e - r)]d{' 

♦/jfdd/df) tn [2 Uol(c . . c) ] dc . (7 ,44b) 

where 

n 

M = Hi L/k = (L/k) I ^(de^/dTJfj (7.44c) 


(7.42a) 

(7.42b) 


Through Integration by parts. It can be shown that Eq. (7.44b) reduces to 
Eq. (7.44a). The parameter M In Eq. (7.44c) Is defined differently than M, 
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in Eq. (5.4e). The nondimensional paraneter M constitutes the radiation- 
conduction interaction parameter for the large path length limit [23, 41]. 
Equation (7.44a) does not appear to possess a closed form solution; a 
numerican solution, however, can be obtained easily. 

In the large path length limit, the transient energy equations, Eqs. 
(7.15a) and (7.15b), reduce to 

0 CC - « T - 3N0 - 12 (C - C 2 ) 

= - (Hjl/k) [ /Ms.t) - $')2 + /* eU',r) d£'/(e' - ^)2] (7.45a) 

e ec - e T - 12 (e - £ 2 ) = (HiL/k) /* (30/3C')de'/(? - S’) (7.45b) 

Since )(£-£ ' ) 2 | *| (I '-£ ) 2 | , Eq. (7.45a) can be written as 
e " 6 t ‘ 3Nfl ' 12 ^ ‘ ^ 

•-(HiL/k) I* 9(C',t) dC'/(C - C') 2 (7.45c) 

0 

Through integration by parts, Eq. (7.45c) can be expressed as 
9 rr - 6 - 3N0 - 12 (e - e 2 ) * (Hj L/k) /* (36/H') dC'/(e - S') (7.45d) 

ts T o 

Equations (7.45a) - (7.45d) represent different forms of the governing equa- 
tions In the large path length limit. With the exception of the term (-3N6) 
on the left-hand side, Eq. (7.45d) Is identical to Eq. (7.45b). Since N 
represents the radiation-conduction interaction parameter only In the 
optically thin limit [23], It should not appear In the governing equation 
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for the large path length limit. Thus, Eq. (7.45b) Is the correct equation 
to use for solution In the large path length limit; the solution of this 
equation is obtained by numerical techniques. 

7.4 Method of Solution 

The solution procedures for both steady and unsteady cases are pre- 
sented in this section. In principle, the same nunerical procedure applies 
to both the general and large path length limit cases. 

7.4.1 Steady- State Solutions 

The general solution of Eq. (7.11a) or Eq. (7.11b) is obtained numer- 
ically by employing the method of variation of parameters. For this, a 
polynomial form for 0(c) is assumed in powers of c as 

0(C) = r a C m (7.46) 

m=0 

By considering a five term series solution (a quartic solution in c) and 
satisfying the boundary conditions 0(0)= 0'(l/2) = 0 and 0'(O) = -0'(1), 
one obtains 


0(C) - a,(C - 2c 3 ♦ C 4 ) + a 2 (C 2 - 2* 3 + £ 4 ) (7.47) 

Thus, 

0 ' (C ) ■ a,(l - 6 C 2 + 4c 3 ) + a 2 (2c - 6c* + 4c 3 ) (7.48) 

A substitution of Eq. (7.48) in Eq. (7.11a) results in 


ai(l - 6c 2 + 4c 3 ) + a 2 (# - 6* 2 + 4* 3 ) - 2(3* 2 - 2* 3 ) + 1 
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= _( L /k ) ^ H llUoi {/« e(€«) A|[-u 01 (e -«•>]*• 

- /W) ^ [- U 01 (c* - e)] d£'} (7.49) 

£ 2 

where expressions for e(£') are obtained from Eq. (7.47). 

The two unknown constants a } and a 2 in Eq. (7.49) are evaluated by 
satisfying the integral equation at two convenient locations (£=0 and £ * 
1/4 In the present case). The entire procedure for obtaining ajand a 2 is 
described In Appendix E from which it follows that 

a! * (llo 2 - 16 aJ/DEN 
aj » (16 a3 - ll ai )/DEN 

where 

DEN = 16 (<*i<i4 - a 2 aj) 

and coefficients ai through 04 are defined In Appendix E. 

Now, with known values of a 3 and a 2 , Eq. (7.47) provides the general 
solution for 6 (£). The expression for the bulk temperature is obtained by 
combining Eqs. (7.7) and (7.47) as 

e b = (1/70) (17 ai + 3a 2 ) (7.51) 

Note that for the case of no radiative interaction « 2 , as, and ai* are equal 
to zero and a t = 1. Thus, % * 0 and a = -1, and Eq. (7.51) gives the 
result of Eq. (7.14). 

The governing equation for the large path length limit is Eq. (7.44a). 


(7.50a) 

(7.50b) 

(7.50c) 
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For this equation also the solution is given by Eq. (7.47) but the values of 
a's are completely different in this case. There are two approaches to 
obtain solutions in the large path length limit. One approach is to make 
use of Eq. (7.44a) and go through the entire numerical procedure described 
in Appendix E for the general solution. Another approach is to work with 
the general solution but evaluate all R. and integrals of Appendix E in 
the large path limit. In the large path length limit, the integrals can be 
evaluated in closed forms. Both procedures are described briefly in Appen- 
dix E. Irt order to distinguish the large path length limit solution from 
the general solution, constants a, and a 2 are replaced by bj and b 2 , and 
coefficients a , through a4 are replaced by 0, through 0 4 . The solution for 
the large path length limit, therefore. Is given by 

e(t) - bjU-253+C 4 ) + bs (C 2 -S 3 * 4 ) (7.52a) 

where 

bj = (Ub 2 - 16fi 4 ) /BOTTOM (7.52b) 

b 2 » (1603 - 1 16 i ) /BOTTOM (7.52c) 

BOTTOM = 16(8, 0 4 - 0 2 8 3 ) (7.52d) 

and coefficients Bi through 6 4 are defined in Appendix E. For this case, 
the bulk temperature is given by 

e b ■ (1/70) (17b, + 3b 2 ) (7.53) 

Note that in this case the value of coeffients b, through b 4 are obtained in 
closed form. 
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7 - 4 -2 Transient Solutions 

n>e governing energy equations for the transient case are Eqs. (7.15a) 
and (7,15b). As in Sec. 6.3, general solutions of these equations are 
obtained numerically by employing the method of variation of parameters. 

For the present problem, a polynomial form for e({,t) is assured as 

,<E ,T> *m'o a ” <T) 5 " (7.54) 

For a quadratic temperature distribution In E (with time dependent coeffl- 
dents), Eq. (7.54) Is written as 

e(Str) = ao( T ) + a^t) £ + a 2 (T) K 2 (7.55a) 

By using the boundary conditions 0(0, r) * 0 and 0 (c=l/2) = 0, this reduces 
to 

«U,t) * g(t) (S- E 2 ) (7.55b) 

where g(t) represents the time dependent coefficient. Consequently. 

V«.t> ’ 9(t> (1-2C); 9 {? ( 5 , t ) . -?,(,); 9 t ( 5 ,,) . ( E . { t )g . (t) ( ;. 56) 

Also, a combination of Eq. (7.6a) and (7.55b) yields the Initial condition 

®(S,0) * g(0) a 0 (7.57) 


Note that essential boundary conditions are used already In obtaining the 
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solution represented by Eq. (7.55b). 

By employing Eqs. (7.55b) and (7.56), Eqs. (7.15a) and (7.15b) are 
transformed in alternate forms which are expressed in a compact form as 


g'(T) + 


Ji (£ ) » 

J 2 (e) . 


g( t ) + 12 * 0 


(7.58) 


where Ji(c) and J 2 (c) are defined in Appendix F. The function Jj (£ ) is used 
for solution of Eq. (7.15a) and J 2 U) Is used for solution of Eq. (7.15b). 
The solution of Eq. (7.58) satisfying the initial conditions of Eq. (7.57) 
is given by 

g(r) = — {exp [-J(C )t ]- l} (7.59) 

J(C) 

The temperature distribution given by Eq. (7.56b) can be expressed now as 


e(C ,t) = — {exp C-0 (c)tJ-1}(c-C 2 ) (7.60) 

J(C) 

The expression for the bulk temperature Is obtained through use of Eq. (7.7) 
as 


= 72 / 1 [(C-C 2 ) 2 /J(C)]{exp[-J(5)T]-l} (7.61) 

0 o 

Note that In Eqs. ( 7 . 59 ) - ( 7 .61 ) . J(e) becomes Ji(£) for solution of Eq. 
(15a) and J 2 (c) for solution of Eq. (15b). 

For a quartlc solution in Eq. (7.54) gives the result Identical to 


54 


Eq. (7.47) which for the transient case is expressed as 

o (C ,T ) * g(r)(5 - 25 3 + 5 4 ) + h(r)(5 2 - 25 3 + 5 4 ) (7.62) 

Thus, 

e ? (C »t ) = (1 - 65 2 + 4c 3 )g(r) + 2(5 - 35 2 + 25 3 )h(r) (7.63a) 

9 ce (5,r) * 12( - 5 + 5 2 ) g(r) ♦ 2(1 - 65 + 65 2 ) h(r) (7.63b) 

9 t (5,t) 3 (5 - 25 3 + 5 4 ) g' (t) + (5 2 - 25 3 + 5 4 )h'(r) (7.63c) 

By substituting Eqs. (7.62) and (7.63) Into Eq. (7.15a), one obtains 

X g'(r) + J 3 (5 ) g(t) + yh' + Jt, (5 ) h(t) = - z (7.64) 

where 

X 3 (5 - 25 3 + 5 4 ); y 2 (5 2 - 25 3 + 5 4 ); z = 12(5 - 5 2 ) 

and functions ^( 5 ) and Jk*(?) are defined in Appendix F. Equation (7.64) 
constitutes one equation In two unknowns, namely g(r) and h(x). However, 
since the equation Is linear In t, the principle of superposition can be 
used to split the solution into two solutions as 

x g' (r ) + J 3 (5) g(r) = -z/2 (7.65) 

yh’ (t) + J 4 U) h( T ) = -z/2 (7.66) 

The Initial condition for this case can be written as 

0(5,0) * g(0) (5 - 25 + 5 3 ) ♦ h(0)(5 2 - 25 3 + 5 4 ) 3 0 (7.67a) 

Consequently, 


g(0) 3 0; h(0) 3 0 


(7.67b) 
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The solution of Eqs. (7.65) and (7.66) satisfying the appropriate 
initial condition of Eq. (7.67b) is given respectively as 

9(r) • [z(c)/2J 3 (c)] ( exp[ -J 3 (C )t / x(C ) ]-l} (7.68) 

h(r) = [z(C)/ 2 A,( ? )] { expC -J 4 (c ) T /y(c ) J-l} (7.69) 

By substituting Eqs. (7.68) and (7.69) into Eq. (7.62), the expression for 
the temperature distribution Is obtained as 

0(C»t) = [6(c-C 2 )(C-2c 3 +c 4 )/J 3 (c)]{exp[-J3(5)T/x(C )]-l} 

+ [ 6 (e-c 2 )(c 2 -2e 3 +* 4 )/J 4 ( 5 )]{exp[-j 4 (e) T /y(e) J-i) (7. 70) 


The bulk temperature in this case is given by 
l 

®b = 36 / 0 He - C 2 )(« - 2C 3 + C 4 )/J 3 (c)]fexp[-J 3 (C)T/x(C)] -l}dC 

+ 36 /^[U - ? 2 )(c 2 . 2 5 3 + c 4 )/Jk.(c)]{exp[-0,(€)T/y(C)]-l}dC (7.71) 

where x and y are defined in Eq. (7.64). 

By substituting Eqs. (7.62) and (7.63) into Eq. (7.15b), there is 
obtained 


*g' + JsU) g(r) + yh* + J 6 (C) h(t) = -z (7.72) 

where again x,y,z are defined in Eq. (7.64) and functions J 5 (c) and J 6 (C) 
are defined in Appendix F. The solution procedure for this equation is 
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identical to that for Eq. (7.64) and the results for temperature distri- 
bution and bulk temperature are given respectively by Eqs. (7.70) and (7.71) 
with J 3 replaced by Jg and J 4 by J 6 . 

In the large path length limit, the two applicable governing equations 
are Eqs. (7.45b) and (7.45d). The solutions of these equations can be 
obtained from the general solutions by evaluating the Integrals in J 
function In the large path length limit. 

Alternately, for a quadratic temperature distribution, Eqs. (7.45d) and 
(7.45b) are transformed respectively to 


g' (x) + 



g(-r) + 12 = 0 


(7.73) 


where J 7 (c ) and Je(S) are defined in Appendix F. The solution of Eq. 
(7.73) is given by Eq. (7.59) and expressions for e(c,r) and can be 

obtained from Eqs. (7.60) and (7.61) respectively. Of course, proper care 
should be taken to use the correct relation for J functions for different 
equations. The large path length limit solutions for a quartlc temperature 
distribution can be obtained in a similar manner. 
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8. PLANS FOR SPECIFIC RESULTS 


Some specific results have been obtained and these are being analyzed. 
The present plans are to obtain extensive results for the following cases with 
varying physical and flow conditions: 


A. Physical Geometries 

1. Parallel Plates: One-Dimensional Radiation 

2. Parallel Plates: Two-DImensonnal Radiation 

3. Diffusing Channel Flow: One- and Two-Dimensional Radiation 

4. Channel Flow: Top Plate Flat, Bottom Plate with a 5-15* ramp 

5. Scramjet Inlet Configurations 

B. Radiative Interaction Cases 


1. Transient Radiative Transfer In Homogeneous Gases 

2. Transient Energy Transfer By Radiation and Conduction In Homogeneous 
Gaseous Systems 

3. Transient Energy Transfer By Radiation, Conduction, and Convection In 
Homogeneous Gaseous Systems 

4. Applications to Flow of Homogeneous Gaseous Mixture 

5. Applications to Flow of Chemically Reacting Gaseous Mixtures 

6. Applications to the Scramjet Inlet Configurations. 

C. Boundary Conditions 

1. Isothermal Black Boundaries 

2. Isothermal Gray Boundaries 

3. Nonlsothermal Boundaries 

4. Boundaries with Uniform Heat Flux 

5. Actual Scramjet Inlet Configurations 

D. Flow Conditions 

1. Incompressible - Various Cases 

2. Compressible - M m * 0.5, 1, 2, 3, 4, and 5 

3. Temperature Range - 300, 500, 1000, 2000, and 5000 K 

4. Pressure Range - 0.1, 1, 2, 5, and 10 atm. 

5. Realistic Conditions for Scramjet Inlet Flows 

E. Participating Medium 

1. CO - One Fundamental Band 

2. CO? - Three Important Bands 

3. H?0 - Five Important Bands 

4. OH - One Fundamental Band 

5. CO 2 + HoO (Different Concentrations) 

6. OH + H 2 O (Different Concentrations) 

7. OH + H? + O 2 + H 2 O (Different Concentrations) 

8. A Realistic Combustion Model for the Hydrogen Air System 
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F. Specific Results 

1. Nongray Solutions Based on Band-Model Correlations 

2. Optically Thin Solutions 

3. Large Path Length Limit Solutions 

4. Gray Solutions 

5. Nongray Solutions for Scramjet Inlet Flows 


The procedure developed In this study has been applied to several realistic 
problems. This has resulted In various publications in forms of technical reports, 
technical papers (presented at national and international conferences) and journal 
articles. Results of selected studies are presented In appendices G through J. 
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9. CONCLUDING REMARKS 

A brief review Is presented on various band models and band model 
correlations that are useful In nongray radiative transfer analyses. 
Different formulations for one-dimensional radiative flux are provided. These 
are used to develop the basic governing equations for transient energy 
transfer In gaseous systems. Limiting forms of these equations are obtained 
In the optically thin and large path length limits. Numerical procedures are 
described to solve the governing equations for different physical and flow 
conditions. The plans for obtaining extensive results for different cases are 
provided. The formulation and numerical procedure presented In this study can 
be extended easily to multi-dimensional analyses. In the near future, the 
Influence of radiative Interactions will be Investigated for the realistic 
flow of hydrogen-air mixture In the scramjet Inlet configuration. 
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APPENDIX A 

EXPONENTIAL INTEGRALS AND EXPRESSIONS FOR RADIATIVE FLUX 
Some important relations for the exponential integrals are given 
Appendix 8 of Ref. 8; and it is noted that 


/ E n (t) dt * - E n + ! (t); E n (0) = l/(n-l) . 

Now, consider the first integral in Eq. (3.1) as 

id). /' J e, (t,.t)n 


By defining x * t - t, dt = - dx, and Eq. (A. 2) becomes 

KU*/° E 2 (x)(-dx) - /* X E 2 (x)dx 
T \ ° 

" - [6s(x)]„ * Ej(0) - Ej (t,) . 

Thus, 

E 3 (t x ) * E 3 (0) - 1(1) = j-/q X E 2 (r x -t)dt 
The second integral in Eq. (3.2) is written as 


1 n 


(A.l) 


(A. 2) 


(A. 3) 


1(2) - L 0X E 2 (t- T2 )dt . 
T X 


(A. 4) 
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By defining x = t - Eq . (A. 4) Is expressed as 

1(2) - / 0 °* * E 2 (x)dx = (Ej(x)I 0 0> * 

* - [E3(x 0 »- t ») - E *«»] 

or 

E 3 ( t 0X -t x ) = M°> - K2) - i- /^° X E 2 (t-T x )dt . (A. 5) 

A substitution of Eqs. (A. 3) and (A. 5) Into Eq. (3.1) results in (for 
the case B lX = ejx and B 2X = e 2 x) 

q RX = e U - e2 x 

+ 2 {/ Q X fe bx (t) - ejX] E 2 (r x -t) dt} 

‘ 2 f e bx (t) ‘ e 2 x]E 2 (t-T x )dt} . (A. 6) 

This equation when expressed In terms of the wave number us Is exactly the 
same as Eq. (3.4). Following a similar procedure, Eq. (3.5) is obtained 
from Eq. (3.2). 
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APPENDIX B 

ALTERNATE FORMS OF RADIATIVE FLUX EQUATIONS 
In radiative formulations. It Is desirable to express the relations for 
Q r and dlv In terms of A and A , and avoid the use of 7 . This 
Is accomplished by expressing the integrals containing A* and f in 
alternate forms through the procedure of Integration by parts. This Is 
performed by using the relation 


/a m dn = (mn)|j - /J n dm 


(B.l) 


Consider now the first Integral in Eq. (3.21) and 


express as 


8I(1) * il F i» (*') A 'r- “ oi («-oj « 


For Integration by parts, let 


( 8 . 2 ) 


m 3 F i“i dn 3 * f- u (C-C')J dc 


Then, 


dm ■ (dF lW /d£') dc' * (de/dC) d?' . 


In order to get n, let 


u * 7 u o du /de' = - 3 u J2 


and 


dn - [dA(u)/du] , fdA(u)/de ' x dj'/du) d{' 

■ - ( 2/3u 0 > [d A(u)/d£ 'J d{' . (. 2/3u 0 )dA(u) . 
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Thus, 

n = -(2/3 u q ) A(u) = -(2/3u 0 )A [| u Q (e-C ' )] . 


Consequently Eq. (B.2) can be written as 


BI(1) - F iu (c)r- — *(0)] - F 1U (0) \-J- J f|u Q 0] 
3u o 3u o 2 


+ — fi fde (C')/d5'] \ % U-e')]de’ 

3u Q 0 “ 


(B.3) 


Note that by the definitions given In Eqs. (2.2) and (2.8), A(0) = 0. In 
the present case, only the definition given in Eq. (2.8) Is acceptable for 
A(0) * 0. Also by definition, Fju>(0) = 0. Thus, Eq. (B.3) reduces to 


BI(1) = (2/3u q ) i\ [de aj (C , )/dC] A [1 u 0 U-e')]d5' 


(B.4) 


The second integral in Eq. (3.21) is written as 


81(2) = (2/3u q ) fl rde M (c') A \l u Q U'-OJdC' 


Let, 


(B.5) 


m * F 2U ; dn = u (c'-c)]dc' 
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Then, 


dm = (dF 2u /dS')d 5 ' = (de /dc')d€' 


Now, to get n, let 


and 


Thus, 


u 3 ± u 0 (C'-C), du/dC’ 3 3u q /2 


dn = [dA(u)/du]d£ ' * [dA(u)/dC* x dc'/dujd?' 


= (2/3u 0 )[dA(u)/dC'j d?' = (2/3u q ) dA(u) 


n = (2/3 u q )A(u) =(2/3u 0 )Afl u Q (e 1 -C)] 


Consequently, Eq. (B.5) is expressed as 


81(2) - F 2 «(l) {(2/3u q ) A f 1 u q (1 -5 )] } -{ F 2 o»(c ) [(2/3u q ) A(0)]} 
- (2/3u q ) /* [de aj (C , )/dC] A [- u Q U'-0] dC 


(B.6) 


Since, by definition, Fp u (l) 3 0 and A(0) = 0, Eq. (8.6) reduces to 



By use of Eqs. (B.4) and (B.7), the Integrals in Eq. (3.21) can be 
expressed in alternate forms and this results in Eq. (3.23). 

Consider now the integrals in Eq. (3.22); the first integral is written 
as 

8I(3) " /o F lo, (r) ^ f | U 0 ‘ (B - 8) 

For integration by parts, let 

m 3 Fiu U'); dn * I (- u (?-?')] dc’ 

2 u 

Then 

dm = (dF lu /d c ') d 5 V = (de^/dc 1 ) d e • 

As before, to obtain n, let 

u 3 j u 0 U-5'). du/dC = - 3u q /2 


and, therefore. 


70 


dn = — f d7T(u)/du] d£ ' 
du 

= {— — fdA(u)/du] de ' 
dt' du 

3 df (-2/3 u 0 )fdA(u)/du]} 

Thus, 

n = (-2/3u q ) A (u) = (-2/3 u Q ) A [i u Q (c-C)] 

and Eq. (B.8) can be written as 

BI(3) = {F lU (-2/3 u Q ) A fl u Q (C-C* )]}J 

- /o H-2/3 u q ) A'fl u 0 (c-{')J }[dFju(C 1 )/d€ ' ] dC 

Since F^u> ( 0 ) = 0 and X (0) = 1, Eq . ( 8.9) reduces to 

8 1(3) = ( 2/3u 0 ) f- F jw(?) + /q f d ( € ' )/d£ ' ] A'fl u Q U-?')]de'} 

Similarly, for the second Integral In Eq. (3.22), one can find 

BI(4) = /* F 2 „ ( € ') a'[1 u q (c'-c)] dc ' 

* ( " 2/3u 0 )^ F 2 (€) + /* (de u U')/dc'] u q (€*-c)]dc‘f 


(B.9) 


( B .10) 


(« - 11 ) 
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A substitution of Eqs. (B.10) and (B.ll) in Eq. (3.22) results in Eq . 
(3.24). Also, a differentiation of Eq. (3.23) with respect to e, by using 
the Leibnitz formula, gives 


dq (C) n 

*1 A oi f 
d? 1 = 1 


de , (c) 


"A(O) - 0 


3u , C de ,(?') , , 

-21 / (—2! — I. [!«., ({-?■)! <*' 

2 „ dr' 1 2 01 




3u . ! de .(?') , , 

- -2i/ t A-,fiu t (f]de'l 

2 c dc' 1 2 01 

Since A(0) =0, Eq. (B.12) reduces to Eq . (3.24). 


( 8 . 12 ) 
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APPENDIX C 

DEFINITION AND EVALUATION OF INTEGRAL FUNCTIONS 
For the convenience and use In the computational procedure, the follow- 
ing definitions are employed In expressing the relations for the integral 
functions: 


b i * 

c i ■ 1/b i - 

c(t) * l/<{-{ 2 ) 

r - 1 /(ot{) 


Various integrals are defined and simplified as follows: 


G i(?) * C(e) ( 3Nj (C -€ 2 ) + 

*7 I 1 ", M ii u oi 1/5 V fb, (e-e')J «- 

+ (c'-r 2 ) fb^e'-o] dc 1 }) 


- C(C) ( 3Nj (5-C 2 ) + 

3 " b i 5 

+ 2 i*l Ml1 Uq1 ^- c i u -« 2+ 2CCiU-(Ciu)2] *(u)du 


♦/ 

o 


b^I-O 


fC+c i u-«*-2|c i u-(c i u) 2 ]-^(u)du}) 


(C.lb) 
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G 2 U) = CU) (3N,(I - *2) 

"T <h (M li U oi f/ o C ' 2 A i"[ b i(« , )l d 5* 

+ /* (C ,2 -l) a’’ [b 1 (C *-C )]d? '}) 

- C(e) ( 3Nj (1 - C 2) 

2 

- j T M L1 u Qi {/^ fc 2 -2c c lU + ( Cl u)2] (u)du 

Ml-0 

+ V U 2 + 2c c i u + (c^) 2 -1] 7^ (u) du } ) 

G 3 (€) = C(C) (i ^ H u \r 0 (i-2e • ) A i [^(e-oldc* 

- /* (1-2C') j\ fb 1 (c'-C)]dC’}) 

n b.c 

- C(c) ( T M H f/ 0 [1-25 + 2 Ci u] ^(u)du 

Ml-5) 

- / 0 f 1-25 - Cju] A^ (u)duf) 

n r 

a,({) »C(?) (-Jt lt„ (/'«' A, [b, (C -t ' )]d? • 

- /{ 5' *i [b, U'-C)]dE’)) 


(C .2a) 


C .2b) 


(C.3a) 


(C.3b) 


(C.4a) 
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* C(0 (- 2 l {/J [C-c 1 u] ^ (u) du 

1*1 


Mi-c) 

+ / C 1 fc + t^u] Aj(u)du}) 


(MO s (3T/8) ^ u Qi H u f/J (e'-C' 2 ) [b^ (5-5 ' )jd5' 


- /J (C'-€' 2 ) ^ fb^C-OldC'} 


n ^4^ 2 _ 

= (r/4) f/ 0 f e -c 1 u -c 2 + 2Cc 1 u-(c i u) ] A! (u) du 


b, (1-5) 2 , 

- / Q [c+c^u-5 2 - Zc^u _ (c^ u) ] ^(u) du} 


n t 2 - 

6g(5) = (3 r/8) J u Qi {/* 5' Aj (5-5 ' )]d5 ' 


+ /* (1 -€ * 2 ) Aj fb i (5'-5)] d5’} 


n M 

* (r/4) y H^{/ 0 ' fC 2 “2C c 1 u +(c 1 u) 2 ] A 1 (u) du 


Ml-5) 

+ f Q [l-5 2 -25 C< u - c iU ) 2 ] A^u) du} 


(C4.b) 


(C.5a) 


(C ,5b) 


( C .6a) 


(C.6b) 
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G,(5) ■ (r„) t H„ (/' (1-25') A, fb, (?-£•)] d«' 

-/; »-*') *, fb, ({'-?)] dC'| 
n b.£ 

= (r/6) t (H 11/ u oi ) {/ 0 1 (l-2c + 2c t u) tyu) du 

Ml-0 

- / 0 1 (1 -2c -2c 1 u) X| (u) du} 

n r 

Ge(0 = (r/4) ^ H X1 {/; («') A 1 [ b 1 (C -e ' ) d?’ 

- z 1 (25') x [bWe'-e)] dc 1 } 

c 

n b.e 

= (r/6) E (H 11 /u Ql )f/ 0 (25 - 2c^u) A 1 (u) du 

1= 1 

Ml-5) 

- J 0 (25 + 2c 1 u) Xj (u) du} 

n l 

Gg * (3T/8) * U(j1 Hl1 f/ Q (C'-C |2 ) A, (b^ 1 ) d5 ' } 

* (r/4) E (/ 0 1 [ u - ( Ci u)2J Xj (u) du} 


(C.7a) 

( C .7b) 
(C.8a) 

(C .8b) 
(C.9a) 


(C.9b) 
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Gio - (3178) l ml u 0 i H n {/ 0 (M 1 ) A| (b^ 1 ) 6V] 


* (r/4) f/ 0 1 (1 -c 1 u) (u) du) 


Gn * (r/4) ^ H u f/ Q (1-2C ' ) A 1 (b^C ’ ) dC 


- (r/6) ^ (H X1 ff 0 ' (1-2 Cl u) ^ (U) du) 


Gl 2 = (r/4) ^ (H^/U^) f/ Q (2C ■ ) A i (b 1 C) d€ 1 } 


= (r/6) (f 0 (2 C 1 u) A i (u) du) 


£ I 

Gi3 (?) * C(c) { 3Ni (c-?2) - M, [/ (c'- c' 2 ) -Jlr- 

(C-C ) 2 


+ / (?'-?' ) — — — ]} (Improper) 


U'-O*' 


G,4 (C) » C(?) { 3N, (_ - ?t) ♦ H, [/ C' 2 ------ 

2 0 (C-O 2 


♦/, (€' 2 - 1 ) 


(«■ - e) : 


(C.lOa) 


(C.lOb) 


(C .11a) 


(C.llb) 


(C .12a) 


( C . 1 2b) 


(C.13) 


(Improper) 


(C.14) 



77 


<hs (€) « M, C(C) /* (1-25 ») -AL 

(€-€') 


= Mi C(C) (2 + (25 -1) tn [(5 -l)/5 ]} 


G 1( («) ’ -2H, C( C )/‘ 

(5-5 ' ) 


2 M! C(5 ) { 1 + 5 5n [(5 -!)/£]} 


G , 7 (*) - (rHj/4) f 1 (? ' "5 ’ 2 ) — ^ 7 - 

U-e') 


= (rHj/4) {-i ♦* - (e - 5 2 ) tn CU-D/c]} 
2 


610 (c) * (rHj/4) f/* c '2 -Jil- / (Improper) 

(c-C) K U-C) 


G 19 


(5 ) = (r/4) l H t1 f/ 5 (1-25 ' ) tn [b^C-5 ■ )]«■ 

1*1 n 


(C .15a) 

(C .15b) 
( C .16a) 

(C .16b) 
( C . 1 7a) 

(C.17b) 

(C .18) 


+ /* (1-25 ’ ) tn [b i (5 '-5)] d5 ’} 

■ (rH,/4) U ! -() tn [»-{)/*]} 

n , 

G 20 (O 3 ( r /4) I H n { r Q 2C t n [b i (5-5 • )] dC 


(C .19a) 
(C.19b) 


a (r/4) ( Hj {^2 t n [ c /(l- c )] - % - l ~) 

+ I H,, in [b. (1-c)]) 

1=1 11 1 

G 21 58 (r^/4) / 0 (£'-£' 2 ) dC'/C’ 3 (rHi/8) 

G 22 - (rHi/4) f 0 (M' 2 ) dC'/C 

3 - (rHj/4) [A. ♦ £n (0)] (Improper) 

2 

n 

G 2 3 « ( r /4) J / Q (1-2C ' ) in (b 1 £') d£' 

3 - (r H! /8) = - Hj / (8 cjTi 3 ) 
n 1 

G 24 * ( r /4) l H u f 0 2£ ' tn (b i £') d£' 

* (r/4) [- ( H, / 2 ) + ? H li tn (b^J 


(C .20b) 

(C.21) 

(C .22a) 
( C -22b) 

(C.23a) 

(C .23b) 
(C .24a) 

( C .24b) 
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APPENDIX D 

INFORMATION FOR NUMERICAL PROCEDURE 

1. Data: 

T 3 300, 500, 1000, 2000, 5000, etc. ~ K 

P 3 0.1, 1, 2, 5, 10, etc. ~ atm 

L = 1, 2, 5, 10, 20, etc. ~ cm 

£ ♦ x = 0.0 ♦ 1.0; x = U( I) 


2. Thermal conductivity of the gas: 


K f = K f (T, P) 


3. The Planck function and its derivative: 


e tJ T > = 


exp (C 2 /T) - 1 


= PF 


de bu Cl C 2 exp ( C 2/ T ) 


dT T2 [exp (C 2 /T)-l) 2 


= PFD 


(D.l) 


( 0 . 2 ) 


(0.3) 


where 

Cj 3 (2Trhc2)ui 3 
C 2 3 (hc/k) <■> 

C 1 C 2 3 (2ir h2 cs/kjo) 1 * 

and u 3 a) = id 3 wave nunber center of the 1th band, 
c o 
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By defining 

TEXP = exp(C 2 /T) 

Eq. (D.3) Is expressed as 

PFO = (C 1 C 2 )(TEXP)/[(T 2 )(TEXP - l) 2 ] ( D.4) 

All values In Eqs. (D.2) - (D.4) should be evaluated for each band. 

4. Information or relation for A„ for each band, A„ = f(T). 

o ’o' 7 

5. Information on C 2 for each band, C 2 = f(T). 

o o 

6. Information on u Q for each band, u Q = C Q 2 PL (nondlmensional) . 

7. Information on 82(1) for each band (nondimensional) . 

8. Equivalent or effective pressure relation for each band. 

p e t = [( p / p 0 ) + (W ( b i ~ !)] n (nondimensional) . (D.5) 

In Eq. (D.5), P Q * 1 atm and, therefore, P^ and P must be In the units 
of atm. Note that P^ Is the partial pressure of the 1th species In a 
gases mixture and P Is the total pressure; b^ (the self broadening 
coefflcent) and n are different for different bands. For a single 
component system, Eq. (D.5) is usually expressed as 

p e ■ (b P/P Q ) n (nondimensional) (0.6) 


9. Line structure parameter for each band: 
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6 = B 2 P g = BETA (D. 7) 

10. Correlation for each band (for Tien and Lowder's correlation): 

f(B) = 2.94 [1 - exp(-2.606)] = F (0.8) 

11. Band absorptance correlation for each band (Tien and Lowder's 
correl ation) : 

7T (u ,6 ) = £n (uF ( U - 2 ■■■) + 1] = AU (0.9) 

u + 2F 

12. The derivative of the band absorptance correlation for each band (Tien 
and Lowder's correlation) 

A ' (u ,8 ) = fF (u 2 + 4 u F + 4F)]/0EN = AUD (0.10) 

where 

DEN = fF (u 2 + 2u + 2) + u](u + 2F) 

These basic relations are used in the governing equations of Section 6 


to obtain nunerlcal solutions for specific gaseous systems. The spectro- 
scopic and correlation quantities needed for these calculations are available 
In [22,24]. 
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APPENDIX E 

EVALUATION OF CONSTANTS FOR STEADY LAMINAR FLOWS 

To determine the constants in Eq. (7.47), Eq. (7.49) Is evaluated at 
5=0 and 5 = 1/4. To avoid excessive writing, the following notations are 
used (some of which are also used In Appendix C) 

£ ' = n , d£ ' = dn 

b i “ 3u 01^ 2; C 1 = ^ b 1 * 2/3u 01 

For £ =0, Eq. (7.49) reduces to 

■, n l 

a, + 1 + i (L/k) z H H u ni f / [aj (n-2n 3 +n 4 ) 

2 1*1 11 U1 o 

+ 32 („2-2 n 3V)] A*(b in )dr,J ■ 0 (E.l) 

By defining u = b^n, Eq. (E.l) Is expressed In an alternate form as 

a if l + (L/k) i H Xi /\c 1 u-2(c i u)3+(c 1 u)-]7Tj(u)du} 

+a 2 f(L/k) r H / ^(c u)2-2(c u)3+( c u)**] P(u) du}=-l (E.2) 

] 3^**0 * * ' ' 

Now, by defining the following Integral functions. 
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R n ■ V u i *; (u i> du t 

r 2 i ' C' u i ^ (u f’ du i 

r 3i ■ v u i (i v du i 

"«1 * ^ "i *i <V dU 1 

Eq. (E.2) is expressed as 


®1 a l + 82 aj = “1 


where 


■I'l* (L/k) ^ H u (c,R lt - 2 c? R 3i * cl R 4i ) 
-»* * (L/k) H n (c^R 21 - 2 c? R 3j » t | R 4i ) 


For C = 1/4, Eq. (7.49) reduces to 


H «, ♦ J. a, + li 

16 16 16 


7 n l/U 1 

1 (L/k) i H li u Qi {/^ 0(n ) J\ [b i (i - n)] dn 


• / 1/4 9(n) TTj [b^n - 1)] dn} 


(E .3a) 
(E .3b) 
(E.3c) 
(E.3d) 

(E.4a) 

(E.4b) 

(E.4c) 


(E.5) 
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By defining u * b^ (— - n) for the first integral and u = b. (n - — ) for the 
1 4 1 4 

second Integral, Eq. (E.5) is written as 


— + il a t + — aj 
16 16 16 


V?! 


(t/k) e H 1 . { / e( — c.u) A! (u) du 
1=1 11 « 4 1 1 


3b./4 . 

/ e(± + c.u) V. (u) du 

o a 1 1 


(E.6) 


By denoting d = c^u, the following relations are obtained from Eq. (7.47): 


and 


e(i- 

d) 

= M 

57 

-lid- 

9 

d 2 

+ d 3 

+ d 1 *) 

4 


256 

16 

8 





+ aj( 

9 

-lid - 

_L_ 

d2 

- d3 

+ d**) 



256 

16 

8 




e(I + 

d) 

= «,(■ 

57 

+ i!d- 

9 

d2 

- d 3 

+ d**) 

4 


256 

16 

8 





+ «2 

<-L 

.♦ J-d 

1 

. d2 

- d 3 

+ d 1 *) 




256 

16 

8 




( E . 7a) 


( E . 7b) 


A combination of Eqs. (E.6) and (E.7) results In 

b,/4 3b. /4 


a i(- E H f ^-[/ (u) du - / X’(u) du] 

K 1*1 11 256 o 1 J o 1 

uA! (u) du] 


11 ..Vi 


3b i /4 


— c,U uA!(u) du ■*■ / ... 

16 1 o 1 o 1 
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9 b 1 /4 

J --9 r f l 


3b. / 4 

c? [/ ' u 2 7T* (u) du - / u 2 A'(u) du] 

Q i o 1 o 1 


bj/4 


3^/4 


+ c 3 [ / u 3 7T'(u) du + / u 5 ft'(u) du] 

1o 1 0 1 


b./4 3b./4 n 

+ c 4 [ / 1 u 4 ff;(u) du - / 1 u-ff;(u) du]}- _) 


" 1 
n 


o 

b./4 


16 


3b. /4 


i II g ^ J / T 1 / T 

+ *2 {- Z H nf t / *i( u ) du - / TWu) du] 

K 1=1 111 256 o 1 ; o i 

, b./4 3b. /4 

- — c [ / u 7T!(u) du + / u 7T’ (u) du] 

16 1 o 1 o 1 

, b./4 3b. /4 

- — c 2 [ / u 2 fi!(u) du - / u 2 J!(u) du] 

8 1 o 1 o 1 

b,/4 3b. /4 

+ c 3 [ / u 3 7T’(u) du + / u 3 7T'(u) du] 

1 o 1 o 1 

b./4 3b. /4 , 

+ c 4 [ / 1 u 4 7T'(u) du - / u 4 ff.’(u) du]} - — •) 

1 o 1 o 1 16 


11/16 


8y noting that for any continuous function F(x) 


3/4 1/4 3/4 

/ F(x) dx * / F(x) dx + / F( x ) dx 

0 0 1/4 


(E.8) 


and defining 
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3b./4 _ 

S n - / Aj( Ul ) di^ (E.9a) 

b^/4 

b./4 

S 21 “ ( 0 U 1 W du l (E.9b) 

3b /4 

S 31 * “i *i<V dU i (E ' 9c) 

3b. /4 

S 4i = / u f *i'( u i ) du i (E.9d) 

b t /4 

b./4 3 

S 5i = L U 1 *)(“() du i (E.9e) 

»,'« 3 _ 

S 61 ' r 0 U 1 du 1 < E - 9f > 

3b./4 ^ 

* / u 1 X|(u.) du. (E.9g) 

b ^ /4 


Eq. (E. 8 ) can be written as 

a i 03 + a 2 a 4 = - 11/16 (E.lOa) 


where 
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n 

a 3 = (L/k) E H Xi f (57/256) + (11/16) c 1 (S 21 +S 3i ) 

- (9/8) S 4i - (S 5 i +S gi ) * o'; S 71 ] + 11/16 (E.lOb) 

n 

a 4 - (L/k) ^E (9/256) S u + (3/16) c 1 (S 2 ,+S 3i ) 

- (1/8) S 4i - c3 (S 5 i +S g1 ) + c*JS 71 1 + 3/16 (E.lOc) 

By solving Eqs. ( E .4 a) and (E.lOa) simultaneously, there is obtained 
the results for constants aj and a 2 as 

a! « (11 o 2 - 16 04 )/DEN ( E .1 la) 

aj = (16 o 3 - 11 ai )/OEN (E.llb) 

where 

DEN = 16 ( 0^04 ~ 02013 ) (E.llc) 

The governing equation for the large path length limit is Eq. (7.44a) 
for which the solution Is also given by Eq. (7.47). For the large path 
length limit, Eq. (7.47) Is expressed In the form of Eq. (7.52) which Is 
represented here as 

0(C) - b, (C - 2C 3 + C 4 ) ♦ b 2 (5 2 - 2C 3 + *") (E.12) 

Thus, 

0'(c) - b, (1 - 6 C 2 + 4 t 3) + b 2 (2c - 6 C 2 + 4 C 3) (E.13) 


\ 
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A substitution of Eqs. (E.12) and (E.13) In Eq. (7.44a) gives 

bj (1 - * 4?3) + b 2 (2e - 6c 2 + 4e 3 ) -2 (3C 2 - 2c 2 ) + 1 

* M (*»(£ ' ) d£ '/(£-£' ) (E.14) 

o 

For £ = 0, Eq. (E.14) reduces to 

b. + 1 « - M f 1 [#(£')/£'] d£' (E.15) 

o 

Upon substituting for 0(£’) from Eq. (E.12) Into Eq. (E.15), the Integrals 
can be evaluated In closed form and there Is obtained 

k 

bj Sj + b 2 b 2 * - 1 

where 

6, = 1 + (7/12) M; Bj - (1/12) N 

For £ * 1/4, Eq. (E.14) reduces to 

— b, +-^b 2 +^-*M f 1 e(£') d£ V(— - £') (E.17) 

16 16 16 o 4 

By substituting for e(£') In this equation, another relation between bj and 
b 2 can be obtained In terms of fl 3 and b 4 . But, this appears to Involve the 
evaluation of a few Improper Integrals. Thus, this approach Is abandoned In 


(E.16a) 

(E.16b) 
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favor of an alternate procedure discussed below. 

The solution in the large path length limit can be obtained from the 
general solution by evaluating the Integrals and in the limit of large 
path length. Since in this limit A ’ (u) = 1/u, the integrals in Eqs. (E.3) 
and (E.9) are evaluated to obtain 

R 11 s b 1 ; R 21 = b V 2; R 31 = b i /3; R 41 = b i /4 (E.18a) 

S u = in (3); S 21 = b 1 /4; S 31 = 3b i /4; S 41 = b*/4 ; 

S 5i = b3/192; S 61 = 9b3/64; S 7i = 5b**/64 (E.18b) 

From Eqs. (E.4) and (E.18a), there is obtained for the large path length 
limit 


bi 6 i + b 2 02 3 " 1 (E.19a) 

0! = 1 + (7/12) M; 02 * (1/12) M (E.19b) 

which is the same result as given by Eq. (E.16). From Eqs. (E.10) and 
(E.18b), one obtains In the large path length limit 

bj 03 + b2 6i* 3 - 11/16 (E.20a) 

where 

0 3 = 11/16 + M [(57/256) *n(3) + 65/192] 

* 11/16 + 0.583 154 559 M (E.20b) 
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B„ = 3/16 + M [(9/256) «n(3) + 17/192)] 
= 3/16 + 0.127 164 755 M 

The solution of Eqs. (E.19a) and (E.20a) yields 

bi =* (lle 2 - 16 6 4 )/B0TT0M 

b 2 = (1663 - 11 01 ) /BOTTOM 

where 

BOTTOM = 16 ( 0 1 8 i* “ B 2 0 3 ) 


With bi and b 2 known, the solution for the temperature distribution 
obtained from Eq. (E.12). 


(E.20c) 

(E.21a) 

(E.21b) 
(E .21c) 
is 
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APPENDIX F 

INTEGRAL FUNCTIONS FOR TRANSIENT LAMINAR FLOWS 

For convenience and use in the computational procedure, the following 
definitions are employed in expressing the relations for the integral 
functions: 


b 1 = 3 u q1 /2; c 1 » l/b 1 ; C(c) = 1/(C-E 2 ) 
x = C - 2E 3 + C 4 ; y = C 2 -2^ 3 + C 4 ; z = 12 (E-E 2 ) 

a(c) = 1 - 6C 2 + 4c 3 ; b(c) 3 2C - 6C 2 + 4C 3 


Various integrals are defined as follows: 


Jj(C) = 3 N + C(C) (2 + 1 (L/k) T Hii uL{/ (C *-E •*) A? [b i (E-E ' ) ] <*C ’ 

4 i=l Al Ul o 11 


+ / (E'-E* ) [b i (E ' -E ) JdE 1 }) 


(P.l) 


MO - C(c) (2 +1 (L/k) z H^u^l/ (1-2* 1 ) [b i (c-e t )]dc’ 

2 i 3 1 n 


- / (1-2C')^ [b 1 (c , -E)]dE , n 
E 


(F.2) 


Me) = z + 3Nx +1 (L/k) E H u2 {/ x( c ’) V' [b .(*-*')] dc’ 

4 1=1 1' oi ,J o 1 1 


I 
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+ / x(C) A? [b.U'-OJd?'} 

5 

9 n * _ 

J 4 U) 3 ( z - 2 ) + 3Ny + - (L/k) r H..U2J/ y (5 ' ) A" [b .(c-e ' ) Jd? 

4 1=1 11 U1 0 11 

1 

+ / y(e') a? [b.(^’)]dr} 

5 

1 n 5 - 

JsU) 3 2 + j (L/k) ^ H n u oi^ 0 a (*’) A i C b i^-e*)] d? 1 

1 

- / a(5 ') A.' [b i U '-5)]d5 '} 

5 


J 6 U) - (z - 2) + - (L/k) r H u If b(C) A! [b.U-e')] dc' 
2 1=1 11 01 0 11 

1 

- / t>u • ) a; [b f (e'-Olde *} 

€ 

1 

•MO * 3N + C(5) [2 + M / (1-25') d5'/(5-5')] 

0 

* 3N ♦ c(c) {2 + M [2 ♦ (25-1) tn 

5 


(F.3) 


(F.4) 


(F.5) 


(F .6) 


(F.7) 


Js(5) = *37 (^ ) - 3N 


(F.8) 
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Abstract 

A general formulation Is presented to 
Investigate the transient radiative Interaction 
In nongray absorbing-emitting species between two 
parallel plates. Depending on the desired 
sophistication and accuracy, any nongray absorp- 
tion model from the line-by-line models to the 
wlde-band model correlations can be employed In 
the formulation to Investigate the radiative 
Interaction. Special attention Is directed to 
Investigate the radiative Interaction In a system 
Initially at a uniform reference temperature and 
suddenly the temperature of the bottom plate Is 
reduced to a lower but constant temperature. The 
Interaction Is considered for the case of 
radiative equilibrium as well as for combined 
radiation and conduction. General as well as 
limiting forms of the governing equations are 
presented and solutions are obtained numerically 
by employing the method of variation of 

parameters. Specific resul ts are obtained for 
CO, C02» HgO, and OH. The Information on species 
H 2 O and OH Is of special Interest for the 

proposed scramjet engine application. The 
results demonstrate the relative ability of 
different species for radiative Interactions. 

Nomenclature 

A band absorptance * A(u,p), cm* 1 

A 0 band width parameter, cm “ 1 

C 0 correlation parameter, atm’ 1 - cm’ 1 

C- specific heat at constant pressure, 

kJ/kg-K * erg/gm-K 

e Planck's function, (W-cm’ z )/cm’ 1 

(1) 

e, t Planck's function evaluated at wave 

hi 

0 number u> 

0 

^ 1^2 emissive power of surfaces with 

temperatures Tj and T 2 * W-cm’ 2 
H I j , H I gas property for the large path length 

limit 

k thermal conductivity, erg/cm-sec-K 

K 11» K 1 9 as property for the optically thin 
limit 

L distance between plates 


^Eminent Professor, Department of Mechanical 
Engineering and Mechanics. AIAA Associate 
Fellow. 

Graduate Research Assistant, Department of 
Mechanical Englneerng and Mechanics. AIAA 
Student Member. 

++ Sen1or Research Scientist, Computational 
Methods Branch, High-Speed Aerodynamics 
Division. AIAA Associate Fellow. 


Mu, Mi large path length parameter, 
nond 1 mens Iona 1 

N 11* N 1 optically thin parameter, 

nond 1 men slonal 

N optically thin radiation-conduction 

parameter * Ni/R, nond 1 mens Iona! 

P pressure, atm 

q conduction plus radiation heat flux * 

Q c * 9R’ w/cm2 

Or total radiative heat flux, w/cm 2 

q c conduction heat flux, w/cm 2 

q Ru spectral radiation heat flux, 

(w-cm" 2 ) /cm" 1 

Q nondlmenslonal radiative heat flux 

Q nondlmenslonal conduction plus 

radiation heat flux 

R nondlmenslonal transient conduction 

parameter 

S Intearated Intensity of a wide band, 

atm* 1 -cm" 2 

t time, Sec (also used as t*) 

t^ characteristic time, sec 

t* nondlmenslonal time ■ t/t^ 

T temperature, K 

Ti # T 2 wall temperature, K; Ti ■ T w 

u nondlmenslonal coordinate * $Py/A 0 

u Q nondlmenslonal path length * SPL/A 0 

y transverse coordinate, cm 

9 nondlmenslonal temperature 

spectral absorption coefficient, cm" 1 
l nondlmenslonal coordinate * y/L ■ u/u Q 

P density, kg/m^ 

a Stefan-Bol tzmann constant, 

erg/(sec-cm 2 -K*) 
w wave number, cm" 1 

w Q wave number at the band center, cm" 1 

Introduction 

The field of radiative energy transfer In 
gaseous systems Is getting an ever Increasing 
attention recently because of Its applications In 
the areas of the earth's radiation budget studies 
and climate modeling, fire and combustion 
research, entry and reentry phenomena, hypersonic 
propulsion and defense-oriented research. In 
most studies Involving combined mass, momentum, 
and energy transfer, however, the radiative 

transfer formulation has been coupled mainly with 
the steady processes 1 " 11 and the Interaction of 
radiation In transient processes has received 
very little attention. Yet, the transient 
approach appears to be the logical way of 
formulating a problem In a general sense for 
elegant numerical and computational solutions. 
The steady-state solutions can be obtained as 
limiting solutions for large times. 

A few studies available on radiative 
Interactions reveal that the transient behavior of 
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a physical system can be Influenced significantly 
In the presence of radiation 12 " 17 . Lick Investi- 
gated the transient energy transfer by radiation 
and conduction through a seml-flnlte medium 12 . A 
kernel substitution technique was used to obtain 
analytic solutions and display the main features 
and parameters of the problem. Steady and 
transient heat transfer In conducting and radia- 
ting planar and cylindrical mediums were analyzed 
In Refs. 13 and 14 according to the differential 
formulation. The analyses based essentially on 
the gray formulation provide some qualitative 
Insight Into the effect of absorption and emis- 
sion on the transient temperature distribution In 
the gas. Ooornlnk and Herlng 1 ® studied the 
transient radiative transfer in a stationary 
plane layer of a nonconducting medium bounded by 
black walls. A rectangular Mllne-Eddlngton type 
relation was used to describe the frequency 
dependence of the absorption coefficient. It was 
found that the cooling of the layer Initially at 
a uniform temperature is strongly dependent on 
the absorption coefficient model employed. 
Larson and Vlskanta 16 Investigated the problem of 
transient combined laminar free convection and 
radiation in a rectangular enclosure. It was 
demonstrated that the radiation dominates the 
heat transfer In the enclosure and alters the 
convective flow patterns slgnlf Icantly. The 
transient heat exchange between a radiating plate 
and a high- temperature gas flow was Investigated 
by Melnikov and Sukhovlch 17 . Only the radiative 
Interaction from the plate was considered; the 
gas was treated as a non-participating medium. 
It was proved that the surface temperature Is a 
function of time and of longitudinal coordinate. 
Some other works on transient radiation and 
related areas are available In Refs. 18-22. 

The goal of this research Is to Include the 
nongray radiative formulation In the general 
unsteady governing equations and provide the 
step-by-step analysis and solution procedure for 
several realistic problems. The specific objec- 
tive of the present study Is to Investigate the 
Interaction of nongray radiation In transient 
transfer processes In a general sense. Atten- 
tion, however, will be directed first to a simple 
problem of the transient radiative exchange 
between two parallel plates. In subsequent 
studies, the present analysis and numerical 
techniques will be extended to Include the flow 
of homogeneous, non homogeneous, and chemically 
reacting species In one- and multi-dimensional 
systems. 


Basic Theoretical Formulation 


The physical model considered for the 
present study Is the transient energy transfer by 
radiation in absorblng-eml ttlng gases bounded by 
two parallel gray plates (Fig. 1). In general, 
Tj and T 2 can be a function of time and position 
and there may exist an Initial temperature 
distribution In the gas. It Is assumed that the 
radiative energy transfer In the axial direction 
Is negligible In comparison to that In the normal 
direction. 

For radiation participating medium, the 
equations expressing conservation of mass and 


momentum remain unaltered, while the conservation 
of energy. In general. Is expressed as 1 

p c p BT * d1v (k 9 rad T) ♦ PT ^ ^ - dlv q R 

( 1 ) 

where p Is dynamic viscosity, p Is the 
coefficient of thermal expansion of the fluid and 
♦ Is the Rayleigh dissipation function. For a 
semi-infinite medium capable of transferring 
energy only by radiation and conduction, Eq. (1) 
reduces to 



where q Is the sum of the conductive heat flux 
q c - - k (5T/&y) and the radiative flux q R . 
For the physical model where radiation Is the 
only mode of energy transfer, the energy equation 
can be written as 


P c 



‘w 
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Use of this simplified equation Is made to 
Investigate the transient behavior of a radiation 
participating medium. 

For many engineering and astrophyslcal 
applications, the radiative transfer equations 
are formulated for one-dimensional planar 
systems. For diffuse nonreflecting boundaries 
and In the absence of scattering, the expression 
for the total radiative flux Is given, for a n- 
band gaseous system, by 1 *®’ 23 


q R (y) ■ ej - e 2 


^ 1*1 L, V 1 ",'*' S “O'- 7 

L 
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where 
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Equation (4) Is In proper form for obtaining the 
nongray solutions of molecular species. In fact, 
this Is an Ideal equation for the line-by-line 
and narrow-band model formulations. However, In 
order to be able to use the wide band models and 
correlations, Eq. (4) Is transformed In terms of 
the correlation quantities as 1 * 7 " 11 * 23 
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q R U) = e l ' e 2 



where 


f, - u/u o - y/L; E' • u'/u o ■ z/L; A - A/A q ; 

u ■ (S/A o ) Py; u q - (S/A q ) PL; PS - ^ d« 

It should be noted that Fj w and in Eq. 

(5) represent the values at tW center <Jf the fth 
band and l'(u) denotes the derivative of A(u) 
with respect to u. Upon performing the Integra - 
tlon by parts, Eq. (5) can be expressed In an 
alternate form as 2 ^ 


q R (E) - - e 2 



1 , 

* f >/<*•] A| [| u o1 (E ,- E)]dE‘} (6) 


A direct differentiation of Eq. (6) provides the 
expression for the divergence of radiative flux 
as 

d<l R U) 3 n 5 

ST~ ’ 7 ' A o1 u o1 { l tde («’)/«•] x 

1*1 0 1 

. 3 1 
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Equations (5) through (7) are the most convenient 
equations to use when employing the band-model 
correlations In radiative transfer analyses. 

Upon defining nondlmenslonal radiative heat 
flux by 


0 <*.t) - q^.O/Ce^t) - e 2 (t)] (8) 


Eq. (5) can be written as 

3 n 1 

Q (E.t) - 1 S u q1 {/ c if (E'.t) x 
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where 


C 1 (? * t) ■ F U 1 U.tJ/Ce^t) - e 2 (t)]A o1 

Equation (9) provides the general expression for 
the radiative flux In the nondlmenslonal form. 

Radiative Interaction 

Specific consideration Is given first to the 
case where radiation Is the sole mode of enerqy 
transfer. * 

By defining « (E.t) - T(E,t)/T o with T 0 
representing some constant reference temperature, 
Eqs. (3) and (7) can be combined to yield the 
energy equation In nondlmenslonal form as 

- ft* (E.t)/8t -2 2 (J 5 ♦ (f.t) x 

£ 1-1 o u1 
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where 


V5.t> * < PS |(T>C» e u< (E.t)/aE]/(p C p T o /t m ) > 

The time t In Eq. (10) Is defined as t* - t/t„, 
with t„ representing some characteristic time 
scale of the physical problem; however, for the 
sake of convenience, the asterisk Is left out 
here as well as In further developments. From the 
definitions of ♦(E.t) and <k ul (E.t), It should 
be noted that Eq. (10) Is a nonlinear equation In 
T(E,t). Equation (10), therefore, represents a 
general case of the transient energy transfer by 
radiation between two semi-infinite parallel 
plates. 

As a special case. It Is assumed that the 
entire system Initially Is at the fixed 
(reference) temperature T 0 . For all time, the 
temperature of the upper plate Is maintained at 
the constant temperature equal to the reference 
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temperature, l.e., T 2 * T 0 . The temperature of 
the lower plate Is suddenly decreased to a lower 
but constant value, l.e., J x < T 2 . The problem, 
therefore. Is to Investigate the transient 
cooling rate of the gas for a step change In 
temperature of the lower plate. 

In many radiative transfer analyses. It Is 
often convenient (although not essential) to 
employ the relation for the linearized radiation 
as 


e Uj m - e Uj <V * < d e u/ dT) T w |H J (11) 

where again the subscript 1 refers to the 1th 
band such that Is the wave number location 
of the band and T w represents the temperature of 
the reference wall which could be either T l or 
T z . For the special case considered, since we 
are Interested In Investigating the transient 
behavior of the gas because of a step change In 
temperature of the lower plate, T* Is taken to be 
equal to Tj. It should be pointed out that for a 
single-band gas, the linearization Is not 
required because the temperature distribution can 
be obtained from Eq. (10) and the radiative heat 
flu* can be calculated from Eqs. (5), (6), or 
(9). However, for the case of multiband gases 
and for systems Involving mixtures of gases. It 
Is convenient to employ the linearization 
procedure In order to use the Information on band 
model correlations. The following definitions 
are useful In expressing the governing equations 
In linearized forms: 


0 ■ (T-T 1 )/a 2 -T l ) (12a) 

N 11 • (P V P c p> K 11 * K 11 * S 1< T > < d 

(12b) 

n 
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M 11 ' ( V L p V h h * H » ' A oi (T)(d e <,i /dT) r 1 

(12d) 

M 1 * (t m /L p V H 1 • H 1 ’ * H n ( 12e) 


M 11 u o1 ’ N 11 • u ol H 11 ’ PL K 11 < 12f > 

where Hj, Kj, Nj end Nj represent the values of 
H, K, N end H evaluated at the temperature T 1# 
As explained In Refs. 1 and 8, these quantities 
represent the properties of the gaseous medium. 


By employing the definitions of Eqs. (11) 
and (12), Eq. (10) Is transformed to obtain a 
convenient form of the energy equation as 2 ** 


aeU.t) 
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(13) 


The parameters In Eq. (13) are Nj and u„. For a 

given gas, the parameters are the gas pressure 
and the temperature of the lower wall. 

The Initial and boundary conditions for the 
physical problem considered are 

e(E.O) » 1 ; e(O.t) - 0 ; 9(1, t) - 1 ( 14 ) 

It Is important to note that the boundary 
conditions given In Eq. (14) are not applicable 
to Eq. (13) because this equation does not 
require a boundary condition. Thus. In this 
case, the temperature of the medium adjacent to a 
surface differs from the surface temperature. 
This Is because the temperature of the medium 
adjacent to a surface Is affected not only by the 
surface but also by all other volume elements and 
surfaces. The radiation slip method Is a means 
of accounting for such temperature jumps and this 
Is discussed In Ref. 1. 

By employing the definitions of Eqs. (8), 
(11), and (12), relations for the radiative flux, 
as given by Eqs. (5) and (6), are expressed as 

q(E.t) ■ 1- (3/8 o T?) I u | H.. <J 5 eu\t) x 
1 1*1 01 11 o 
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A i C 7 u o1 ( ^‘^«‘ * / Cl - 8{f,t)] x 

A i c ! u oi<* , -«J« , » (15a) 


and 
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1 1*1 11 o 85 
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It should be pointed out that Eq. (15a) Is a 
convenient form for the optically thin and 
general solutions while Eq. (15b) Is useful for 
solutions In the large path length limit. Once 
the solutions for 0(5* t) are known from the 
energy equation, the appropriate relations for 
the heat flux can be obtained from Eqs. (15). It 

should be noted that the quantity H./(aT,) In 
Eqs. (15) Is nondlmenslonal . 1 1 

Radiation and Conduction 


For this case where conduction heat transfer 
takes place simultaneously with radiative 
transfer, the energy equation is given by Eq. 
(2). Thus, a combination of Eqs. (2), (11) and 
(12) results In 


as » r a 2 e . 3 " , f 5 90 ( 5 ', t) 


*1*1 u o1 <5-5’)]d5’ ' J — x 
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(16) 


where 


R * k t /(p C (/) 

m p 


The nondlmenslonal parameter R In Eq. (16) Is 
analogous to the Fourier number. For R * 0, Eq. 
(16) reduces to Eq. (13). Since the presence of 
conduction Implies continuity of temperatures at 
the boundaries, the boundary conditions for Eq. 
(16) are those given In Eq. (14). The quantity 
N|/R can be expressed as N ■ Nj/R * (PlvkJKj. 
The nondlmenslonal parameter N denotes the rela- 
tive Importance of radiation versus conduction In 
the gas. For particular values of P and L, It Is 
actually the dimensional gas property Nj/(RPL^) * 
Kj/k that represents the relative Importance of 
radiation versus conduction. 

For the case of no radiation, Eq. (16) 
becomes 


be 



(17) 


The separation of variables results in a general 
solution of Eq. (17) as 


0 - exp(-X 2 Rt) (Bj sin X5 + B 2 cos X 5) (18) 

2 

where X Is the separation parameter. The 
particular solution of Eq. (18) can be obtained 
by satisfying the boundary conditions. 
Alternately, by defining a similarity variable 

n * 5/^Rt, Eq. (17) can be written as 

5-J+2n^* 0 (19) 


The solution of Eq. (19), with Initial and 
boundary conditions given by Eq. (14), Is found 
to be 


0(5*t) » (2//5) erf n (20) 


This solution Is applicable for relatively large 
separations between the plates. 

In the case of simultaneous conduction and 
radiation heat transfer, the nondlmenslonal heat 
flux Is defined as 


0 - tq c (5*t) + qaU.tn/tejU) - e 2 (t) ] (21) 


Expressions similar to Eqs. (15a) and (15b) can 
be obtained easily In this case also (Ref. 24). 

Method of Solutions 


The solution procedures for the radiative 
equilibrium and radiation with conduction cases 
are available In Ref. 24. For the sake of 
brevity, only the solution procedure for the 
radiative equilibrium case Is given here. 

For the general case of radiative 
equilibrium, the temperature distribution Is 
obtained from the solution of the energy 

equation, Eq. (13). Once 0(5. t) Is known, the 
radiative heat flux Is calculated by using the 
appropriate form of Eq. (15). Before discussing 
the solution procedure for the general case, 
however. It Is desirable to obtain the limiting 
forms of Eqs. (13) and (15) In the optically thin 
and large path length limits and Investigate the 
solutions of resulting equations. 

Optically Thin Limit 

In the optically thin limit X(u) * u and 

A'(u) - 1, and therefore, Eq. (13) reduces 
^ 1 , 8,23 


♦ 3 Nj 0(E,t) - | Kj - 0 (22a) 


From an examination of Eq. (22a) along with the 
definitions given In Eq. (12), it is evtdient 
that in the optically thin limit the temperature 
distribution In the medium is Independent of the 
E-coordinate for the case of pure radiative 
exchange. This Is a charateristic of the 

optically thin radiation in the absence of other 
modes of energy transfer. Thus, Eq. (22a) can be 
written as 


ONj 8(t) - | Mj - 0; 0(5,0) - 1 

(22b) 

Since gas properties are evaluated at known 
reference conditions, Nj Is essentially constant, 
and solution of Eq. (22b) Is found to be 
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0(t) - j£l + exp(- 3 Njt)] (23) 

In the optically thin limit both forms of 
Eq. (15) yield the same final relation for the 
radiative flux as 2 ® 


Q(5.t) * 1 - [3/(8 <jTj)] (PLKj) [(1-0 + 

♦ (25 - 1) 0(t)] (24) 


of variation of parameters. For this, a 
polynomial form for 8(5,t) Is assumed In powers 
of 5 with time dependent coefficients as 

0(5.t) - Z c (t) 5* (27) 

m«o 

By considering only the quadratic solution In 

c 5, ..V,* 1 sa *1sfy1ng the boundary conditions of 
Eq. (14), one finds 


It should be pointed out that In Eq. (24) the 
quantity (PLKj/oT®) is nondlmenslonal . The 

relation for 0( t) In Eq. (24) Is obtained from 
Eq. (23). Thus, evaluation of the temperature 
distribution and radiative heat flux In the 
optically thin limit does not require numerical 
solutions. 


Large Path Length Limit 

In the large path length limit (l.e., for 
u 0 j » 1 for each band), one has l(u) - 

• *n(u), J'(u) » 1/u, and J"(u) » - 1/u 2 . 
Thus, In this limit, Eq. (13) reduces 


80(5, t) 

at 
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(25) 


An analytical solution of Eq. (25) may be 
possible, but numerical solution can be obtained 
quite easily. 


In the large path length limit, Eqs. (15a) 
and (15b) reduce respectively to 


0(5. t) 


1 - (l/a4Tj) 
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‘ ^ TTTT 1 (26a) 
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The expressions for dimensionless radiative heat 
flux from or to the wall are obtained by setting 
5-0 In Eqs. (26). 


0(5. t) - 5 2 + g(t) (5-5 2 ) 


(28) 


where g(t) represents the time dependent 
coefficient. At t » 0, a combination of Eqs. 
(14) and (28) yields the result 


9(0) - (l-5 2 )/(5-5 2 ) (29) 


Equations (28) and (29) are used to obtain 
specific solutions of Eqs. (13) and (25). Once 
the temperature distribution Is known, the 
expressions for heat flux are obtained by 
numerically Integrating the corresponding 
equations. The entire numerical procedure Is 
described In detail In Refs. 23 and 24. The 

numerical procedure is similar for higher order 
solutions In 5 of Eq. (27), but computational 
resources required are considerably higher. 

Physical Conditions and Oata Source 

For the physical problem considered (Fig. l) 
four specific absorblng-eml ttlng species were 
selected for an extensive study; these are CO, 
C0 2 , OH and H 2 0. The species CO was selected 
because It contains only one fundamental 
vibration-rotation (VR) band and all spectral 
Information are easily available In the 

literature. It Is a very convenient gas to test 
the numerical procedure without requiring 
excessive computational resources. Species OH 
and are the primary radiation participating 
species for the pressure and temperature range 
anticipated In the combustor of the scramjet 
engine. Species C0 2 , and combinations of C0 2 and 
H 2 O are Important absorbing-emitting species In 
many other combustion processes. 

In radiative transfer analyses. It Is 
essential to employ a suitable mode 1 to represent 
the absorption-emission characteristics of 
specific species under Investigation. Several 
line-by-line, narrow-band, and wlde-band models 
are available to model the absorption of a VR 
band (Refs. 7-11). However, It Is often 
desirable to use a simple correlation to 
represent the total absorption of a wide band. 
Several such correlations are available In the 
literature (Refs. 7-11). The relative merits of 
these correlations are discussed In Ref. 11. In 
this sludy, the correlation proposed by Tien and 
Lowder' Is employed and this Is given by 


Numerical Solutions of Governing Equations 

The general solutions of Eqs. (13) and (25) 
are obtained numerically by employing the method 


A(u) * In (uf O) ty ^ + (30) 
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where 


f(P) = 2. 94 [ 1 - exp(- 2*60 p)] 


and p represents the line structure parameter. 

The spectral Information and correlation 
quantities needed for CO, C0 2 and H 2 0 were 
obtained from Refs. 7-9. The spectral data for 
OH were obtained from Ref. 25 and correlation 
quantities were developed In Ref. 24. The 
specific VR bands considered for each species 
are: CO ( 4. 7u fundamental), OH (2.8p 

fundamental), CO2 (15 p, 4.3n and 2.7p), and H2O 
(20|i, 6.3»i, 2.7n, 1.87}*, and 1.38^). 

For the specific problem considered, the 
dependent variables are 0, Q, and 3. and 

Independent variables are * and t. The 
parameters for general solutions are Ti, Tg/Tj, 
u Q , and t^. For the radiative equilibrium case, 
9 and Q depend only on t and Nt In the 
optically thin limit and on t, and In the 
large path length limit. For the case of 

combined radiation and conduction, 0 and 0 depend 
on 5, t, and Mj/R In the large path length limit. 
The parameters for specific solutions for differ- 
ent species are Ti, T 2 /T!, p< and V* 
Extensive results, therefore, can be obtained by 
varying these parameters. For parametric 
studies, however, only certain values of 
pressure, temperature, and plate spacing were 
selected and results were obtained for the 
general as well as limiting cases. Unless stated 
otherwise, specific results were obtained for To 

■ 2 v 

Results and Discussion 

Results have been obtained for different 
radiation participating species for both cases, 
the radiative equilibrium and radiation with 
conduction. Selected results are presented here 
and extensive results are provided In Ref. 24. 
It should be realized at the outset that, 
according to the physics of the problem, the gas 
Initially Is at a high temperature Ti * T 2 . At t 
* 0, the temperature Ti Is lowered to a constant 
value. The energy excnange then occurs and the 
gas cools down In time until a steady-state 
condlton Is reached. At this time, a certain 
temperature profile Is established and a fixed 
amount of energy exchange occurs Irrespective of 
the time. The rate of cooling of the gas layer, 
therefore. Is dependent on the nature of the 
participating species and on the physical 
parameters of the problem. 

Some limiting solutions that are Independent 
of any participating species are presented first 
In Fig. 2 for the radiative equilibrium case. 
The temperature distribution In the channel Is 
plotted as a function of the optically thin 
parameter Nj for different times. These results 
show an exponential decay with time reaching the 
steady-state value of 0*1/2 for t+®. The 
temperature distribution for the large path 
length limit Is shown In Fig. 2 (with broken 
lines) as a function of the large path length 
parameter Mj and for different times. Although 


the numerical values are entirely different, 
these results also show the exponential decay 
with time and reach the limiting value of 
0 * 1/2 for t + ®. It should be noted that 
while the optically thin solutions are 

Independent of the ^-coordinate, the large path 
length solutions do depend on l and they have 
been obtained for E * 0.5. In the case of 
simultaneous radiation and conduction, both 

optically thin and large path length solutions 
for temperature distribution depend on £. These 
results, however, can be expressed In terms of 
the radiation-conduction parameter N*Ni/R In the 
optically thin limit and M - Mj/R In the large 
path length limit. 

The radiative equilibrium temperature 
distribution for CO are shown In Fig. 3 for three 
different characteristic times and for P * 1 atm, 
L * 10 cm, and T w ■ 500 K. For small t^, t* ■ 
t/t„, becomes large and, therefore, 0 varies 
slowly with t*; the reverse Is true for large 
values. This trend Is evident clearly from the 
simplest case of the optically thin solution 
given by Eq. (23). Similar trends In results 

were observed also for other species for 
different values of P, T w and L. Thus, to 

demonstrate typical transient trends, other 
results presented In this study were obtained for 
an Intermediate value of the characteristic time 
tn, * 0.00001 sec. 

The centerline temperature variations with 
time are Illustrated In Fig. 4a for CO and In 

Fig. 4b for OH. General and limiting solutions 
are shown for pure conduction, radiation, and 
radiation with conduction for P * 1 atm, T w * 500 
K, and L * 10 cm. It Is noted that for both 

gases the optically thin solutions approach the 
steady-state conditions faster than the large 
path length and general solutions. For the 
physical conditions considered, the energy Is 
transferred faster by conduction than by 
radiation, and the steady-state conditions are 
reached earlier by the combined radiation and 
conduction process. Although OH Is a relatively 
better heat conducting gas, CO Is seen to be more 
effective In the radiative transfer. For the 
same physical conditions as In Figs. 4, the 
radiation and radiation with conduction results 
are compared for CO, OH, H2O, and C0 2 In Fig. 5. 
It Is seen that H2O Is most effective and OH Is 
least effective In transferring the radiative 
energy. The ability of a gas to transfer 
radiative energy depends on the molecular 
structure of the gas, band Intensities and 
physical conditions of the problem. Thus, H2O 
with five strong VR bands Is a highly radiation 
participating species and the steady-state 
conditions are reached quickly for H2O than for 
other species. However, CO with one fundamental 
band Is seen to be a better radiating gas than 
CO2 with three VR bands. This Is because for the 
given physical conditions, the optical thickness 
of C0 2 Is sufficiently large and In the large 
path length limit C0 2 Is relatively less 
effective In transferring the radiative energy 
(Ref. 8). Further results for CO and OH are 
Illustrated In Figs. 6 for different wall 
temperatures. It Is seen that while radiation Is 
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less effective than conduction at T * 500 K, It 
Is highly effective at T w * 1,000 K. This, 
however, would be expected because radiation 
becomes considerably Important at higher tempera- 
tures. The steady-state condition Is reached 
quicker for T w * 1,000 K than for T w * 500 K. In 
fact for the characteristic time considered (t- - 
0.00001 sec.), the steady-state condition Is 
reached quickly for all species for temperatures 
higher than T w * 1,000 K. Results for the pure 
radiation case are Illustrated In Figs. 7 for CO 
and OH for L * 1 cm and 10 cm. It Is seen that 
while the general and large path length solutions 
depend on the plate spacing the optically thin 
solutions are Independent of the spacing. This 
fact was pointed out earlier In the method of 
solution. In the presence of other modes of 
energy transfer, the optically thin solutions 
also depend on the plate spaclngs. As would be 
expected, for the same physical conditions, the 
steady-state condition Is reached quicker for the 
lower plate spacing. 

The temperature variations within the plates 
are shown In Figs. 8-11 for different species and 
for P - 1 atm and L - 5 cm. In the absence of 
molecular conduction, temperature jumps (radia- 
tion slips) occur at the boundaries and, there- 
fore, the general solutions for the case of 
radiative equilibrium are not presented. 
However, general as well as limiting solutions 
are presented for the case of radiation with 


conduction. 

It Is noted 

, In 

general , 

that 

for 

the case of 

radiation 

with 

conduction. 

the 

steady-state 

conditions 

are 

reached 

for 

all 

species at t 

>0.1 and 

for 

1 > 500 

K. 

For 

the case of 

pure radiation, 

the steady-state 


conditions are reached at relatively longer 
tiroes. The optically thin results are seen to be 
independent of the {-coordinate for the case of 
radiative equilibrium and are seen to vary slowly 
In the central portion of the plates for the case 
of radiation with conduction. This Is because, 
In this limit, the gas Interacts directly with 
the boundaries and conduction Is predominant near 
the walls. Specific results for CO are 

illustrated in Figs. 8 for T^ * 1,000 K. For the 
case of radiation with conduction, general and 
limiting solutions are compared In Fig. 8a; and 
for both cases, the radiative equilibrium and 
radiation with conduction, limiting solutions are 
compared In Fig. 8b. The steady-state results 
for pure conduction are also shown In Fig. 8b for 
comparative purposes. The results demonstrate 
the typical trends for limiting and general 
solutions, l.e., a lower temperature gradient 
implies a higher rate of energy transfer. 
Specific results for OH are Illustrated In Figs. 
9 for T w * 500 K. General and limiting solutions 
are shown In Fig. 9a; and limiting solutions are 
compared in Figs. 9b and 9c. For the case of 

radiation with conduction, the limiting and 
general solutions are seen to be In good 
agreement for all times (Fig. 9a). This Is 

because for the conditions of the Illustrated 
results, conduction dominates the energy transfer 
process in OH. The typical trends In results for 
the optically thin and large path length limits 
are shown In Figs. 9b and 9c, respectively. 
Figure 9b clearly shows that for all times the 
radiative equilibrium results are independent of 
the {-coordinate In the optically thin limit. 


Figure 9 shows that at earlier times the rate of 
energy transfer Is higher In the presence of 
conduction. Specific results for H 2 0 are 
Illustrated In Figs. 10 for the case of radiation 
with conduction. It is seen clearly that the 
rate of cooling Is slgnficantly higher In the 
large path length limit (Fig. 10a), and the 
steady-state conditions are reached at relatively 
longer times for lower T* values (Fig. 10b). For 
the case of combined radiation and conduction, a 
comparison of results for different species Is 
shown In Fig. 11 for t - 0.01 and 0.1. The 
results for t * 0.1 essentially correspond to the 
steady-state conditions. For t - 0.01, the 
variation in temperature is seen to be relatively 
small between {-0.2 and 0.9. The centerline 
temperature Is found to be the lowest for H 2 0, 
and this Is followed by OH, CO, and C0 2 . 
However, It is noted that OH is very effective In 
transferring the net energy In comparison to the 
other species. As discussed earlier, this Is 
mainly due to relatively higher conductive 
ability of OH at T w - 500 K. 

The centerline temperature distributions are 
shown In Figs. 12-15 for different gases as a 
function of the spacing between the plates. In 
most figures, results are presented for both 
cases, the radiative equilibrium and radiation 
with conduction. In selected figures, results 
for the case of pure conduction are Included also 
for comparative purposes. For a particular gas, 
specific results are presented for various times 
to demonstrate the radiative nature of the gas 
under different pressure and temperature 
conditions. 

The results for CO are presented In Figs. 12 
for different cases. For P • 1 and T w - 500 K, 
the results Illustrated in Fig. 12a show that the 
time required to reach the steady-state condition 
Increases with increasing plate spaclngs. For a 
particular plate spacing, the centerline tempera- 
ture Is lower for the case of radiation with 
conduction than for pure radiation for all times. 
For P - 1 atm and T w - 1,000 K, results presented 
In Fig. 12b show that the large path length 
solutions are closer to the general solutions for 
L > 20 cm; and the results for pure radiation and 
radiation with conduction are Identical for 
t > 0.5. The centerline temperature variations 
are shown In Fig. 12c for t - 0.5 T - 500 K, and 
different pressures. It Is noted that while the 
heat transfer by conduction Is Insensitive to the 
change in pressure, the radiative heat transfer 
Is strongly dependent on It. The rate of 

radiative interaction increases with increasing 
pressure until the large path length limit Is 
reached for sufficiently large values of L. For 
the case of pure radiation, the results for P - 
0.1 atm differ considerably from other results. 
This Is due to use of the Tien and Lowder's 
correlation which is suitable only at relatively 
higher pressures (Ref. 11). The centerline 
temperature variations are shown In Fig. 12d for 
t « 0.5, P - 1 atm, and different values of T w . 
As would be expected, both conductive and radia- 
tive interactions Increase with Increasing 
temperatures, although the Increase In radiative 
transfer is comparatively higher. It should be 
noted that for T w - 300 K, T 2 - 2 T w - 600 K, for 
T w - 500 K, T 2 « 1,000 K, and so on. Thus, for a 
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higher value of T = Tj, the energy Interactions 
occur at a sufficiently large temperature differ- 
ence between the upper and lower plates. At 
these temperatures. If the plate spacing Is 
small, the energy Is transferred quickly and the 
steady-state condition Is reached at relatively 
shorter times. This fact was pointed out also In 
the discussion of results of Figs. 7. 

The centerline temperature variations for OH 
are Illustrated In Figs. 13 for different condi- 
tions. The results presented In Fig. 13a for P * 
1 atm and T w * 500 K show the similar trend as CO 
In Fig. 12a, although the extent of energy trans- 
fer by simultaneous radiation and conduction Is 
relatively higher. This Is because at T w * 500 
the energy transfer In OH Is dominated by the 
conduction heat transfer. General and limiting 
solutions for radiative equilibrium are shown In 
Fig. 13b and for radiation with conduction In 
Fig. 13c. These results clearly demonstrate the 
typical radiative Interaction trends for 
different times. The results show that the 
optically thin solutions are Independent of the 
plate spacing In the case of pure radiation but 
depend on the spacing when molecular conduction 
Is Included. The large path length results are 
seen to be valid only for large values of L for 
the case of pure radiation (Fig. 13b), but they 
appear to be valid In the entire range for the 
case of radiation with conduction (Fig. 13c). The 
results for pure conduction, pure radiation, and 
conduction with radiation are Illustrated In Fig. 
13d for P * 1 atm and T w » 1,000 K. For this 
temperature, the results for pure radiation and 
radiation wlh conduction are found to be 
identical. This Indicates that at higher temper- 
atures, OH becomes a highly radiation participat- 
ing gas. The results for variation of the 
centerline temperature for OH with pressure and 
temperature are given In Ref. 24 and they show 
the same general trend as the results for CO 
shown in Figs. 12c and 12d. 

Extensive results of e(£ ■ 0.5) versus L 
have been obtained for H^O and CO 2 for different 
conditions and these are discussed In Ref. 24. 
In general, these results show similar trends as 
exhibited by the results for CO and OH but the 
extent of radiative Interactions is entirely 
different. For example, a comparison of results 
presented In Fig. 14a for H^O with the results of 
Fig. 12a for CO and Fig. 13a for OH for Identical 
conditions reveals that HpO Is a highly radiation 
participating gas even for shorter times. For 
HpO, results presented In Fig. 14b demonstrate 
that the large path length solutions are closer 
to the general solutions for both cases, the 
radiative equilibrium and radiation with 
conduction. 

The centerline temperature variations are 
compared for different gases In Fig. 15 for P * 1 
atm, T w « 500 K, and t - 0.05. For the case of 
radiative equilibrium. It Is noted that OH is the 
least effective and H 2 0 Is the most effective gas 
In transferring the radiative energy for plate 
spaclngs greater than two centimeters. When 

molecular conduction Is Included, OH becomes more 
effective because of Its relatively higher 
conductive ability. These points were noted also 
In earlier discussions. The story, however, can 
be entirely different for other physical 


conditions because of the radiative/conductive 
nature of participating species (Refs. 8 and 
24). This fact Is partially evident from the 
steady-state results, for the case of combined 
radiation and conduction, presented In Fig. 16 
for two different temperatures, T w ■ 300 K and 
500 K. For example, for T w * 30o K and L * 10 

cm, the temperature values for CO and C0 2 are 
about the same, for H 2 0 It Is lower, and for OH 
it Is the lowest; however, for plate spacing 
greater than L * 20 cm, the trend Is entirely 
different. Also, It should be noted that the 
steady-state (t * 0.5) results for T w * 500 K In 
Fig. 16 show different trend than the results for 
the same temperature In Fig. 15 for t * 0.05. 
Thus, In order to predict the relative ability of 
a gas for radiative interactions. It Is very 
Important to consider the exact physical 
conditions for all species. These predictions 
may not be applicable If physical conditions of 
the problem are changed. 

Extensive results for the variation In heat 
transfer can be presented analogous to the 
variation of temperature for different condi- 
tions. However, this should not be necessary 
because the heat transfer variation follows the 
trend of the temperature variation In a reverse 
manner. If the temperature differences are 
higher, the rate of heat transfer will be higher 
and the steady-state conditions will be reached 
at earlier times. The results for heat transfer 
variations have been obtained for selected 

conditions and these are available In Ref. 24. 
Only a few results are presented here to show the 
trend In cooling rates for different wall 
temperatures and plate spaclngs. 

For P » 1 atm, the results^for Q and Q are 
Illustrated, as a function of t , In Figs. 17 for 
H 2 0 and In Fig. 18 for COp. The results for t * 
0.0 - 1.0 are shown in Fig. 17a for HpO. 

However, It Is found that for the t^ value 

selected In this study, the steady-statg 
conditions are reached In most cases at about t 

* 0.2 Consequently, the results In Figs. 17b^ 
17c, and 18 are presented only In the range of t 

* 0.0 - 0.2 to demonstrate the rate of cooling at 
different temperatures. As would be expected, 
the results show that for a given plate spacing 
the gas layer reaches the steady-state condition 
faster at higher values of T w because of stronger 
radiative Interactions. It should be noted that 
the rate of energy transfer increases with time 
for a gas layer closer to the upper wall 

(£ * 0.75) and decreases with time for a gas 
layer closer to the lower wall U * 0.25) until 
the steady-state conditions are reached. The 
rate of cooling 1$ entirely different If the 
plate spacing Is changed (Fig. 17c). From a 
comparison of results of Figs. 17 and 18, It Is 
noted that the trend and rate of energy transfer 
are different for different species. This, 
however, would be expected because of the rela- 
tive ability of different species to participate 
in the radiation-conduction Interaction process. 

Conclusions 

The problem of transient radiative Interac- 
tion In nongray absorbing-emitting species has 
been formulated In a general sense such that 
sophisticated absorption models can be used to 
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obtain accurate results if desired. Results have 
been obtained for the special case of radiative 
Interactions In a plane gas layer bounded by two 
parallel plates when the temperature of the 
bottom plate Is suddenly reduced to a lower but 
constant temperature. The energy transfer by 
pure radiation, and by simultaneous radiation and 
conduction were considered and specific results 
have been obtained for CO, C0 2 , H 2 0, and OH by 
employing the Tien an Lowder's correlation for 
band absorption. It Is noted that the extent of 
radiative Interaction Is dependent on the nature 
of the participating species and parameters Tj, 
T 2 /Tj, P, L, and t^. The steady-state conditions 
are reached at relatively longer times for radia- 
tive equilibrium than for radiation with conduc- 
tion. For a particular value of P and T|, the 
time required to reach the steady-state condition 
Increases with Increasing plate spacing. For a 
fixed plate spacing, the energy Is transferred 
quickly for higher Tj values because of large 
temperature differences between the plates. The 
rate of radiative Interaction Increases with 
Increasing pressure until the large path length 
limit is reached. The radiative equilibrium 
solutions are found to be Independent of the 
plate spacing In the optically thin limit. In 
the case of simultaneous radiation and conduc- 
tion, both optically thin and large path length 
solutions depend on the y-locatlon between the 
plates. At moderate temperatures, OH Is a poor 
radiating but better heat conducting gas. For 
most conditions, H 2 0 Is found to be highly 
radiation participating species, and the steady- 
state conditions are reached quickly for H?0 than 
for other species. The Information on radiative 
Interactions for OH and H?0 Is useful In the 
analysis of the flow field In the scramjet 
engine. 
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Optically thin 



Fig. 2 Limiting solutions for the radiative 
equil Ibrlum case. 



Fig. 3 Temperature variation for transient 
radiation for different characteristic times 
(CO, P=1 atm, T =50C K,L=10cm, and £=0.5). 


General Gas : CO 


Large path length E = 0. 5 



Fig. 4a Centerline temperature variation with 
time for CO with P*1 atm, T *500 K, and 
L=10 cm. * 
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Fig. 1 Physical model and coordinate system. 


General Gas ; OH 



Hg. 4b Centerline temperature variation with 
time for OH with P* 1 atm. T • =500 K, and L*1C 
cm. " 
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Fig. 6a Centerline temperature variation with 
time for CO with P=1 atm and L =1 0 cm. 
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Fig. 7a Centerline temperature results for 
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Abstract 

Analysis and numerical procedures are 
presented to Investigate the transient radiative 
Interactions of nongray absorblng-eml ttlng 
species In laminar fully-developed flows between 
two parallel plates. The particular species 
considered are OH, CO, CO 2 * and HoO and different 
mixtures of these species. Transient and steady- 
state results are obtained for the temperature 
distribution and bulk temperature for different 
plate spaclngs, wall temperatures, and pressures. 
Results, In general, Indicate that the rate of 
radiative heating can be quite high during 
earlier times. This information Is useful In 
designing thermal protection systems for 
transient operations. 

Nomenclature 

A band absorptance * A(u»p), cm 1 

A Q band width parameter, cm -1 

C Q correlation parameter, atm -1 - cm* 1 

C D specific heat at constant pressure, 

kJ/kg-k * erg/gm-k 

e Planck's function, (W~cnr*)/cm -1 

a) 

e Planck's function evaluated at wave 

w o number w 

0 

ej,e 2 emissive power of surfaces with 

temperatures Tj and T 2 , W-cnf* 

Hji,Ht gas property for the large path length 
limit 

k thermal conductivity, erg/cm-sec-K 

Kj gas property for the optically thin limit 

distance between plates 

M large path length parameter, 

nondlmenslonal 

N optically thin parameter, nondlmenslonal 

P pressure, atm 

q D total radiative heat flux, w/cm* 
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q n spectral radiation heat flux, 

q wall heat flux, w/cm* 

^ 0 ) 

S Integrated Intensity of a wide band, 

atm* 1 -cm - * 

t time, sec 

T temperature, K 

Tj wall temperature, K; Tj * T w 

T b bulk temperature, K 

u nondlmenslonal coordinate * $Py/A 0 

u Q nondlmenslonal path length * SPL/A Q 

y transverse coordinate, cm 

p line structures parameter 

0 nondlmenslonal temperature 

0^ dimensionless bulk temperature 

k spectral absorptln coefficient, cm -1 

hi 

1 nondlmenslonal coordinate - y/L - u/u Q 

p density, kg/m* 

t nondlmenslonal time 

(i> wave number, cm -1 

u> wave number at the band cener, cm -1 

o 

Introduction 

The radiative energy transfer In participat- 
ing medium has received special attention in 
recent years because of Its applications in the 
areas of the remote sensing, earth's radiation 
budget studies and climate modeling, fire and 
combustion research, entry and reentry phenomena, 
hypersonic propulsion and defense-oriented 
research. In most studies Involving combined 
mass, momentum, and energy transfer, however, the 
radiative transfer formulation has been coupled 
mainly with the steady processes [1-13] and the 
interaction of radiation In transient processes 
has received very little attention. However, the 
transient approach appears to be the logical way 
of formulating a problem In general sense for 
elegant numerical and computational solutions. 
The steady-state solutions can be obtained as 
limiting solutions for large times. 

The limited number of studies available on 
the transient radiative transfer In gaseous 


1 


systems [14-24] are reviewed critically In [25- 
27], The literature survey reveals that the 
transient behavior of a physical system can be 
Influenced significantly In the presence of 
radiation. 


and u Is given by the well-known parabolic 
profile as 


U * 6 (5 - c 2) ; 5 ■ y/L (3) 


The goal of this research Is to Include the 
nongray radiative formulation In the general 
governing equations and provide the step-by-step 
analysis and solution procedure for several 
realistic problems. The basic formulations are 
presented In [25] and the specific case of 
transient radiative exchange In nongray gases 
between two parallel plates Is investigated In 
[26, 27]. The objective of this study Is to 
Investigate the Interaction of transient radia- 
tion In fully-developed laminar flows between two 
parallel plates. Thus, special attention Is 
directed to Include the nonsteady phenomenon only 
In the energy equation. Homogeneous as well as 
nonhomogeneous mixtures of absorblng-eml ttlng 
species are considered. In subsequent studies, 
the present analysis and numerical techniques 
will be used to investigate the unsteady flow of 
compressible and chemically reacting species In 
one- and mul tl -dimensional systems. 


where ^ represents the mean fluid velocity. 
Also, for the flow of a perfect gas with uniform 
heat flux, 5T/dx Is constant and Is given by 


6T/6x - (2oq w )/(u m L/k) (4) 


Mow, by combining Eqs. (2) - (4), the energy 
equation Is expressed In nondlmenslonal form as 


06 

ST 


♦ 12 (C - 


< 2 > 


d 2 e 






where 

T » «t/L 2 ; 9 « (T - T^/^L/k) 


(5) 


Basic Formulation 

The physical system considered Is the energy 
transfer In laminar. Incompressible, constant 
properties, fully-developed flow of absorbing- 
emitting gases between parallel plates (Fig, 1). 
The condition of uniform surface heat flux for 
each plate is assumed such that the temperature 
of the plates varies in the axial direction. 
Extensive treatment of this problem Is available 
in the literature under steady state conditions 
[B, 11]. The primary motivation of this study Is 
to Investigate the extent of transient radiative 
Interaction for high temperature flow conditions. 


By assuming that the Initial temperature 
distribution In the gas is some uniform value T 0 
■ Tj, the Initial and boundary conditions for 
this problem can be expressed as 


9(5. 

0) - 0 

(6a) 

9(0. 

i) " 9 (1, %) ■ 0 

(6b) 

V' 

- 1/2) - 0 ; 6^(5 - 0) - - 0^(5 - 1) 

(6c) 


For the physical problem considered, the 
energy equation can be expressed as [1] 



+ n(f£) 2 - d1vq R (1) 


where u and v denote x and y components of 
velocity, respectively. In deriving Eq. (1) It 
has been assumed that the net conduction heat 
transfer In the x direction Is negligible 

compared with the net conduction In the y direc- 
tion. This represents the physical condition of 
a large value of the Peclet number. By an 
analogous reasoning, the radiative heat transfer 
In the x direction can be neglected in comparison 
to that transferred In the y direction. If, In 
addition. It is assumed that the Eckert number of 
the flow Is small, then Eq. (1) reduces to 


dT 

ST 


♦ 




5 2 T 


1 


(2) 


where a - (k/pC ) represents the thermal 
dlffuslvlty of tWrluld. 


For a steady fully-developed flow, v » 0, 


It should be noted that all the boundary 
conditions given In Eqs. (6) are not Independent 
and any two convenient conditions can be used to 
obtain solutions. Also, the Initial temperature 
distribution can be any specified or calculated 
value of 6(t, 0) - f(0. 

As discussed In Refs. 25-27, the radiative 
transfer equations are formulated for one- 
dimensional planar systems for many engineering 
and astrophyslcal applications. For diffuse 
nonreflecting boundaries and in the absence of 
scattering, the expression for the total 
radiative flux 1$ given, for an n-band gaseous 
system, by [1, 8, 25] 


q R (y) - - e 2 


* * A L, Fl “< ,! ’ S eKpt ' ^ 

L , 

- / F^U) exp [- { ^(z-yHdzJd^ (7) 
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where 


F i W| (z) " \ {z) ' e i.y f 2 Wi (zI * \ {z) ' e 2uj 

Equation (7) is in proper form for obtaining the 
nongray solutions of molecular species. In fact, 
this is an Ideal equation for the line-by-line 
and narrow-band model formulations. However, In 
order to be able to use the wide band models and 
correlations, Eq. (7) Is transformed In terms of 
the correlation quantities as [1, 7-13, 25] 


Equations (8) through (10) are the most 
convenient equations to use when employing the 
band-model correlations In radiative transfer 
analyses. 

For the present physical problem, ej * e 2 
and F. * F 0 . Thus, for the case of 

luij 2 Wj 

linearized radiation, a combination of Eqs. (5) 
and (8) results In [8, 18] 

0^ - 9 t - 3M0 - 12 (5*5 2 ) 


**R * e l ' e 2 


*1 * A d u o« { ' F i* w *; [ l u oi 
1=1 0 1 


/ *i [ ! u o1 (8) 


where 


5 « u/u Q = y/L; V « u 7 u q “ z/L; A * A ^ A 0 ; 

u * (S/A ) py; u * (S/A ) PL; PS = / k dw 
0 ° Aw 

It should be noted that and Fg In Eq. (8) 

represent .the values at the center of the. 1th 
band and A'(u) denotes the derivative of A(u) 
with respect to u. Upon performing the integra- 
tion by parts, Eq. (8) can be expressed In an 
alternate form as [25] 


q R (PJ * ej - e ? 


♦ Z A . {/ [de (5')/d5'] A. [I u ol (5-5‘i 
1-1 01 o “1 12 01 




+/ [de^ (C’)/d5’] u of <5’-5)]dt'} (9) 

A direct differentiation of Eq. (9) provides the 
expression for the divergence of radlatve flux as 


dq D U) i n £ 

R - 3 s A„ 4 u „, {/ [de {?•)/«•] X 
0 1 


— * 7 * xl A 0l u o1 


a; [| u Ql (5-5')]d5' - J [de^(5')/d5'] x 
S; [f U o1 (5'-O]d5’) (10) 


- | (L/k) Z H,. u 2 . {/ 5 0(5’, x) x 
1*1 o 

, 1 
A|’ u ol (5“t‘ )]d5* + / 0(5', x) x 

£ 

A i' [ 7 U oi u '" 0]dt ' } (lla) 

where 

N - (PL 2 /k) K. - (Pl 2 /k) " SdT) (de,, /dT). 

1 1-1 1 “l T 1 


H 11 ■ A o1 (T) (de o,1 /dT)T l • H 1 “ ^ H 11 

The dimensionless gas property N characterizes 
the relative importance of radiation versus 
conduction within the gas under optically thin 
conditions. Also, by combining Eqs. (5) and (10) 
another form of the transient energy equation Is 
obtained as 


v - 9 x - 12 <^ 2 > 


■2 (L/k) Z H., u . {/ 00/85' ) x 
1 1-1 11 01 o 


A 1 [| u 0 i(^')]d5’ - / 00/85') x 
£ 


A* t| u©i (5'-5)] d 5‘) (Hb) 

Note again that Eq. (lib) can be obtained 
directly by Integrating the right-hand side of Eq 
(11a) by parts. Quite often, Eq. (lib) Is the 
convenient form to use In radiative transfer 
analyses. 

For flow problems, the quantity of primary 
Interest Is the bulk temperature of the gas, 
which may be expressed as [11] 


3 


(16) 


1 ? 

e b * <T b' T l )/<q w L/k) * 6 / S(5,t)(e-5 )d5 (12) 

The heat transfer q w Is given by the expression, 
q w * h c (Tj-Tjj), where h c Is the convective heat 
transfer coefficient (W/cm 2 -K). In general, the 
heat transfer results are expressed In terms of 
the Nusselt number Nu - h c D h /k. Here, D h 
represents the hydraulic diameter, and for the 
parallel plate geometry It equals twice the plate 
separation, l.e., D h * 2L. Upon eliminating the 
convective heat transfer coefficient h c from the 
expressions for and Nu, a relation between the 
Nusselt number and the bulk temperature is 

obtained as 


a; [| u o1 U'-0]d5'> 

It should be pointed out that by combining Eqs. 
(9) and (15) another useful form of the energy 
equation can be obtained for the steady case. 

For the case of negligible radiation, Eq. 
(16) reduces to a very simple form and utilizing 
the boundary condition 6(0) * 0, the solution 
of the resulting equation Is found to be 


0(0 ■ £(2 5 2 - t 3 - 1) (17a) 


The result for the bulk temperature is found by 
combining Eqs. (12) and (17a) as 


Nu = 2 Lq w /k(T r T b ) - -2/0 b (13) 


9 b - 17/70 


(17b) 


The heat transfer results, therefore, can be 
expressed either In terms of Nu or 6^. 

limiting Cases and Solutions 

Before discussing the method of solution for 
the general case, It Is advisable to explore the 
various limiting cases. Quite often, closed form 
solutions can be obtained for some of these 
cases. Specifically four limiting cases are 
considered here and these are the steady laminar 
flows, the case of negligible radiation, the 
optically thin limit, and the large path length 
limit. 

Steady Laminar Flow 

For steady-state conditions bQ/bx * 0 and 
Eqs. (11) provide two forms of the energy 
equation for this case. Another convenient form 
Is obtained by letting d6/ax « Q In Eq. (5) such 
that 

6“ - 12 <£-£ 2 ) * (1/qJ dq R /d£ (14) 

By integrating Eq. (14) once and using the 
conditions that at 5 * 1/2, q R U) and 
(d8/d£) are equal to zero, one obtains 

e' - 2(3£ 2 - 2C 3 ) + 1 « (15) 

A combination of Eqs. (8) and (15) results In 
8‘ - 2(3£ 2 - 2C 3 ) + 1 


, n £ 

* , (L/k ) £ H,. u ni {/ 6(5') x 

1 1«1 11 01 0 

, 1 
Aj [| u^U-OJdC* - J e(c') X 


The results provided by Eqs. (17) are useful In 
determining the extent of radlatve contributions. 

Negligible Radiation 

For the case of negligible radiation, N ■ 0 
and both forms of Eq. (11) reduce to 

- e t - 12 (c-c 2 ) (18) 

By employing the product solution procedure, the 
solution of Eq. (16) can be obtained and the 
result can be expressed In terms of the bulk 
temperature through use of Eq. (12). 

Alternately, the solution of Eq. (18) Is assumed 
to be of the form 


* g(0 + h( 5 ,t) (19) 


From Eqs. (18) and (19), there is obtained two 
separate equations as 

g" - 12(£-£ 2 ) (20) 


- 0 ( 21 ) 

Solutions of Eqs. (20) and (21) are found 
separately from which the complete solution for 
the temperature distribution Is obtained as [25] 

e(c.t) « c(2 c 2 - c 3 - i) 

" 2 
♦ £ C n sin (a l) exp(-a x); a - nx 

" ml (22a) 

where 


C - (4/a 5 ) [ (12 - 12a 2 + a 4 ) cos(a) 
n 
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- 24], n = 1. 2.... 


where 


Finally, the expression for the bulk temperature 
Is found by using Eq. (12) as 

CD 

9. = -17/70 + 6 £ C [(1/a) 
b n-1 n 

+ (4/a 3 )] exp(-a 2 x) (22b) 


C n'° 

• 32 (3N + a 2 )/(3N 2 a 3 ) 


» for n even (28b) 


+ 2a/[3N 2 (3N + a 2 )] 


, for n odd (28c) 


Equations (22) are useful In determining the 
extent of radiative contributions for the 
transient case. 

Optically Thin Limit 

In the optically thin limit, the expression 
for the bulk temperature for the steady case Is 
found to be [11, 25] 

9. = G 1/ ( 3N ) 3 ] {576(3N)' 1/2 (NEXP) - 21. 6N 2 
D 

+ 72N - 288} (23) 


where 


NEXP - {1-exp [-(3N 1/2 ]> /{ 1 + exp[-(3N) 1/2 ]} 


Both forms of the transient energy equation, 
Eqs. (11a) and (lib), reduce to a simplified form 
in the optically thin limit as [25] 

9 „ - 9^ - 3N0 = 12U-E 2 ) (24) 

Assuming a solution of the form given by Eq. 
(19), Eq. (24) Is written as 

- h t - 3Nh * - g^ + 3Ng «■ 12(5 - 5 ? ) (25) 

Consequently, 

g" - 3Ng ■ 12(5 - 5 2 ) (26) 

and h - h - 3Nh ■ 0 (27) 

From the solution of Eqs. (26) and (27), the 
solution for Eq. (24) Is obtained as 


By combining Eqs. (12) and (28), the expression 
for the bulk temperature Is obtained as 

9 b ■ 6 <{l6/3N 2 )[s(nh(-/3R/2)/s1nh(/3JT)] x 

[ ( 1/3N) cosh(/3ff/2) - (4 + /3IT) (3N)' 3/2 x 

s1nh(/3IT/2)] + (4/N) [-1/30 + 1/(9N) ] 

» 

♦ £ C (1/a + 4/a 3 ) exp[-(3N + a 2 )x]} (29) 

n*l " 

Large Path length Limit 

In this limit, the steady-state energy 
equation, Eq. (16), reduces to [8, 11, 25] 

? i 1 

0 ' - 2(35 - 25 ) + 1 * M / 9 ( 5 ' ) d 5'/(5 - 5 ') 

0 (30) 

where 


n 

M - H. L/k * (L/k) z Am (de ,/dT ) T 
1 ol wl Tj 

The nondlmenslonal parameter M constitutes the 
radiation-conduction Interaction parameter for 
the large path length limit. Equation (30) does 
not appear to possess a closed form solution; a 
numerical solution, however, can be obtained 
easily. 

In the large path length limit, the 
transient energy equations, Eqs. (11a) and (lib), 
reduce to 


6^ - 0 x - 3M0 - 12(5 - 5 2 ) 



0(5*. x) d5'/(5 - 5* ) 2 


0(5, x) = ( 16/3N 2 ) [slnh ( -/W2)/s)nh(/3tf) ] 
cosh [/3JT (5-1/2)] + (4/N)(5 2 - 5 + 2/3N) 

* 7 

+ E C n sin (aFj exp [ - ( 3N + a ) x]; a * nn 
n=1 (28a) 


♦ / 0(5* ,x)d5'/(5‘ -5) 2 ] 
5 

- 9 t - 12 (5 * 5 2 ) 


Ola) 


5 


1 

* m / (3e/ae , )d57(e-c‘) uib) 

o 

Since |(5-C*) 2 |-|tt'-C) 2 |. Eq. (31a) can be 
written as 


0 U - e t - 3N8 - 12(5 * C 2 ) 

l , 

* - M / 9(5'.x) <157(5 - 5'r (31c) 

0 

Through integration by parts, Eq. (31c) can be 
expressed as 


- e t - 3N6 - 12(5 - 5 2 ) 

l 

= m / (ae/ac'jds 7(5-5') Old) 

o 

Equations Ola) - Old) represent different forms 
of the governing equations in the large path 
length limit. With the exception of the term 
(-3N0) on the left-hand side, Eq. (31d) Is 
identical to Eq. (31b). Since H represents the 
radiation-conduction Interaction parameter only 
in the optically thin limit [8], it should not 
appear In the governing equation for the large 
path length limit. Thus, Eq. (31b) is the 
correct equation to use for solution In the large 
path length limit; the solution of this equation 
is obtained by numerical techniques. 

Method of Solution 

The solution procedures for both steady and 
unsteady cases are presented in this section. In 
principle, the same numerical procedure applies 
to both the general and large path length limit 
cases. 

Steady-State Solutions 

The general solution of Eq. (16) is obtained 
numerically by employing the method of variation 
of parameters. For this, a polynomial form for 
9(C) Is assumed In powers of £ as 

e( 5 ) - £ a m 5 m (32) 

m*0 ■ 

By considering a five term series solution (a 
quartic solution in £) and satisfying the 
boundary conditions 9(0) * 9' (1/2) * 0 and 

9'(0) * -9* ( 1 ) , one obtains 


8(5) « 4j(5 - 25 3 + 5 4 ) + a 2 (5 2 - 25 3 + 5 4 ) (33) 

9* (5) * Sjd - 65 2 + 45 3 ) + a 2 C 25 - 65 2 + 45 3 ) 

(34) 

A substitution of Eq. (34) in Eq. (16) results in 


aj(l - 65 2 + 45 3 ) + a 25 * 65 2 + 45 3 ) 


- 2(35 2 - 25 3 ) + 1 - | <L/k) Hjj u q1 x 


5 , 

</ 8(5') a; [| u ol (5 - 5')]d5* 


I , 

- / 8(5') [| u of (5' - 5)] d5‘) (35) 

where expressions for 9(5') are obtained from 
Eq. (33). 

The two unknown constants aj and a 2 In Eq. 
(35) are evaluated by satisfying the Integral 
equation at two convenient locations (5«0 and 
5-1/4 In the present case). The entire 
procedure for obtaining aj and a 2 Is described In 
[25], With known values of aj and a 2 , Eq. (33) 
provides the general solution for 9(5). The 
expression for the bulk temperature Is obtained 
by combining Eqs. (12) and (33) as 


8 b » (1/70) (17a j + 3a 2 ) (36) 

It should be noted that for the case of no 
radiative Interaction a 2 » 0 and aj « -1, and Eq. 
(36) gives the result of Eq. (17). 

The governing equation for the large path 
length limit Is Eq. (30). For this equation also 
the solution Is given by Eqs. (33) and (36) but 
the values of a's are completely different In 
this case [25]. 

Transient Solutions 

The governing energy equations for the 
transient case are Eqs. (11a) and (11b). General 
solutions of these equations are obtained also 
numerically by employing the method of variation 
of parameters. For the present problem, a 
polynomial form for 0(5, t) Is assumed as 


n 


8(5. t) - £ a (x) 5 m 

m*0 m 

(37) 

For a quadratic temperature distribution In 5 
(with time dependent coefficients), Eq. (37) Is 
written as 

8(5. x) » a 0 (x) ♦ aj(x) 5 + a 2 (x) 5 2 

(38a) 

By using the boundary conditions 9(0, x) 
8^ ( 5*1/2) « 0, this reduces to 

* 0 and 

8(5. x) - g(x) (5-5 2 ) 

(38b) 
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where g(t) represents the time dependent 
coefficient. Consequently, 


B U t x) = g ( x ) (1-2^); 

* -2gU); 9 T (5.t) * (e-E 2 ) g’(t) (39) 

Also, a combination of Eq. ( 6 a) and (38b) yields 
the Initial condition 


e(S, 0 ) - g( 0 ) - 0 (40) 


Note that essential boundary conditions are used 
already In obtaining the solution represented by 
Eq. (38b). 

8 y employing Eqs. (38b) and (39), Eqs. (11a) 
and (lib) are transformed In alternate forms 
which are expressed in a compact form as 


VO, 

g*(x) + g(x) + 12 = 0 (41) 


where Jj(C) and VO are defined in [25]. The 
function VO Is used for solution of Eq. (11a) 
and * s usec * ^ or solution of Eq. (lib). The 

solution of Eq. (41) satisfying the Initial 
conditions of Eq. (40) Is given by 

g(x ) = j|^ j (exp [“J(Ox]~l} (42) 

The temperature distribution given by Eq. (38b) 
can be expressed now as 

0 (£,x) = y ( exp[- J( £ ) x ]”1) ( ) (43) 

The expression for the bulk temperature is 
obtained through use of Eq. (12) as 

6 . = 72 / [(C-5 2 ) 2 /J(5)](exp[-J(5)T]-l) (44) 

D 0 

Note that In Eqs. (42)-(44), J(E) becomes J,(E) 
for solution of Eq. (11a) and J-(E) for solution 
of Eq. (lib). 

For a quartlc solution In E. Eq. (37) gives 
the result similar to Eq. (33) which for the 
transient case Is expressed as 

9(E.x) * g(-c) (E - 2e 3 + E 4 ) 

+ h(x) (l 2 - If} + E 4 ) (45) 


By substituting Eq. (45) Into Eq. (11a), one 


obtains 


x g'(x) + O^U) + yb'(x) + J 4 U) h(x) ■ - z 

(46) 

where 

X » (E - 2E 3 + E 4 ); y * (E 2 - 2E 3 + E 4 ) ; 
z = 12(E - E 2 ) 


and functions J 3 U) an< * are defined In 

[25]. Equation (46) constitutes one equation In 
two unknowns, namely g(x) and h(x). However, 
since the equation Is linear In x, the principle 
of superposition can be used to split the 
solution Into two equations as 


x g* (x) + J 3 U) * -z/2 (47) 

yh'(x) + J 4 (0 h(x) * -z/2 (48) 

The Initial condition for this case can be 
written as 

9(E,0) - g(0) (E - 2E + E 3 ) 

+ h( 0 ) (E 2 - 2E 3 ♦ E 4 ) * 0 (49a) 


Consequently, 


g( 0) « 0; h( 0) - 0 (49b) 


The solution of Eqs. (47) and (48) 
satisfying the appropriate Initial condition of 
Eq. (49b) Is given respectively as 


g(x) « [zU)/2J 3 U)] {exp[-J 3 U)x/xU)]-l} (50) 

h ( x) - [zU)/2J 4 U)3 (exp[-J 4 (^)x/y(^)]-l} (51) 

By substituting Eqs. (50) and (51) Into Eq. (45), 
the expression for the temperature distribution 
1 $ obtained as 

9(E.x) - [6(E-E 2 )(E-2E 3 H 4 )/J 3 (E)] x 
(exp[-J 3 (E)T/x(5)]-l> + [6(E-E 2 ) x 
(E 2 -2E 3 +E 4 )/J 4 (E)] (exp[-J 4 (E)x/y(E)]-l) (52) 

The bulk temperature In this case Is given by 
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e h = 36 / [(e-5 2 )U-25 3+ 5 4 )/>M£)] X 

0 0 J 

1 ? 

(exp[-J (5)x/x(5)]-l)d5 + 36 / [(5-5 ) x 
J 0 

(5 Z -25 3 H 4 )/J 4 (5)](exp[-J 4 (5)x/y(5)]-l)d5. (53) 

where x and y are defined in Eq. (46). 

By substituting Eq. (45) Into Eq. (lib), 
there is obtained 

xg‘ + J 5 U) gU) + yh + J 6 (t) h(x) - -z (54) 


employed and this is given by 


A(u) » In (uf(p) [ ~ gf(g) J + 1> (55) 

where 

f(0) * 2.94[1 - exp(- 2.60 0)] 


and 0 represents the line structure parameter. 

The spectral information and correlation 
quantities needed for different species were 
obtained from Refs. 7-9, 26 and 28. The specific 
VR bands considered for each species are: CO 

(4. 7p fundamental), OH (2.8m fundamental). C0 2 
(15m. 4.3m and 2,7m), and H 2 0 (20m, 6.3m, 
2.7m, 1.87m, and 1.38m). 


where again x.y.z are defined In Eq. (46) and 
functions J^U) and Jg(£) are defined in 
[25]. The solution procedure for this equation 
is Identical to that for Eq. (46) and the results 
for temperature distribution and bulk temperature 
are given respectively by Eqs. (52) and (53) with 
J 3 replaced by J 5 and J 4 by J 6 . 

In the large path length limit, the two 
applicable governing equations are Eqs. (31b) and 
( 3 Id ) . The solutions of these equations can be 

obtained from the general solutions by evaluating 
the Integrals in J function In the large path 
length limit [25]. 

Physical Conditions and Data Source 

As discussed In [25-27], four specific 

absorbing emitting species were selected for an 
extensive study; these are CO, C0 2 , OH and H 2 0. 
The species CO was selected because It contains 
only one fundamental vlbra tlon-rota tlon (VR) band 
and all spectral Information are easily available 
In the literature. It is a very convenient gas 
to test the numerical procedure without requiring 
excessive computatlnal resources. Species OH and 
H 2 0 are the primary radiation partlclpa ting 
species for the pressure and temperature range 
anticipated in the combustor of the scramjet 
engine. Species C0 2 , and combinations of C0 2 and 
H 2 0 are Important absorbing-emitting species In 
many other combustion processes. Different 
mixtures of various species (such as C0 2 + H 2 0, 
OH + HoO, and C0 Z ♦ H 2 0 + OH) were selected for 
parametric studies. Triermophyslcal properties of 
these species are given In [26] for different 
tempera tures. 

In radiative transfer analyses. It Is 
essential to employ a suitable model to represent 
the absorption-emission characteristics of 
specific species under Investigation. Several 
line-by-line (LBL), narrow-band, and wlde-band 
models are available to model the absorption of a 
VR band [7-12]. However, It Is often desirable 
to use a simple correlation to represent the 
total absorption of a wide band. Several such 
correlations are available In the literature [7- 
12). The relative merits of these correlations 
are discussed In [12]. In this study, the 
correlation proposed by Tien and Lowder [7] Is 


In a mixture of several species, spectral 
lines and bands overlap In certain spectral 
regions. The total absorptance In such regions 
cannot be calculated simply by adding the 
contributions of different bands and corrections 
should be made to account for the partial 
overlapping. If line-by-line or narrow band 
models are employed In the general formulation of 
the physical problem, then there Is no need for 
such corrections [12, 29, 30], The solution of 
LBL formulation, however, requires considerably 
large computational resources. Use of narrow 
band models offers some computational relief but 
certain spectral Information needed are not 
available for many species for temperatures 
higher than about 600 K. A relatively easier 
procedure (called the block method) Is suggested 
by Edwards [13] and Is useful in calculating the 
total emlsslvlty of a mixture of several species. 
Another method suggested by Penner and Varanasi 
[31] Is probably the most convenient method to 
use In the frame work of the radiative flux 
formulatln expressed In terms of the wld-band 
model absorptance and correlations. 


For a homogeneous path, the total 
absorptance of a band Is given by 

0» 

A(y) * / [1 - exp(- k y)] du (56) 

o 

where both k and <*> have units of cm" 1 . If In 

U) 

a spectral range Aw^, there are contributions 
from bands of different species, then for a 
homogeneous path the transmittance Is given by 


T 


Auk 


N 

exp(- E k y) 

J*1 U 1J 


where N represents the number of participating 
species In the gaseous mixture. Consequently, 
Eq. (56) can be expressed as 


A. » / [1 - exp(- E k y)] d<J (58) 

1 Au < j-1 u 1j 


If two bands of different (or same) 


species are 
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occupying approximately the same spectral range 
Au>j, then Eq. (58) reduces to 


A. * / {1 - exp(ic + r )y]} du (59) 

1 “ll “l2 


By employing the relations for the exponentials, 
Eqs. (59) can be expressed as [31] 


A. » / [1 - exp( - < y) ] du 

1 Auj “11 

+ J [l - exp( - < )] du> - 6A (60a) 

6a), “1Z 

or 

A * A i + A 2 - 6A (60b) 


where 

6A = / {[1 - exp(- k y)] x 

6 ( 0 , “11 


[1 - exp(- 



y) ] ( do) 


(60c) 


Use of Eqs. (60) has been made by Felske and Tien 
[32] to calculate absorptances of homogeneous and 
non homogeneous mixtures of CO 2 and H?0 In the 2.7 
p region for different pressure and temperature 
conditions. A similar procedure Is used In this 
study to account for the overlapping effects of 
different species [28]. 


For the physical problem considered, the 
dependent variables are 9 and 9^ (or Nu) and 
Independent variables are x and 1 !. The para- 
meters, In general, are T^, P, and L. The large 
path length and optically thin limits are 
characterized respectively by parameters M and N. 
Radiative and thermophysical properties of 
participating species are evaluated at different 
specified pressures and temperatures. 

Results and Discussions 


Extensive results have been obtained for 
variation of 9 and 0. for different conditions. 
Selected results are 0 presented here to compare 
solutions of quadratic and quartlc formulations 
and demonstrate the variation of 9 with x and l 
and of 0 h with t for single component systems 
and homogAeous mixtures. Some other results are 
available In [28]. For all results presented 
here, a lower value of 0 (or 9.) In the figures 
Indicates a higher value of temperature In the 
medium; this. In turn. Implies a relatively 
higher ability of the gas to transfer radiative 
energy. 

Results of quadratic and quartlc formula- 
tions are compared In Figs. 2-5 for different 
species. The centerline temperature variations 
with nondlmenslonal time are compared In Fig. 2 
for P * 1 atm, T w * 500 K, and L * 5 cm. The 


results show that the steady-state conditions are 
reached at an earlier time for H 2 0; and this Is 
followed respectively by C0 2 , CO, and OH. 
Significant differences are noted between the 
quadratic and quartlc solutions for larger 
times. For the conditions of results presented 
In the figure, the difference are found to be 
greatest for OH and lowest for HoO. The results 
for 0 versus £ are compared for P * 1 atm, L 

■ 10 cm, and T w * 500 K In Fig. 3 and T w * 1,000 

K In Fig. 4. The results demonstrate that 

considerable differences In solutions can occur 
at different locations In the channel and that 
the differences are larger for the lower wall 
temperature. Results for 0. versus L presented 
In Fig. 5 show that the differences In two 
solutions are relatively larger at lower plate 
spaclngs and that quartlc solutions approach the 
correct limiting solution for the case of no 

radiative Interaction (0^ * -0.243). From the 
results presented In Figs. 2-5 and In Refs. 26 

and 27 It Is concluded that while quadratic and 
quartlc solutions are Identical for radiative 
equilibrium and radiation with conduction cases, 
they differ significantly for the case of 
combined conduction, convection and radiation. 
As such, all other results for this study were 
obtained by using the quartlc formulation. 

The results for temperature variations with 
nondlmenslonal time are presented In Figs. 6-9 
for various species and for different physical 
conditions. The centerline temperature 

distribution for general and limiting cases are 
Illustrated In Fig. 6 for T w - 500 K, P * 1 atm, 
and L ■ 5 cm. The results show that the steady- 
state conditions are reached at about x * 0.5 
for all species. As noted earlier, these results 
also demonstrate that H 2 0 Is a highly radiation 
participating gas as compared to C0 2 , CO, and 
OH. In comparison to other species considered, 
OH takes relatively longer times to reach the 
steady state and Is least effective In 
transferring the radiative energy. For the 
specified physical conditions, the large path 
length solutions are closer to the general 
solutions and optically thin solutions provide 
higher rate of energy transfer. For all species, 
optically thin solutions reach the steady state 

faster than other solutions. Also, In the 
optically thin limit, C0 2 Is more effective In 

the radiative transfer process than other 
species. The reasons for such trends are given 
in [8, 26, 27]. The results for T w » 500 K and 
1,000 K are compared In Fig. 7, and they simply 
Indicate that the rate of energy transfer is 
higher at the higher temperature. The results 
for l * 0.25 and 0.5 are compared In Fig. 8 for 

T w * 1,000 K, P * 1 atm, and L - 5 cm. It Is 

noted that the rate of energy transfer Is higher 
at earlier times and at locations closer to the 
wall. The centerline temperature variations for 
L * 5 cm and 10 cm are Illustrated In Fig. 9 for 
T w * 1,000 K and P * 1 atm. As would be 
expected, the rate of energy transfer Is seen to 
Increase with the Increasing path length. 

The results for temperature variations 
within the plates are presented in Figs. 10-14 
for different conditions. Since the temperature 
profiles are symmetric, most of these results are 
Illustrated only for £ * 0 to £ * 0.5. The 
stady-state solutions for various species are 
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compared in Fig. 10 for P 3 1 atm* L * 10 cm, and 
different wall temperatures. Results again 
demonstrate the relative importance of different 
species for energy transfer. It is noted that 
the results for C0 2 at * 500 K are exactly the 
same as for OH at T • 1,000 K, and the results 
for CO and H 2 0 are about the same at T w 3 1,000 K 
and 2,000 K. This obviously Is a coincidence for 
the physical case considered. The general and 
limiting solutions for OH, CO, C0 2 , and H 2 0 are 
presented, respectively, in Figs. 11 through 14a 
for T w 3 500 K, P = 1 atm, and L 3 10 cm. These 
results also show that the steady-state condi- 
tions reach earlier for H 2 0 and C0 2 than for CO 
and OH. In general, the differences between 
limiting and general solutions are found to be 
small at earlier times; the naxlmum difference 
occurs at the steady-state conditions. In each 
case, the large path length solutions are closer 
to the general solutions, but the optically thin 
solutions are found to deviate considerably. 
This Is because for P 3 1 atm and L * 10 cm, the 
pressure path length is sufficiently high and the 
optically thin limit is not the correct limit for 
the physical case considered. It should be noted 
that for CO, C0 2 , and H 2 0, the optically thin 
solutions are identical for x > 0.05. Additional 
results given In [28] reveal that the differences 
between general and large path length solutions 
are Insignificant at higher temperatures for all 
species. The general solutions for the 
temperature variation across the entire duct are 
illustrated in Fig. 14b for H 2 Q for the physical 
conditions of T w 3 1,000 K, P * 1 atm, and L = 5 
cm. The parabolic nature of the transient 
profiles Is clearly evident and, In this case, 
the steady-state is reached at x > 0.5. 

The bulk temperature results as a function 
of the distance between the plates are presented 
In Figs. 15-19 for different times. General as 

well as limiting solutions are illustrated in 
these figures. As mentioned earlier, the 

limiting value of 9. = - 0.243 corresponds to 

negligible radiation. For all species, the 
results presented In the figures are for a 

pressure of one atmosphere. However, for any 
particular L, the large path length results 
essentially represent the limiting solutions for 
high pressures. The results. In general, demon- 
strate that the rate of energy transfer is higher 
at earlier times, the effect of radiation 
increases with Increasing plate spacing, and the 
radiative transfer Is more pronounced at the 
higher wall temperature. 

General as well as limiting solutions for 
the bulk temperature are Illustrated in Figs. 15- 
18 for individual species. It Is seen that for 
all species the general solutions for x = 0.5 
and 1.0 are essentially the same for all plate 
spaclngs, and the large path length results are 
valid for spaclngs greater than L * 10 cm for all 
times. The results for OH are presented In Figs. 
15a and 15b for T w » 500 K and 1,000 K, 
respectively. It is noted that optically thin 

results provide the correct limiting solutions 
for plate spaclngs upto L * 3 cm for all times. 
The results for T w * 500 K show only slight 
difference between general and large path length 
solutions for x = 0.5, and no significant 
difference was noted at earlier times (Fig. 15a). 
This, however, is not the case for the results 


presented in Fig. 15b at T w = 1,000 K. This 
trend In results for OH was noted also in [26, 
27], The results for CO are Illustrated in Fig. 
16 for T w = 500 K and P * 1 atm. In this case, 
the optically thin solutions are seen to be valid 
only upto L * 1.5 cm. For T w = 1,000 K and P 3 1 
atm, the results presented In Fig. 17 for C0 2 and 
In Fig. 18 for H 2 0 show the same general trend 
but the extent of radiative interaction Is 
entirely different. For C0 2 , there Is a 
considerable difference In general and optically 
thin solutions for all times. For H 2 0, however, 
the optically thin results are closer to the 
general solutions for spaclngs upto L = 2.5 cm. 
The results presented In Figs. 15-18 clearly 
reveal that for a fixed spacing between the 

plates the rate of radiative heating will be 

considerably higher at earlier times than at the 
steady state. Thus, In a particular physical 
system, the extent of radiative heating can be 
very Intense during the Initial stages of opera- 
tion. It Is also Important to note that both the 
optically thin and large path length results 
overestimate the Influence of radiation. Since 

these solutions can be obtained with less 

numerical complications, they can be utilized to 
assess whether or not, for a given gas, the 
Interaction of radiation Is going to be 

Important, 

A comparison of the general band absorptance 
results for the four gases Is shown In Fig. 19 
for a pressure of one atmosphere and a wall 
temperature of 1,000 K. The results clearly 

demonstrate the relative ability of the four 
species for radiative transfer at different path 
lengths. For a plate spacing of greater than L = 
3.0 cm, the results show the same trend as noted 
In Figs. 2-10. For lower plate spaclngs and 

relatively higher tempera tures, however, C0 2 
shows a significantly higher ability for radia- 
tive transfer than other species. This is a 
typical distinguishing feature of the C0 2 under 
optically thin conditions [8, 26, 27]. 

For steady-state conditions, bulk tempera- 
ture results were obtained for mixtures of 
different absorbing-emitting species under 
various conditions and some of these are 
presented in Figs. 20-27. The relative amount of 
each species In the mixture and the relative 
ability of the species for radiative transfer In 
a given physical condition determine the extent 
of radiative interaction. 

As mentioned before. In many combustion 
processes involving fossil fuels, the predominant 
products of combustion are carbon dioxide and 
water vapor. Extensive studies have been 
conducted in the literature to determine the 
total emissl vlty of C0 2 and H 2 0 for homogeneous 
and non homogeneous conditions. However, only 
limited studies are available involving mixtures 
of C0 2 and H 2 0 for realistic and Important 
physical conditions. The bulk temperature 
results for three different mixtures of C0 2 + H 2 0 
are presented in Figs. 20-23 for different 
temperatures and pressures. As would be 
expected, the results. In general, show that the 
extent of radiative Interaction Increases with 
Increasing temperature, pressure, and path 
length. It Is seen that the radiative Interac- 
tion Is stronger for higher amount of H 2 0 in the 
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mixture. However, for the optically thin 
conditions, the radiative contribution is seen to 
Increase with Increasing amount of CO 2 at 
relatively higher temperatures. 

The bulk temperatures results for mixtures 
of OH and H 2 O are Illustrated In Figs. 24-26 and 
for mixtures of C0 2 , H^O, and OH In Fig. 27 for 
different pressures ana temperatures. All these 
results clearly demonstrate that the radiative 
ability of a gaseous mixture essentially depends 
on the amount of highly radiation participating 
species In the mixture. For example, a 

comparison of results presented In Figs. 24 and 
27 for P * 1 atm Indicates that the rate of 
radiative transfer Is significantly higher with 
the Inclusion of 20? C0 2 In the mixture of OH and 
HoO. It Is further noted that OH becomes a 
highly radiation participating species at higher 
temperatures and pressures. 

Concluding Remarks 

Analytical formulations and numerical 
procedures have been developed to Investigate the 
transient radiative Interaction of absorbing- 
emitting species in laminar fully-developed 
flows. Extensive results have been obtained for 
OH, CO, CO 2 and HoO for different physical 
conditions. Illustrative results for the 
temperature distribution and bulk temperature are 
presented for different pressures and wall 
temperatures. In these results, a lower value of 
temperature Implies a relatively higher ability 
of the gas to transfer radiative energy. 

Comparative results of quadratic and quartlc 
formulations confirm the need to use the quartlc 
formulation In the numerical procedures for the 
case of combined conduction, convection, and 
radiation. The results. In general, demonstrate 
that the steady-state conditions are reached at 
about t * 0.5 for all species, H 2 O Is a highly 
radiation participating species (as compared to 
C0 2 , CO, and OH), the rate of energy transfer Is 
higher at earlier times and at locations closer 
to the wall, differences between the limiting and 
general solutions are small at earlier times, the 
effect of radiation increases with Increasing 
plate spacing, and the radiative transfer Is more 
prounounced at higher wall temperature and 
pressure. Similar conclusions can be drawn from 
the results presented for various mixtures. The 
results clearly show that for a given physical 
condition the radiative ability of a gaseous 
mixture depends essentially on the amount of 
highly radiating species In the mixture. From 
the results presented In this study, the extent 
of total heating can be determined for different 
times. This Information Is essential in design- 
ing thermal protection systems for operations 
during the Initial stages of Intense heating. 
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q w = CONST. 

Fig. 1 Physical model and coordinate 
systems for flow of radiating gases between 
parallel plates. 
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Fig. 4 Comparison of quadratic and qua r tic 
results for t e mp e rature variations across 
the duct; P=1 atm, 1-10 cm, and T*»1,000 K. 


Fig. 2 Comparison of quadratic and quartlc 
results for the centerline temperature 
variations with time; 1**500 K, P*1 atm. and 
L=5 cm. 



Fig. 3 Comparison of quadratic and quartlc 
results for temperature variations across 
the duct; P-1 atm, L*10 cm, and T**500 K. 



Fig 5 Comparison of quadratic and quartlc 
results for the bulk temperature; P-1 atm 
and T*-500 K. 
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Fig. 6 Coaparlson of general and Halting 
results for the centerline temperature 
variations with tlae; T w *500 K, P-1 ata, and 
L-5 ca. 


Fig. a Coaparlson of teaperature variations 
with tlaa for t-0.25 and 0.5; T w -1.000 K, 
P-1 ata, and L-5 ca. w 
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Fig. 7 Coaparlson of centerline teaperature 
variations with tlae for T w «500 K and 
1.000 K; P-1 ata and L-5 ca. 


C *at/L' 

Fig. 9 Coaparlson of centerline teaperature 
variations with tlae for L-5 ca and 10 ca; 
V1.000 K and P-1 ata. 
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Fig. 10 Comparison of temperature 
variations across the duct for T w *500 K, 
1.000 K, and 2,000 K; P«1 atm and L=10 cm. 


Fig. 12 Comparison temperature variations 
across the duct for CO; T w *500 K» P-1 atm, 
and L-10 cm. 
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Fla 11 Comparison of temperature 
vari»tlons»cross the duct for OH; V500 K. 
p-1 atm, and L-10 cm. 


Fig. 13 Comparison of temperature 
variations across the duct for CO? 
V500 K, P-1 atm, and L-10 cm. 
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Fig. 14a Comparison of temperature 
variations across the duct for HjO; 
V500 K, P-1 atm. and L-10 cm. 


Fig. 15a Variation of bulk temperature with 
pltte spacing for OH; T>500 K and P-1 ata. 
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Fig. 14b Comparison of temperature 
variations across the duct for H 2 0; 
VI. 000 K. P-1 atm. and L-5 cm/ 


Fig. 15b Variation of bulk temperature with 
plate spacing for OH; VI, 000 K and 
P-1 atm. 
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Fig. 16 Variation of bulk temperature with 
plate spacing for CO; T “500 K and P*1 ata. 
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Fig. 17 Variation of bulk temperature with 
plate spacing for C0 2 ; 7**1,000 K and 
P-1 atm. 



Fig. 18 Variation of bulk temperature with 
plate spacing for 8*0; T*-1.000 K and 
P*1 atm. 
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Fig. 19 Comparison of bulk temperature 
results for T w -1,000 K and P=1 atm. 
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Fig. 20 Bulk teaperature results for 
C0 2 ♦ H 2 0 for P-0.1 at*. 
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Tig. 21 Bulk temperature results for 
C0 2 ♦ H 2 0 for atm. 
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Fig. 22 Bulk teaperature results for 
^®2 + **2° for F*10 ata. 



Fig. 23 Bulk teaperature results for 
C0 2 + H 2 0 for T w «2,000 K. 
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INCOMPRESSIBLE AND COMPRESSIBLE INTERNAL FLOWS 


S. N. Tiwari 1 , D. J. Singh 2 and P. A. TVivedi 3 
Old Dominion University, Norfolk, Virginia 23529-0247 

Abstract 

Analyses and numerical procedures are presented to investigate the radiative interactions of 
gray and nongray absorbing-emitting species between two parallel plates and in a circular tube. 
Laminar fully developed incompressible as well as entrance region subsonic flows are considered. 
The participating species considered are OH, CO, CO2, CH4, and H2O. Results obtained for 
different flow conditions indicate that the radiative interactions can be quite significant in fully 
developed incompressible flows. For subsonic flows, however, the flowfield is not changed 
significantly due to radiative interaction. 

Nomenclature 


A 

A 0 

Co 

C p 

eu> 

e u Jo 

ej ,e2 

k 

L 

P 

Pi 

Pt 

Pr 

ar 

q c 

qRu> 

r 

r 0 (Xl 


band absorptance = A(u,/?), cm 1 

band width parameter, cm' 1 

correlation parameter, atm' 1 - cm' 1 

specific heat at constant pressure kJ/kg-K = erg/gm-K 

Planck’s function, (W-cm'^/cm 1 

Planck’s function evaluated at wave number u> 0 

emissive power of surfaces with temperature Ti and T2, W-cm' 2 

thermal conductivity, erg/cm-sec-K 

distance between the plates 

pressure 

partial pressure 

stagnation pressure 

Prandtl number 

total radiative heat flux, J/m 2 -s 
conduction heat flux, J/m 2 -s 
spectral radiation heat flux, (w-cm'^/cm* 1 
physical coordinate for circular tube 
radius of the tube 
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Materials Inc., Hampton, Virginia. 

3 Graduate Research Assistant, Dept of Mechanical Engineering and Mechanics. 
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s 

integrated intensity of a wide band, atm' 1 -cm 

T 

temperature, K 

t,,t 2 

wall temperature, K; Ti = T w 

T b 

bulk temperature, K 

u 

nondimensional coordinate = SPy/A 0 

Uo 

nondimensional path length = SPL/A 0 

0 

flow angle 

Ku 

spectral absorption coefficient, cm' 1 

Kp 

Planck mean absorption coefficient, cm' 1 

Z’V 

computational coordinate 

p 

density, kg/m 

a 

Stefan-Boltzmann constant, erg/(sec-cm-K 4 ) 

to 

wave number, cm" 1 

W 0 ,U>i 

wave number at the band center, cm 1 


I. Introduction 

There is a renewed interest in investigating various aspects of radiative energy transfer in 
participating mediums. Radiative interactions become important in many engineering problems 
involving high temperature gases. Recent interest lies in the areas of design of high pressure 
combustion chambers and high enthalpy nozzles, entry and reentry phenomena, hypersonic 
propulsion, and defense oriented research. 

Basic formulations on radiative energy transfer in participating mediums are available in 
standard references [1-8]. The review articles presented in [9-15] are useful in understanding 
the radiative properties of participating species and the nature of nongray radiation. The validity 
of radiative transfer analyses depends upon the accuracy with which absorption-emission and 
scattering characteristics of participating species are modeled. There are several models available 
to represent the absorption-emission characteristics of molecular species and these are reviewed 
in [12, 13]. 

The purpose of this study is to present analyses and solution procedures for infrared 
radiative energy transfer in molecular gases when other modes of energy transfer simultaneously 
occur. Attention is directed to investigate radiative interactions in laminar incompressible and 
compressible internal flows between two parallel plates and in a circular tube. Radiative 
interactions in incompressible duct flows have been investigated extensively in the literature 
with certain inherent simplifying assumptions [16-35]. The main thrust of this research is to 
develop procedures for investigating radiative interactions in subsonic compressible entrance 
region flows and provide a parametric study for different participating species. However, basic 
formulations, numerical procedures, and certain results of incompressible fully-developed duct 
flows are also presented for comparative purposes. 
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II. Basic Formulation 

Basic governing equations and boundary conditions are provided here for the two physical 
problems considered in this study. These are fully developed laminar incompressible flow 
between two parallel plates and within a circular tube, and developing laminar compressible 
flow in ducts. 


A. Fully-Developed Incompressible Flow 

The physical problem considered is the energy transfer in laminar, incompressible, constant 
properties, fully developed flow of absorbing-emitting gases between two parallel plates (Fig. la) 
and within a circular tube (Fig. lb). The condition of uniform surface heat flux is assumed such 
that the surface temperature varies in the axial direction. Extensive treatment of this problem 
for the parallel plate geometry is available in the literature. However, only limited attentions 
have been directed in obtaining nongray solutions for the circular tube geometry [19, 20, 23, 
26], One of the objective of this study is to obtain extensive nongray solutions for the circular 
tube geometry for different flow and physical conditions. 

For the parallel plate geometry (Fig. la), the energy equation can be expressed as [1] 


dT dT _ d 2 T 1 dg R 
U dx + V dy ° dy 2 pC p dy 


( 1 ) 


where a = (k Jp Cp), represents the thermal diffusivity of the fluid. In deriving Eq. (1), it has 
been assumed that conduction as well as radiation heat transfer in the x direction are negligible 
compared to that transferred in the y direction. In addition, it has been assumed that the Eckert 
number of the flow is small. These assumptions are consistent with the formulations presented 
in [4, 6, 8, 10, 1 1]. 


The justification for a fully developed flow in the presence of radiative interaction is provided 
in [4, 8] and in cited references. For a fully developed flow, v = 0,and u is given by the well- 
known parabolic profile u = 6 u m (£ - £ 2 ), where ( = y/L and u m represents the mean fluid 
velocity. Also, for the flow of a perfect gas with uniform heat flux, dT /dx is constant and 
is given by 


dT /dx = (2aq w )/(u k) 


( 2 ) 


A combination of Eqs. (1) and (2), therefore, results ,n 


d? z q w d£ 


(3) 


where 

0 = (T — T 1 )/(q w L/k); T t = T w 

Equation (3) is the governing energy equation in nondimensional form for the parallel plate 
geometry. The boundary conditions for this problem can be expressed as 


6 = ( 0 ) = 0 ( 1 ) = 0 ; 0 ( (t = 1 / 2 ) = 0 ; 0 ( (( = 0 ) = - 0 <(£ = 1 ) 


(4) 
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It should be noted that all boundary conditions given in Eq. (4) are not independent; any two 
convenient conditions can be used to obtain specific solutions. By integrating Eq. (3) once and 
using the conditions that at £ = 1/2, q R (0 and (d0/d£) are equal to zero, one obtains 

de/d( - 2(3£ 2 - 2£ 3 ) + 1 = q R (0/qw (5) 

Either Eq. (3) or Eq. (5) can be used to investigate the effect of radiative interaction. 


For the circular tube geometry (Fig. lb), the energy equation, as given by Eq. (1), can 
be expressed as 


dT a d ( dT\ lid, , 

U dx r dr \ dr ) pC p r dr f<lR 

where u is given by the relation u = u m [1 — (r 0 ) 2 ]. For a uniform wall heat flux and 
developed heat transfer, in this case, dT/dx is given by 


( 6 ) 

fully 


dT/dx = (2aq w )/(u m r 0 /k) 
Consequently, Eq. (6) may be expressed as 



where 


(7) 

( 8 ) 


0 _ (T T w )/(q w r 0 /k),£ — r/r D 

It is important to note the difference in the definitions of 9 given in Eqs. (3) and (8). Upon 
integrating Eq. (8) and noting that dO/df, = 0 and q R = 0 at ( = 0, there is obtained 

de/di + £ 3 - 2£ = q R /q w (9) 

The boundary condition for this equation is given as 0(1) = 0. 


For incompressible flow problems, the quantity of primary interest is the bulk temperature 
of the gas. For a fully-developed flow between parallel plates, this is expressed as 

0 b = (T b -Ti)/(q w L/k) = 6 [ 0(()(Z-( 2 )d( (10) 

Jo 

where q w = h (Tj - T b ), and h represents the equivalent (or effective) heat transfer coefficient 
(W/cm 2 - K). 


For the fully-developed flow in a circular tube, the expression for the bulk temperature is 
given by 

«'b(Tb-T„)/(qr 0 /k) = 4 f 0«)(«-( 3 )d{ (11) 

where q w h (T w - T b ). 

In general, the heat transfer results are expressed in terms of the Nusselt number N u = h 
Dh/k. Here, D b represents the hydraulic diameter. For the parallel plate geometry, Dh equals 
twice the plate separation, i.e., D b = 2L. For the circular tube geometry, Dh represents the actual 
diameter of the tube, i.e., Dh = 2r 0 . Upon eliminating the heat transfer coefficient h from the 
expressions for q w and Nu, a relation between the Nusselt number and the bulk temperature 
is obtained for both geometries as Nu = - 2/0 b . The heat transfer results, therefore, can be 
expressed either in terms of Nu or 0 b . 
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B. Laminar Compressible Flow 

The physical problem considered for basic understanding of radiative interaction in com- 
pressible flows is two-dimensional variable property laminar flow between two parallel plates. 
For this model, two-dimensional Navier-Stokes equations in fully conservative form are used to 
describe the flow field. These equations, in physical domain, can be written as (Fig. lc) 


dU dF dG 
dt + dx dy 


= 0 


(12) 


where vectors U, F and G are expressed as 



’ p 


> 

u = 

pn 

, F = 

pn 2 + P - r xx 


pv 


pnv - T xy 


.PH-P. 


. (E + P)u - ur xx - vr xy + q Cx + q Rx 


P 


pnv - r xy 

PV 2 + P - Tyy 


,P = pRT 


L(E + P)v 

ur xy vr yy + qCy T qRy 


The viscous stress terms appearing in the definitions of F and G are given in [36]. The relations 
for conduction heat transfer in x and y directions are given by 


qcx = ~ k fa' qcy = (13) 

The terms qR X and qR y represent radiative fluxes in x and y directions, respectively; expressions 
for these are provided in the next section. The total energy flux in a given direction is given by 
the corresponding last term in the definition of F or G. Consequently, this formulation involves 
all kinds of energy interaction including frictional (aerodynamic) heating. The coefficient 
of viscosity is evaluated by using the Sutherland’s formula and the coefficient of thermal 
conductivity is calculated by using a constant value of the Prandtl number. 

Equation (12) can be used to obtain solutions for all kinds of compressible flows. However, 
boundary conditions and numerical procedures for different flows are quite different. For 
subsonic flows the treatment of the inflow conditions is guided by the theory of characteristic. 
A locally one-dimensional flow has four characteristic equations with slopes u, u+c, u and u-c. 
If the flow is subsonic at the inflow, then the u-c characteristic has a negative slope and it 
propagates informations from the interior upstream to the inflow boundary. In this case, only 
three quantities can be specified at the inflow and the fourth quantity must be allowed to vary 
as the solution progresses. 
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In this study, the stagnation pressure, stagnation temperature and flow angle are specified 
at the inflow. These quantities are related to the state pressure and state temperature by the 
following equations 


Pt 

P 

Tt 

T 



2 

7 - 


-y _ 1 \ T/(7-l) 

7 ^M 2 ) ; 

~M 2 ^; u = T an ^ IF ) 


(14) 


Equations (14) is a system of three equations in four unknowns, P, T, u and v. To complete 
the system of equations, a zero order extrapolation is used for the pressure at the inflow. The 
outflow boundary is also calculated based on the theory of characteristic. For subsonic flow 
at the outflow, the u-c characteristic propagates information upstream from the boundary to the 
interior, i.e., only one quantity can be specified at the outflow. The state pressure is specified 
at the outflow while u, v and T are calculated using a zeroth-order extrapolation. Along the 
surfaces, following boundary conditions were applied 


u = 0, v = 0, <9P/dy = 0, T = T w (15) 

The density is obtained from the equation of state using the computed surface pressure and 
prescribed surface temperature. 


III. Radiative Transfer Models 

Evaluation of the energy equation presented in Eqs. (3), (5), (8), (9) and (12) requires an 
appropriate expression for the net radiative flux. A suitable radiative transport model is needed 
to represent the true nature of participating species and transfer processes. In this section, a 
brief discussion of various absorption models is given and essential equations for the radiative 
flux are presented. 


A. Absorption Models 

Several models are available in the literature to represent the absorption-emission character- 
istics of molecular species. The total band absorptance of a vibration-rotation band is given by 


A = 



— exp (— « w X)d(u; — u; 0 )] 


(16) 


where is the volumetric absorption coefficient, ui is the wave number, u> 0 is the wave number 
at the band center, X = Py is the pressure path length, and the limits of integration are over the 
entire band pass. Various models are used to obtain the relation for A in Eq. (16). 

The gray gas model is probably the simplest model to employ in radiative transfer analyses. 
In this model, the absorption coefficients is assumed to be independent of frequency, i. e., k w 
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is not a function of u>. A convenient model to represent the average absorption coefficient of a 
gray gas is the Planck mean absorption coefficient « p which is defined as [1] 

= f k w ebu,(T)du;/eb(T) (17a) 

J O 

For a multiband gaseous system, this is expressed as 


n 

«p = [Pj/t^T" 1 )] 5Z e “'( T ) Si(T) < 17b > 


where Pj is the partial pressure of jth species in a gaseous mixture, e Wi (T) is the Planck function 
evaluated at the ith band center, and Sj(T) is the integrated band intensity of the ith band. 

Several other models for the mean absorption coefficient are available in the literature [1, 37]. 
Since these models account for detailed spectral information of molecular bands, this approach of 
radiative formulation is referred to as the “pseudo-gray formulation.” The gray gas formulation 
for radiative transport is very useful in parametric studies. 

There are several nongray models available in the literature to represent the absorption- 
emission characteristics of vibration-rotation bands. These are classified generally in four classes, 
(1) line-by-line (LBL) models, (2) narrow band models, (3) wide band models, and (4) band 
model correlations. A complete discussion on usefulness and application of these models is 
provided in [12, 13]. For many engineering applications, wide band model correlations provide 
quite accurate results. The most commonly used wide band model correlations are due to 
Edwards [11, 16] and Tien and Lowder [9]. The Tien and Lowder correlation for the total band 
absorptance is a continuous correlation and is given by the relation 

A(u, 0) = A(u, 0)/A o = In {uf(t)[(u + 2 )/(u + 2 f(t))] + 1 } (18) 


where 

f(t) = 2.94[1 - exp (— 2.60t)], t = /?/2 

and u = SX/A 0 is the nondimensional path length, 0 = is the line structure parameter, 7 

is the line half width, S is the integrated band intensity, and A 0 is the band width parameter. 
Equation (18) provides accurate results for pressures higher than 0.5 atmosphere [12, 13]. 

Spectral properties and correlation quantities for various radiation participating species are 
available in [5, 9, 11]. These are useful in gray as well as nongray radiative formulations. 
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B. Radiative Flux Equations 

For many engineering and astrophysical applications, the radiative transfer equations are 
formulated assuming one-dimensional planar systems. For diffuse nonreflecting boundaries and 
in absence of scattering, the expression for the spectral radiative flux is given by [1, 38] 


where 



exp 
k w exp 


-jMs-y) 
- y) 


dz 


dz 


(19) 


F w (z) ej^; F 2 u>(z) — e w (z) — e 2 o» 

It should be pointed out that the exponential kemal approximation has been used in obtaining 
Eq. (19). The total radiative flux in a given direction is expressed as 


fCO 

q R = / qR W dw 


( 20 ) 


A combination of Eqs. (19) and (20) provides a proper form of total radiative flux equation 
for obtaining nongray solutions of molecular species. Any convenient absorption model can be 
used to obtain nongray results. 


For a gray medium, Eq. (20) reduces to a simpler form and upon differentiating the resulting 
equation twice, the integrals are eliminated and there is obtained a nonhomogeneous ordinary 
differential equation as [1, 39] 


1 d 2 q R (y) 
k 2 dy 2 


9 3 de(y) 

- tor y = — —■ 

4 k dy 


(21) 


where k - « p . Equation (21) is a second order differential equation and, therefore, requires two 
boundary conditions. For nonblack diffuse surfaces, these are given as 


(^-j)[ qR (y)U-^ 


1 

dqR 

3 K 

. d y . 


= 0 


y=0 


fdqR 


dy 


= 0 


y=L 


(22a) 


(22b) 


Equation (21) along with boundary conditions can be used to obtain the energy equation for gray 
gas radiative interaction. In this formulation, there is no need to linearize the radiative interaction. 
For linearized radiation, however, T 4 - Tj = 4T 2 (T - TO and Eq. (21) can be written as 


d 2 qR/d £ 2 - -r 2 q R = 7 iqu , dO/d( 

where 

71 = 3r 2 /N; N = k k p /(4(tTi); r D = « P L 


( 23 ) 
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For black walls and Tj = T2, the boundary conditions for Eq. (21) become 


9r( 1/2) = 0; ~qR.(0) = — (dq R /d£)£ =0 

l T 0 


( 24 ) 


Equation (23) is useful in obtaining closed form solutions of simplified problems such as fully- 
developed flows. The parameters for this equation are ft and r Q ; ft characterizes the relative 
importance of radiation versus conduction for a gray gas. 

For a black circular tube, the spectral radiative heat flux is given by the expression [23] 


n 

/ F u ,(r')« u ,aexp 

/ t\ 

(r r ) 

Jo l 

. Jt sin 7 

cos 7 


dr' 


- J F w (r , )«: tJ aexp 

r° 

/ F w (r')«: a ,aexp 

Jt sin 7 


b/c. 


cos 7 
b/c w 
cos 7 


* V - r) 


dr' 


(r + r — 2r sin 7) dr 1 




d7 


(25) 


where F w (r') = e aJ (r') - eo;(T w ), and constants a and b have values of unity and 5/4, respectively. 
A combination of Eqs. (20) and (25) provides a convenient form of the total radiative flux for 
nongray analyses. 

For a gray medium, the expression for the total radiative flux can be obtained from differential 
approximation as [1, 23] 


d^ 

dr 


1 i, . 

7d7 (rqR ) 
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-/c 2 q R = 3 <r/e 


dT 4 

dr 


(26) 


For linearized radiation, dTVdr = 4 Tj, (dT/dr) and Eq. (26) is expressed as 


d 




9 _o dQ 

4 r 0 q R = 7 2 qw^- 


(27) 


where 

72 = 3 f 0 2 /ft; ft = k /c p /(4aT 2 ); f 0 = /c p r 0 ; ( = r/r 0 

The difference in definition for r 0 for parallel plate case and f Q for the circular tube case should 
be noted. For a black tube, the boundary conditions for Eqs. (26) and (27) are found to be 



~« qR) 


;qa(0) = 0 

<= 1 


(28) 


Equation (26) is used for general one-dimensional gray gas formulation and Eq. (27) for the 
case of linearized radiation. 
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IV. lYansformed Equations and Limiting Forms 

Depending upon the physical problem considered, a convenient form of the total radiative 
flux should be used in the governing energy equations. For nongray formulations, in order to 
be able to use the band model correlations, the equations for radiative flux must be expressed in 
terms of the correlation quantities. These equations are provided in this section. 

For a nongray gas consisting of n-molecular bands, a combination of Eqs. (5), (19) and (20) 
provides the governing equation for black parallel plate geometry as (17, 35, 40] 


d#/d£ — 2(3 £ 2 — 3£ 3 ) + 1 


= HiU « 


»«')*! 




df' 


i = l 


/V) a; 


2 U oi(^ - £) 


*} 


(29) 


where 

t = u / u o = y/L; ? = u'/uo = z/L; A = A/A 0 ; 
u = (S/A 0 )Py; u 0 = (S/A C )PL; PS = f AC^dw; 

J Aw 

Hi = A oi (T)(de^/dT) Tw 

and A'(u) denotes the derivative of A(u) with respect to u. Equation (29) is the most convenient 
equation to use when employing the band model correlations in nongray radiative transfer 
analyses. The limiting forms and their solutions, in the optically thin and large path length 
limits, are provided in [35, 40]. 


A combination of Eqs. (9), (20) and (25) provides the nongray form of the energy equation 
for the circular tube geometry as [23] 


69 

d? 


+ r - 2£ 




+ 


/V)A| 

r °(m 

sin y 


oi /^/ 


COS 7 

bu oi 
cos 7 


(e'-O 


dC 


(£ + ('- 2 ^ sin 


in 7 ) d^'j 


d7 


(30) 


where 


( = u/u 0 = r/r 0 ; u = (S/A 0 )Pr,u o = (S/A 0 )Pr 0 

The limiting forms of Eq. (30) are available in [23]. The large path length limit (large u oi limit) 
constitutes a very useful limit for analyzing the radiative transfer capability of gases like C0 2 
at room temperature and at one atmospheric pressure. For gases such as CO and C4, the results 
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obtained in this limit constitute an upper bound in radiative transfer capability. In the large path 
length limit, Eq. (30) reduces to 


5 + 


4r 0 

?rbk 


n /"V 2 ( r 1 

|cos 7 / 0(t') 

i = j Jo \ -/(sin7 


*£_ + _ 
K-e t+f 


je 11 

' -2^sin7j / 


d 7 


(31) 


Equations (30) and (31) are useful in determining the extent of radiative interaction of partic- 
ipating species in a circular tube. 

For linearized gray gas radiation, a combination of Eqs. (5) and (23) provides the equation 
for the radiative flux for the parallel plate geometry as 


d 2 qR/d£ 2 - ^ 7 1 + ^ r o^qR = 7iqw(6£ 2 - 4£ 3 - l) (32) 

The solution of Eq. (32) is obtained first, then the solution for 9 is obtained from Eq. (5), and 
finally the relation for the bulk temperature is obtained from Eq. (10) as [17] 


where 


9 h = Ci [24 - 12Mi + M 3 + (M 3 - 12Mi - 24)e“ M '] 
_ 12 7. 17^ _ 17 

5 Mf + 70 M? 70 


r> _ 71 

1 Mf 


48 — 3t 0 Mj + 36r 0 
3r 0 (l — e _Ml ) + 2Mj(l + e -Ml ) 


(33) 



Limiting forms (optically thin and optically thick) of Eqs. (33) are available in [17]. 


From a combination of Eqs. (9) and (27) the relation for linearized gray gas radiative flux 
for the circular tube is obtained as 


+ - 0 + M & 2 )qR = 72qw(2^ 3 - £ 5 ) (34) 

where 

M 2 = + 72 

Equation (34) is solved first for qR, then the solution for 9 is obtained from Eq. (9), and finally 
the bulk temperature result is obtained from Eq. (11) as [23] 


9 h = C 2 [M 2 (8 - M2)I 0 (M 2 ) - 16 I,(M 2 )] 
11 72 _ 8 j2_ _ n 

+ 24 3 M< 24 


( 35 ) 
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where 

^ _ 72 T 3f Q M2 - 24f 0 — 32 
2 “ M® [2M 2 I 0 (M 2 ) + 3f 0 I 1 (M 2 ). 

Limiting forms of Eq. (35) are available in the cited reference. 

No attempt is made to simplify the radiative flux equation for use in Eq. (12). In this study, 
however, the nonlinear form of Eq. (21) is used to obtain solutions for laminar compressible 
flow between two parallel black plates. 


V. Method of Solution 

The solution procedure for the fully-developed incompressible flow is completely different 
than the solution procedure for the laminar compressible flow. These procedures are discussed 
here briefly. 

A. Fully Developed Incompressible Flow 

For both the parallel plate and circular tube geometries, the general nongray solutions are 
obtained numerically by employing the method of variation of parameters. For both geometries, 
it was essential to assume quartic formulations for #(£) to obtain converged solution. In principle, 
the same numerical procedure applies to both the general and large path length limit cases. The 
entire solution procedure for the parallel plate geometry is provided in [38, 40] and for the 
circular tube in [41], 


B. Laminar Compressible Flow 

The governing equations, Eqs. ( 12 ) and (13), are transformed form the physical plane (x, 
y, t) to computational domain (£, 77 , t) to facilitate the treatment of general geometry. Equation 
( 12 ) is expressed in the computational domain as 


dV dF dG 

dt d£ drj 


= 0 


(36) 


where 

U = UJ; G = Gx*-Fy f ; 

F = Fyrj - Gx,; J - x f y, - y ( x. 

Equation (36) is solved by a time-asymptotic two step, explicit MacCormack method [42]. This 
method is second order accurate in space. If a solution to Eq. (36) is known at some time t = 
nAt, then the solution at next time step t = (n+l) At can be calculated from 


ur/ 1 = HAtjufj 

for each grid point (i, j). The operator L consists of predictor and corrector steps. For this study, 
the code developed by Kumar [43] was modified to include the radiation model. The details of 
the solution procedure are available in [ 44 ], 
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The radiative flux equations, Eq. (21), is discretized using the central differencing scheme 
as [39, 44] 


Ayf(l+ft) 


+ 


<U-i 


Ay?(l + Wj 


^yj(i + # 
qj+i = i-5«j 


/ 1 1 \ 9 2 


9j 


e j+l e j . e j e j-l 


P'Ayi 


Ayj 


(37) 


where 


Ayj = yj - yj— i ; ft = 


yj+i ~ yj 

yj - yj-i 


Equation (37) along with Eq. (22) forms a tridiagonal system of equation which can be solved 
efficiently by the Thomas algorithm. 


In the nongray gas formulation, the divergence of the radiative flux is evaluated using 
a central differencing scheme and is treated as radiative source term in the energy equation. 
Since the nongray formulation involves an integro-differential equation, the radiative flux term 
is uncoupled and treated separately [39]. 


VI. Results and Discussion 

Extensive results have been obtained fro fully developed incompressible laminar flows 
between two parallel plates and within a circular tube for participating species OH, CO, CH4, 
CO 2 and H 2 O. For the case of compressible flow between two parallel plates, however, results 
have been obtained only for pure H 2 O as participating species and different mixtures of H 2 O 
and air. Selected results for both cases are presented and discussed in this section. 

A. Fully- Developed Laminar Flows 

A complete discussion (and physical interpretations) of the various parameters entering into 
the present problem is given in [1, 10, 17, 23]. Numerical solutions were obtained in terms 
of nondimensional temperature and bulk temperature. Specific results were obtained for OH 
(2.8 n fundamental), 4.7p fundamental + 1st overtone bands), CO 2 (15/i, 4.3/z and 2.1 pi bands), 
H 2 O (rotational, 6.3 pi, 2.7 fi, 1.87/i and 1.38/i bands), and CH 4 (7.6/z and 3.3 pL bands) for which 
spectral information were obtained from [5, 9, 11, 16]. 

For CO, CO 2 , H 2 O, and CH 4 , the results obtained by employing the Tien and Lowder’s 
correlation for band absorptance are illustrated in Figs. 2-6 for the parallel plate geometry. The 
limiting value of Ob = -0.243 corresponds to negligible radiation, and the effect of radiation 
increases with increasing plate spacing. As would be expected, radiative transfer is more 
pronounced for higher pressures and wall temperatures. Also shown in Figs. 2-5 are the limiting 
solutions for large Uo (LLU). It is seen that, for a given wall temperature, the large u 0 limit can 
be obtained either by going to large values of L or to high pressures. These results also indicate 
that, for a particular wall temperature, the large u 0 limit for CO 2 is achieved at a relatively lower 
pressure than for other gases. As a matter of fact, for most practical purposes involving CO 2 
at room temperature, the results for the four gases is shown in Fig. 6 for a pressure of one 
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atmosphere and a wall temperature of 1,000 °K. As discussed in [1, 10, 17], the relative order 
of the four curves, for small values L is characteristic of the interaction parameter for optically 
thin radiation, and for large L is characteristic of the interaction parameter for the large u 0 limit. 

For the circular tube geometry, preliminary results have been obtained for different gases 
under varying conditions and selected results are presented in Figs. 7-12. The results for 
temperature variations within the tube are presented in Figs. 7 and 8 for different conditions. 
Since the temperature profiles are symmetric, the results are illustrated only for £ = 0 to £ = 0.5. 
The general band absorptance and large u 0 solutions for OH are illustrated in Fig. 7 for T w = 
500K and 1000K. It seems that the LLU limit is achieved for OH at T w = 500K and P = 1 atm. 
General results for diffemt species are compared in Fig. 8 for T 2 = 300K, P = 1 atm, and r 0 
= 5 cm. These results demonstrate the relative importance of radiative interactions in a circular 
tube. This is analogous to the case of parallel plate results presented in [40]. 

The bulk temperature results for the circular tube geometry are presented in Figs. 9-12 
for different gases. The results are expressed in normalized form as ^b(NR), where 0 b (NR) 
represents the case with no radiative interaction. The results, in general, exhibit the same 
trend as presented in Figs. 2-6 for the parallel plate geometry. However, the extent of radiative 
interactions is entirely different. This is because the circular geometry provides additional degrees 
of freedom for radiative interactions [23], This is the first effort where general nongray band 
absorptance results have been obtained for the circular tube geometry. 

B. Laminar Compressible Flow 

A computer code was developed to solve the two-dimensional Navier-Stokes equations for 
radiating subsonic laminar flows between two parallel black plates. The dimensions of the 
duct were taken as 3cm x 15cm. The radiative interaction was considered only in the normal 
direction. Extensive results have been obtained for pure H 2 0 as participating species and different 
mixtures of H 2 0 and air flowing laminarly between the plates. These are available in [39, 44] 
and selected results are presented in this section. 

The results for subsonic flows were obtained for two specific Mach numbers, = 0.3 
and Mqo = 0.8. Most results presented here are for = 0.3; however, certain results for 
Moo = 0.8 are also presented for comparative purposes. For Mqo = 0.3, the free stream (inflow) 
conditions considered are, P, = 1.064 P^, T t = 1.018 T*,, P^ = 1 atm. Too = 500 K, f Hj0 
= 0.5, and f«r = 0.5. The freestream conditions corresponding to Moo = 0.8 are; P, = 1.524 
Poo* T t = 1.128 Too, Poo = 1 atm. Too = 500K, f^o - 0.5, and f„} r = 0.5 For both flows, 
the wall temperature was maintained at T w = 1,500K and the pressure at the channel exit was 
taken to be one atmosphere. 

The variation in axial velocity across the channel is shown in Fig. 13a for different x 
locations. It is clearly evident that the fluid velocity in the inviscid core increases along the axial 
length due to increase in the boundary layer. A fully-developed flow has not been achieved at 
the channel exit. Figure 13b shows the variation of normal component of velocity across the 
channel at various axial locations. The two-dimensional effects are clearly evident from the 
results x = 0.015 m which is predominantly in the entrance region. The magnitude of v velocity 
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decreases as x increases; it is positive in the lower half and negative in the upper half of the 
channel due to the symmetry of the problem. For a fully-developed flow, the v component of 
velocity should be zero. 

The variation in temperature across the channel is shown in Fig. 14 for different x locations. 
The temperature decreases from the walls to the center of the channel indicating the development 
of the thermal boundary layer. The flow is not thermally developed at the exit of the channel. 

The variations in conductive and radiative fluxes along the length of the plate are shown in 
Fig. 15 for different y locations. The fluxes are zero along the centerline of the plate because 
of physical symmetry. The conductive flux increases slowly with increasing x and becomes 
essentially constant at the channel exit. The radiative interaction is strong in the first ten percent 
of the length and then it slowly decreases and reaches about the same value as conductive flux 
at the exit. 

The variations in qc y , qRy and qty with y are shown in Figs. 16, 17 and 18, respectively. It 
is noted that the variation in conductive flux across mainly with the first 20% distance from the 
wall (Fig. 16). The situation, however, is not the same with respect to the radiative transfer (Fig. 
17). The radiative interaction occurs in the entire width of the channel. The results presented 
in Fig. 18 demonstrate that the rate of total energy transfer in the y direction is highest closer 
to the entrance (x = 0.015 m) than any other x locations. This is a result of high enthalpy flow 
and strong radiative interaction in this region. 

Figure 19 shows the variation of radiative and conductive flux across the channel. The 
radiative flux peaks at some distance away from the wall while conductive flux peaks at the 
wall. This is because at the lower wall the effect of positive heat flux is partially cancelled by 
the negative flux from the layers of hot gases next to the wall. At a small distance away from 
the wall, however, the positive flux from both the wall and gas combine to give a maximum 
heat flux. Farther into the gas, the flux from the wall and hot gas is attenuated bya the cooler 
gas and is partially cancelled. At the center of the channel, the fluxes from both modes cancel 
each other. This effect prevails throughout the length of the channel as is clearly evident from 
the results of Figs. 19a and 19b and in Ref. 44. 

The variation of the radiative, conductive and normal component of the total flux across the 
channel is shown in Fig. 20 for x = 0.06 m and 0.15m. The total flux includes the convective 
(flow), radiative and conductive fluxes. It is increasing to note that except very near the boundary, 
the total flux is about an order of magnitude higher than the radiative flux and much higher than 
the conductive flux. In other words, according to the definition used in this study, the convective 
(flow) energy is the predominate mode of energy transfer. This trend is result continues for 
all axial stations [44]. 

The results for M m = 0.8, in general, show the same trend as the results for = 0.3. 
Certain definite changes, however, are noted specific cases. For = 0.8, the magnitude of 
velocity variations obviously are relatively higher but the velocity boundary layers are thinner 
[44], The temperature variations shown in Fig. 21 for = 0.3. This is because Too and T w 
are the same for both flows. The thermal boundary layer, however, is thicker for = 0.3. 
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A comparison of = 0.8 and 0.3 results for qc y , qRy and qty is shown in Fig. 22 for x 
= 0.12 m. The results demonstrate that while conductive and total energy fluxes are influenced 
considerably by changes in the Mach number, the radiative energy transfer is insensitive to such 
changes. This is because radiation is influenced primarily by specified temperature and pressure 
conditions and these conditions are the same for both cases. 

VII. Concluding Remarks 

Analytical formulations and numerical procedures have been developed to investigate the 
radiative interaction of absorbing-emitting species in laminar fully-developed flows between 
parallel plates and within a circular tube. The general nongray results for the circular tube 
geometry have been obtained for the first time. The results demonstrate the relative ability of 
various participating species for radiative interactions. 

Two-dimensional compressible Navier-Stokes equations have been used to investigate the 
influence of radiative energy transfer on the entrance region flow under subsonic flow conditions. 
Computational procedures have been developed to incorporate gray as well as nongray formu- 
lations for radiative flux in the general governing equations. Specific results have been obtained 
for different amount of water-vapor in water vapor-air mixtures. Results demonstrate that the 
radiative interaction increases with an increase in pressure, temperature and the amount of water 
vapor. This can have a significant influence on the overall energy transfer in the system. Most 
energy, however, is transferred by convection in the flow direction. As a result, the radiative 
interaction does not alter the flow field significantly. Further parametric studies are needed to 
make definite recommendations. 
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INVESTIGATION OF RADIATIVE INTERACTIONS 
IN HIGH-SPEED ENTRANCE REGION FLOWS 


S. N. Tiwari and D. J. Singh 
Abstract 

The influence of radiative energy transfer on the entrance region flow is investigated under 
supersonic flow conditions in a channel. Two-dimensional compressible Navier-Stokes equations 
are solve numerically in conjunction with the radiative flux equations. The channel walls are 
assumed to be black. Nongray as well as pseudo gray gas models are used to represent the 
absorption-emission characteristic of the medium. The participating species considered are 
different amounts of water vapor in water vapor-air mixtures. Results obtained for different 
flow conditions indicate that the radiative interaction can moderately influence the overall energy 
transfer, but the flowfield is not changed significantly. 

1. INTRODUCTION 

The analysis of combined modes of heat transfer involving conduction, convection and 
radiation has been the subject of considerable research in the last several years. The problem of 
combined modes of energy transfer is a very complicated phenomenon. While conduction and 
convection processes can be described by differential equations, radiation is described by integral 
equation. The integrations are with respect to frequency, solid angle, and spatial coordinates. At 
present, there is no analytical solution available for governing equations describing the combined 
modes of energy transfer, even for simple geometries. However, with the availability of fast 
computers, it is possible to solve these complicated equations numerically. 

In many engineering problems involving high temperature gases such as reentry, hypersonic 
propulsion, design of combustion chambers for high pressure spacecraft engines, chemical 
transfer vehicles, heavy lift launch vehicles, and furnaces, the radiative interaction becomes very 
important. Under certain conditions, the radiation is the predominant mode of heat transfer. 
Basic formulations on radiative transfer in participating mediums are available in standard 
references [1-8]. The review articles presented in [9-15] are useful in understanding the radiative 
properties of participating species and the nature of nongray radiation. In most of the studies 
involving radiation, certain simplifying assumptions such as gray gas [16-18], fully developed 
flow [19-22], constant thermophysical properties [23-25] and linearized radiation [10, 20-22, 
26-30] are invoked. The validity of radiative transfer analysis depends upon the accuracy with 
which absorption-emission characteristics of participating species are modeled. There are several 
models available in literature; and these are reviewed critically in [12, 13]. 

In general, there has been a lack of studies on radiative interaction in the entrance region of 
ducts. In this region, momentum and energy equations cannot be decoupled. Some simplified 
studies on this problem are available in [5, 31-33] and the cited references. Most of these 
studies are restricted to constant properties and incompressible flows. Im and Ahluwalia [34] 
analyzed the combined conduction, convection and radiation in rectangular ducts by solving 
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simultaneously the fluid dynamical and the radiation transport equations. The flow was assumed 
as parabolic thereby excluding the possibility of flow recirculation. Soufiani and Taine [35] 
studied the thermal and dynamical entrance region for a steady laminar flow between two parallel 
plates. The walls were considered diffuse and fluid as homogeneous absorbing, emitting and 
nonscattering medium. The viscous dissipation effects and axial diffusion were neglected in 
the energy equation. The results showed a significant difference in temperature and velocity 
profiles by inclusion of radiation. Mani et al. [36] studied the chemically reacting and radiating 
supersonic flow between parallel plates in a channel with a ten degree ramp. They showed that 
in the case of flow without chemical reaction, most of the energy is transferred by convection. 
As a result, the flow field is not affected significandy. 

The objective of this study is to investigate the influence of radiative interactions on the 
subsonic and supersonic region flows. To accomplish this, a mixture of water vapor an air 
flowing laminarly between parallel plates is considered. The walls are maintained at constant 
temperatures. The strong temperature gradient between walls and the fluid requires the use of 
temperature dependent fluid properties. The flow and energy conservation equations, therefore, 
are solved simultaneously. The effects of various parameters such as Mach number, pressure, 
amount of participating species, and the plate separation distance are investigated. 

2. BASIC THEORETICAL FORMULATION 


2.1 Flow Equations 

Two-dimensional Navier-Stokes equations in fully conservative form are used to describe 
the flow field. These equations, in physical domain, can be written as (Fig. 1) 


dU_ 5G 

dt + dx + dy 


where vectors U , F and G are expressed as 


( 2 . 1 ) 


U = [p pu pH — P] 

F = [p pu 2 +p- T xx puv - T xy (Et + p)u - UT XX - VT X y + q Ct + qRx ] 

G = [p PUV - T xy pv 2 +P~ Tyy (E t -I- p)v - UT X y ~ VTyy + + tfflj,] 

The viscous stress terms are given by 
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where qc x and qn y are component of conductive flux in x and y directions, respectively, and 
are expressed as 


dT 

qcx = -A — 

OX 

(2.3a) 

dT 

qcy -~ K dy 

(2.3b) 


The terms qp x and qn y in the definitions of vector F and G represent the radiative flux in x and 
y directions, respectively. The relations for the radiative flux are developed in Sec. 3. The total 
energy flux in a given direction is given by the corresponding last term in the definition of F and 
G. The coefficient of viscosity p is evaluated using the Sutherland’s formula and the coefficient 
of thermal conductivity is calculated by using a constant value of the Prandtl number equal to 
0.72. In order to complete the system, the equation of state is used as 


P = pRT 


(2.4) 


The above equations can be used for subsonic as well as supersonic flows. However, the 
boundary conditions are numerical procedures for the two flows are quite different. 

2.2 Inlet and Boundary Conditions 

The treatment of the inflow boundary conditions is guided by the theory of characteristic. 
A locally one-dimensional flow has four characteristic equations with slopes u, u+c, u and u-c. 
If the flow is subsonic at the inflow, then the u-c characteristic has a negative slope and it 
propagates information from the interior upstream to the inflow boundary. In this case, only 
three quantities can be specified at the inflow and the fourth quantity must be allowed to vary 
as the solution progresses. 

In this study, the stagnation pressure, stagnation temperature and flow angle are specified 
at the inflow. These quantities are related to the state pressure and state temperature by the 
following equations 

M i+ 1 f L m2 ) (2 - 6) 

— = tan ( &ip ) (2.7) 

u 

Equations (2.5)-(2.7) are a system of three equations in four unknown, P, T, u and v. To 
complete the system of equations, a zero order extrapolation is used for the pressure at the 
inflow. The outflow boundary is also calculated based on theory of characteristic. For subsonic 
flow at the outflow, the u-c characteristic propagates information upstream from the boundary to 
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the interior, i.e. only one quantity can be specified at the outflow. The state pressure is specified 
at the outflow while u and v and T are calculated using a zeroth-order extrapolation. Along the 
surfaces, following boundary conditions were applied 

« = 0, u = 0, d ? = 0, T = T W (2.8) 

dy 

The density is obtained from the equation of state using the computed surface pressure and 
prescribed surface temperature. Fourier’s heat conduction law is used to calculate the conductive 
heat flux. 

For supersonic flows, the inflow properties are assumed to be chemically frozen. The 
outflow conditions are calculated from adjacent nodes using a first order extrapolation. The 
pressure at walls are calculated by a second order one-sided differencing procedure. The free 
stream conditions considered for each case are indicated in discussing the results. 

3. RADIATIVE TRANSFER MODELS 

Radiative transport is a quite complicated phenomenon, as an element not only exchanges 
energy with its neighbors, but also from all other elements. This results in an integral expression 
with integration with respect to frequency, solid angle, and spatial coordinates. The applicability 
and validity of the radiative transfer in a physical problem depends on the accuracy with which 
the radiation is modeled. In most studies, the “tangent slab approximation” [8, 10, 36] is invoked; 
it treats the gas layer as one-dimensional slab in evaluation of the radiative flux (Fig. 1). In recent 
years, there have been few studies [36-39] which include the multidimensional effects. In this 
section, a brief discussion of various absorption models is given and essential equations for the 
radiative flux are presented. 

3.1 Absorption Models 

The gray gas approximation is the simplest model to calculate the radiative flux. Although 
unrealistic, this model has yielded accurate prediction of heat flux burning rates in moderate size 
fires [40]. It has also been useful for estimating the total radiative flux in furnaces [41]. In this 
model, the absorption coefficient is assumed to be independent of frequency. It overestimates the 
radiative flux because participating gas is actually transparent over large region of spectrum, and 
may be practically opaque in other regions. Thus, an attempt to represent the entire spectrum 
with a mean absorption coefficient will lead to overestimation of radiant flux in the transparent 
and opaque regions. The model, however, is very useful as first approximation in providing 
insight into the parametric trends to be expected in the system and in estimating the significance 
of radiation. 

There are various nongray models available which relax the gray gas assumption. A few of 
these are line by line models, narrow band models and wide band models. These are discussed 
here briefly. 

The line by line models [12, 13] by far are the most accurate models available. They 
account for each line of the absorbing species. However, these models are very expensive to 
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use. Since most applications require knowledge of only the properties averaged over spectral 
range, the narrow band models which are computationally fast and accurate are widely used. Four 
commonly used narrow band models are the elasser, statistical, random elasser and quasi random 
[12, 13]. Next in order are the wide band models which provide correlations that are valid over 
the entire band pass. These correlations are quite accurate for many engineering applications. 
The most commonly used wide-band model correlations are the Edwards exponential wide band, 
Tien and Lowder, and Cess and Tiwari [9-13]. The application of a model to a particular case 
depends on the nature of the absorbing emitting species. 

3.2 Radiative Flux 

For a nonscattering medium and diffuse boundaries, the expression for the radiative flux is 
given as [1, 42] 


qRu>{y ) = 2e iw E 3 (K w y) - 2e 2a ,£ 3 [/c w (L - y)] 
+ 2 {/ e ^(0«w^2[«a;(y - OR 

L 

- j tw{0 K uE2[Kw{( - y)R| 


(3.1) 


where 


and 


Eiiv) = 2 ex P 


EM 


(-?) 
= r xp ( J i) 


The total radiative flux is given by 


fOO 

QR — I QRuj du) 

Jo 


(3.2) 


For a gray medium, k w is independent of the wave number and an expression for the radiative 
flux is obtained from Eqs. (3.1) and (3.2) as 


qR = (y) = e! - e 2 + 

’■L 


w: 


[e(0 - e^exp |--/c(y -0 


K,d£ 


~ J [e(0- e 2]exp -^(y-0 j 


(3.3) 


Differentiating Eq. (3.3) twice, the integrals are eliminated and there is obtained a nonhomoge- 
nous ordinary differential equation as [1, 36] 


1 d 2 qR(y) 
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[2 dy 2 -!««(») = r 


3 de(y) 
k dy 


(3.4) 
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This is a second order differential equation; hence it requires two boundary conditions. For 
nonblack diffuse surfaces, these are given as 

In order to solve Eq. (3.4), one needs an expression f 

r e " ,( ’ 

Here, k represents the Planck mean absorption coefficient and is function of temperature and 
species partial pressure P r 

An expression for the divergence of radiative flux in the optically thin limit is obtained from 
Eq. (3.1) as [1]. 


“HR I 

[dy\ 

[dy 


= 0 


y=o 


= 0 


(3.5a) 

(3.5b) 


I y—L 


i>r «, which is given as 


r )Si(T) 
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dqR„ 

dy 


= 2/c w [2e ftu ,(i/) - ei - e 2 ] 


In the optically thick limit, Eq. (3.1) reduces to [1] 


(3.7) 




4 defo 

3«u; dy 


(3.8) 


Equations (3.7) and (3.8) can be used to obtain limiting solutions for the gray model. 

For a nongray model consisting of n-molecular (vibration-rotation) bands, a combination of 
Eqs. (3.1) and (3.2) results in [1, 42] 


<7rt(0 = ei - e 2 
+ 


l E j F t „ - o] 


where 

^iwi(0 = «<*(£) - = e Wi (£) - e 2a)i .,£ = y!L = u/u a 

and A'(u) denotes the derivative of A(u) with respect to u. Equation (3.9) is the most convenient 
equation to use when employing the band model correlations in nongray radiative transfer 
analyses. Spectral properties and correlation quantities for various participating species are 
available in [9, 11], 
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In the optically thin limit [1, 10], A(u ) = u and A'(u) = 1. In this limit, the divergence of 
radiative flux is obtained by differencing Eq. (3.9) as 


< 3 “» 

y i=i 

In the large path length (i.e., for u 0 , >> 1 for each band), /4'(u) = ln(u) and A'(u) = 1/u. 
Consequently, in this limit, one obtains a simple relation from Eq. (3.9) as [1, 10] 


Mr _ y- . de w , M' 

d( ~ °'Jo d? 


(3.11) 


Equations (3.10) and (3.11) are used to obtain limiting solutions for nongray radiation. 

For nongray radiative interactions, the continuous correlation proposed by Tien and Lowder 
[9] is employed in this study. This correlation is relatively simple and provides accurate results 
for pressures higher than 0.5 atmosphere. 


4. METHOD OF SOLUTION 

The governing equations, Eqs. (2.1) through (2.3), are transformed from the physical plane 
(x, y, t) to computational domain (£, tj, t) to facilitate the treatment of general geometry. 
Equations (2.1) and (2.2) are expressed in the computational domain as 


where 


dU_ dF dG_ (] 

dt + di + drj ~ 0 


U = UJ 


(4.1) 


F = Fy v — Gx n \ G = Gx ^ = Fy { 

J = x^y v - y^x v 

Equation (4.1) is solved by a time-asymptotic two steps, explicit MacCormack method [43]. 
This method is second order accurate in space. If a solution to Eq. (4.1) is known at some time 
t = nAt, then the solution at next time step t = (n + 1)A< can be calculated from 

J )”/ 1 = 

for each grid point (*, j). The operate L consists of predicter and corrector steps. For this study, 
the code developed by Kumar [44] was modified to include the radiation model. The details of 
the solution procedure are available in [44, 45]. 
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The radiative flux equations, Eq. (3.4), is discretized using the central differencing scheme 
as 


AyjO + Pj) 


<7j-l 




L At/;(1 + 

q J+ i = 1.5 k j 


( 1 1 \ 9 j 
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e i+l ~ e j , e j ~ e j -l 


0j A y } 


A y } 


(4.2) 


where 

A y } = y> - y>-i 


p. _ y>+i - 
J V] - yy-i 

Equation (4.2) along with Eq. (3.5) forms a tridiagonal system of equation, which can be 
efficiently solved by the Thomas algorithm. 

In the nongray gas formulation, the divergence of the radiative flux is evaluated using 
a central differencing scheme and is treated as radiative source term in the energy equation. 
Since the nongray formulation involves an integro-differential equation, the radiative flux term 
is uncoupled and treated separately [45]. 


5. RESULTS AND DISCUSSION 

Based on the theory and computational procedure described in the previous sections, a 
computer code was developed to solve the two-dimensional Navier-Stokes equations for radiating 
supersonic laminar flows between two parallel black plates. A similar but different code was 
developed for radiating subsonic flows. The dimensions of the channel were taken as 3 cmx 
10 cm for the supersonic flow case and 3 cmx 15 cm for the subsonic case. The radiative 
interaction was considered only in the normal direction. Extensive results have been obtained 
for pure H 2 O as participating species and different mixture of H 2 O and air flowing laminarly 
between the plates. These are available in [36, 45] for the supersonic case and in [46] for the 
subsonic case. Selected results for both case are presented and discussed in this section; results 
for supersonic flows are presented first. 

For the case of supersonic flow, a comparison of the divergence of radiative flux for general 
(nongray), gray, and their optically thin limit models is presented in Fig. 2 for two different 
y-locations (y=0.2 and 1.5 cm). The inflow conditions for this case are Poo = 1 atm, Too = 
1,700 K, Moo = 4.3, f^o - 0.5, f 0j =0.1, and = 0.4. The gray formulation is based on 
the Planck mean absorption coefficient which accounts for the detailed information on different 
molecular bands. As such, this approach is referred to as the “pseudo-gray formulation.” The 
range of optical thickness calculated in [45] was found to be between 0.003 and 0.4. Thus, for 
the physical model and inflow conditions considered, the radiative interaction is essentially in 
the optically thin range. No significant difference in results is observed for the two y-locations; 
this is a typical characteristic of the optically thin radiation [1, 10]. The solution of the gray 
formulation requires about ten times less computational resources in comparison to the solution 


9 


of the nongray formulation [45]. As such, for basic parametric study, all other results presented 
in this section were obtained by using the pseudo-gray gas formulation. 

Another case considered for the supersonic flow corresponds to the inflow conditions of Poo 
= 1 atm, = 1,700 K, Uoo = 2574 m/s (Moo = 3.0), and various amounts of water vapor and 
air mixtures. The results for radiative flux are illustrated in Figs. 3 and 4 as a function of the 
nondimensional y-coordinate. For P = 1 atm, the results presented in Fig. 3 for different water 
vapor concentrations indicate that the radiative interaction increases slowly with an increase 
in the amount of the gas. The results for 50% H 2 O are illustrated in Fig. 4 for two different 
pressures (Poo = 1 and 3 atm) and x-locations (x = 5 and 10 cm). It is noted that the increase 
in pressure has dramatic effects on the radiative interaction. The conduction and radiation heat 
transfer results are compared in Fig. 5 for P = 3 atm and for two different x-locations (x = 
5 and 10 cm). The results demonstrate that the conduction heat transfer is restricted to the 
region near the boundaries and does not change significandy from one x-location to another. 
The radiative interaction, however, is seen to be important everywhere in the channel, and this 
can have significant influence on the entire flowfield. The results presented in Figs. 3-5 should 
be physically symmetric; but, due to the predictor-corrector procedure used in the McCormack’s 
scheme, they exhibit some unsymmetrical behavior. 

The results for subsonic flows were obtained for two specific Mach numbers, Moo = 0.8 and 
Moo = 0.3. Most results presented here are for Moo = 0.8; however, certain results for Moo — 
0.3 are also presented for comparative purposes. 

The free stream (inflow) conditions corresponding to Moo = 0.8 are, Pt = 1.524 Poo, Tt = 
1.128 = 1 atm. Too, Poo = 1 atm, Too = 500 K, f Hj o = 0.5, and f* = 0.5. The wall temperature 
was maintained at T w = 1 ,500 K. The pressure at the channel exit was taken to be one atmosphere. 

The variation in axial velocity across the channel is shown in Fig. 6a for different x locations. 
It is clearly evident that the fluid velocity in the inviscid core increases along the axial length 
due to increase in the boundary layer. A fully-developed flow has not been achieved at the 
channel exit. Figure 6b shows the variation of normal component of velocity across the channel 
at various axial locations. The two-dimensional effects are clearly evident from the results of x = 
0.015 m which is predominantly in the entrance region. The magnitude of v velocity decreases 
as x increases; it is positive in the lower half and negative in the upper half of the channel 
due to the symmetry of the problem. For a fully-developed flow, the v component of velocity 
should be zero. 

The variation in temperature across the channel is shown in Fig. 7 for different x locations. 
The temperature decrease from the walls to the center of the channel indicating the development 
of the thermal boundary layer. The flow is not thermally developed at the exit of the channel. 

The variations in conductive and radiative fluxes along the length of the plate are shown 
in Fig. 8 for different y-locations. The fluxes are zero along the centerline of the plate because 
of physical symmetry. The conductive flux increases slowly with increasing x and becomes 
essentially constant at the channel exit. The radiative interaction is strong in the first ten percent 
of the length and then it slowly decreases and reaches about the same value as conductive flux 
at the exit. The variation in total energy flux (qty) with x is illustrated in Fig. 9 for different 
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y-locations. A comparison of results presented in Figs. 8 and 9 demonstrates that the rate of 
convective (flow) energy is considerably higher than the rate of conductive or radiative energy 
transfer. 

The variations in qcy, qRy and qty with y are shown in Figs. 10, 11 and 12, respectively. 
It is note that the variation in conductive flux occurs mainly within the first 20% distance from 
the wall (Fig. 10). The situation, however, is not the same with respect to the radiative transfer 
(Fig. 11). The radiative interaction occurs in the entire width of the channel. This is similar to 
the results for supersonic flow presented in Fig. 5, although the extent of interaction is entirely 
different. The results presented in Fig. 12 demonstrate that the rate of total energy transfer in 
the y-direction is highest closer to the entrance (x = 0.015 m) than any other x location. This is 
a result of high enthalpy flow and strong radiative interaction in this region. 

Figure 13 shows the variation of radiative and conductive flux across the channel. The 
radiative flux peaks at some distance away from the wall while conductive flux peaks at the 
wall. This is because at the lower wall the effect of positive heat flux is partially cancelled by 
the negative flux from the layers of hot gases next to the wall. At a small distance away from 
the wall, however, the positive flux from both the wall and gas combine to give a maximum 
heat flux. Farther into the gas, the flux from the wall and hot gas is attenuated by the cooler 
gas and is partially cancelled. At the center of the channel, the fluxes from both modes cancel 
each other. This effect prevails throughout the length of the channel as is clearly evident from 
the results of Figs. 13a and 13b and Ref. 46. This trend in result is somewhat similar to the 
trend for the supersonic flow presented in Fig. 5. 

The variation of the radiative, conductive and normal component of the total flux across 
the channel is shown in Fig. 14 for x = 0.06 m. The total flux includes the convective (flow), 
radiative and conductive fluxes. It is interesting to note that except very near the boundary, the 
total flux is about an order of magnitude higher than the radiative flux and much higher than the 
conductive flux. In other words, according to the definition used in this study, the convective 
(flow) energy is the predominant mode of energy transfer. This trend in result continues for 
all axial stations [46]. 

For Moo = 0.3, the free stream conditions considered arc, P t = 1.064 Poo, T t = 1.018 Too, 
Poo = 1 atm. Too = 500 K, fH,o = 0.5, and f^ = 0.5. In this case also, the wall temperature 
was maintained at T w = 1,500 K and the channel exit pressure was one atmosphere. 

The results for Moo = 0.3, in general, show the same trend as the results for Moo = 0.8. 
Certain definite changes, however, are noted for specific cases. For Moo = 0.3, the magnitude 
of velocity variations obviously are relatively lower but the velocity boundary layers are thicker 
[46]. The temperature variations shown in Fig. 15 for Moo = 0.3 are about the same as shown 
in Fig. 7 for Moo = 0.8. This is because Too and T w are the same for both flows. The thermal 
boundary layer, however, is thicker for Moo 0.3. 

A comparison of Moo = 0.8 and 0.3 results for qc y , qRy and qty is shown in Fig. 16 for x 
= 0.12m. The results demonstrate that while conductive and total energy fluxes are influenced 
considerably by changes in the Mach number, the radiative energy transfer is insensitive to such 
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changes. This is because radiation is influenced primarily by specified temperature and pressure 
conditions and these conditions are same for both cases. 

6. CONCLUSIONS 

Two-dimensional compressible Navier-Stokes equations have been used to investigate the 
influence of radiative energy transfer on the entrance region flow under supersonic and subsonic 
flow conditions. Computational procedures have been developed to incorporate gray as well 
as nongray formulations for radiative flux in the general governing equations. Specific results 
have been obtained for different amounts of water-vapor in water vapor-air mixtures. Results 
demonstrate that the radiative interaction increases with an increase in pressure, temperature and 
the amount of water vapor. This can have a significant influence on the overall energy transfer 
in the system. Most energy, however, is transferred by convection in the flow direction. As 
a result, the radiative interaction does not alter the flow field significantly. Further parametric 
studies are needed to make definite recommendations. 
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NOMENCLATURE 


A band absorptance = A(u, /?), cm' 1 

A 0 band width parameter, cm' 1 

Co correlation parameter, atm' 1 - cm' 1 

C p specific heat at constant pressure kJ/kg-K = erg/gm-K 

eu > Planck’s function, (W-cm' 2 )/cm' 1 

e Wt> Planck’s function evaluated at wave number u> 0 

2 

ei,e 2 emissive power of surfaces with temperature T| and T 2 , W-cm' 

k thermal conductivity, erg/cm-sec-K 

P pressure 

Pi partial pressure 

P t stagnation pressure 

Pr Prandtl number 

qR total radiative heat flux, J/m 2 -s 

qc conduction heat flux, J/m 2 -s 
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9Ru> 

spectral radiation heat flux, (w-cm'^/cm 1 

S 

integrated intensity of a wide band, atm" 1 -cm 

T 

temperature, K 

t,,t 2 

wall temperature, K; Ti = T w 

u 

nondimensional coordinate = SPy/Ao 

Uo 

nondimensional path length = SPL/Ao 

e 

flow angle 

K u> 

spectral absorption coefficients, cm' 1 


computational coordinate 

P 

density, kg/m 

a 

Stefan-Boltzmann constant, erg/(sec-cm 2 -K 4 ) 

uj 

wave number, cm' 1 

UJo 

wave number at the band center, cm' 1 
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Radiative flux vs. y at the channel exit, = 3.0. 
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Fig. 8 Variation of q CT 
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Fig. 12 Variation of 
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Fig. 16 Temperature variation across the channel at various x - stations. 









